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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC) , under 
program  element  65807F,  The  results  were  obtained  by  ARO,  Inc.,  AEDC 
Division  Ca  Sverdrup  Corporation  Company) , operating  contractor  for 
the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  research  was 
done  under  ARO  Project  No.  V33A-A8A.  Elton  R.  Thompson  was  the  Air 
Force  project  manager.  The  manuscript  was  submitted  for  publication 
on  September  2,  1977. 
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1.0  INTRODUCTION 

This  document  describes  a computer  program  which  allows  greatly 
simplified  Input,  for  certain  common  axlsymmetrlc  and  nonaxlsymmetric 
nose-vehicle  geometries  of  general  interest,  to  the  more  general  Grumman 
QUICK-geometry  system.  QUICK,  a Fortran  computer  program  developed  by 
Vachrls  and  Yaeger  (Refs.  1 and  2),  provides  the  user  with  a simple  and 
well  organized  method  of  constructing  and  Interrogating  a mathematical 
model  of  an  arbitrary,  three-dimensional  flight  vehicle.  This  geometry 
system  has  primary  application  to  the  numerical  computation  of  the  flow 
fields  over  arbitrary  aerodynamic  configurations.  QUICK  was  developed 
in  two  phases,  the  QUICKDEF  phase,  which  provides  facilities  to  build 
and  check  the  geometry  model,  and  the  QUICKLOK  phase,  which  provides 
surface  coordinates  and  derivatives  to  user  programs  without  requiring 
of  the  user  a detailed  knowledge  of  how  the  model  was  constructed. 
Current  applications  of  the  QUICKLOK  phase  Include  the  time- dependent 
blunt  body  program  BLUNT  (Ref.  3);  the  supersonic,  three-dimensional, 
external,  inviscid  flow-field  code  STEIN  (Refs.  3 and  4);  and  the 
ffiultishocked,  three-dimensional,  supersonic  finite  difference  method  of 
Kutler,  Reinhardt,  and  Warming  (Ref.  5).  The  QUICKDEF  phase,  of  primary 
concern  here,  is  a highly  user-oriented  program  which  greatly  simplifies 
the  geometric  input  by  providing  general  curve  fitting  services  to 
define  arbitrary  bodies  from  a minimum  of  logical  and  numerical  input. 

In  essence,  QUICK  requires  the  user  to  define  each  unique  cross  section 
of  the  vehicle  in  a plane  normal  to  the  coordinate  axis  of  the  vehicle 
by  merely  identifying  the  types  of  curves  (lines,  ellipses,  parabolas, 
...)  which  form  the  segmented  cross  section.  The  cross  sections  then 
become  numerically  specific  when  the  user  further  defines,  logically  and 
numerically,  certain  key  control  points  which  determine  the  intersec- 
tions and  slopes  of  the  curves  comprising  the  cross  sections.  The  only 
two  limitations  Imposed  by  QUICK  upon  the  vehicle  geometry  are  (1)  that 
the  vehicle  must  have  a vertical  plane  of  symmetry  and  (2)  that  each 
cross  section  must  be  definable  in  terras  of  a single-valued  radius 
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vector  swung  about  a single  point  within  the  cross  section.  While  these 
limitations  prevent  complete  treatment  of  certain  vehicle  geometries 
(such  as  engine  inlets,  the  oblique  wing  supersonic  transport,  or  the 
double  vertical  stabilizer  of  the  F-15),  QUICK  is  an  important  advance 
in  the  description  of  vehicle  geometry  and  makes  practical  the  solution 
of  flows  over  such  complicated  bodies  as  the  space  shuttle  orblter. 

In  contrast,  QUICK  is  also  well  suited  to  the  treatment  of  rela- 
tively simple  bodies  such  as  sphere-cones,  biconics,  and  slab  deltas, 
which  are  of  basic  interest  to  aerodynamlcists.  It  was  therefore  deemed 
logical  to  develop  yet  another  sequential  computer  program  to  set  up 
input  to  QUICK  for  certain  such  simple  geometries,  thus  reducing  user- 
prepared  input  to  a matter  of  two  or  three  cards  for  the  programmed  ' ' 
geometries.  This  fourth  program  FREKWIK  is  documented  in  a previous 
user's  manual. 

During  the  development  of  PREKWIK,  the  need  was  recognized  to' study 
more  complicated  nose  shapes  than  the  simple  spherical  cap  treated  by 
PREKWIK.  Interest  centered  around  the  effects  of  ablated  nose  shapes 
(Fig.  1)  upon  the  aerodynamics  of  otherwise  simple  bodies  such  as  cones 
and  biconics  with  slices  on  the  conic  portions  of  the  vehicle.  A feature 
which  most  of  the  symmetrical  noses  have  in  common  is  that  they  can  all 

be  modeled  (with  varying  degrees  of  accuracy)  as  a stack  of  conical  and 

1 

ogival  sections  which  intersect  in  different  fashions.  Thus,  it  appeared 
practical  to  program  the  logic  needed  to  model  these  geometries  and  set 
up  the  input  to  QUICK.  Initially,  the  nose  code  was  envisioned  as  a ' 
modification  to  PREKWIK,  but  the  need  for  a fundamentally  different  ' 
programming  approach  soon  forced  the  nose  geometry  problem  into  a 
separate  computer  program.  The  essential  goal  of  PREKWIK  was  to  receive 
a minimum  of  geometric  input  data  via  very  few  input  options  for  a given 
geometry,  perform  a few  relatively  trivial  intersection  computations, 
and  then  generate  the  card  image  input  required  by  QUICK.  For  the 
programmed  geometries,  PREKWIK  eliminated  the  tedious  chore  of  preparing 
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a new  QUICK  deck  for  every  minor  quantitative  change  in  geometry  (such 
as  a different  bluntness,  cone  half-angle,  or  vehicle  length).  In 
contrast,  the  nose  geometry  problem  was  perceived  as  logically  more  com- 
plicated, Not  only  would  the  nose  program  have  to  compute  the  inter- 
sections of  more  complicated  functions,  but  the  capability  to  round  or 
fillet  sharp  corners  was  deemed  essential.  However,  the  analytic  geo- 
metry of  rounds  and  fillets  - determlng  the  equation  and  intersection 
points  of  a circle  tangent  to  two  Intersecting  curves  - can  result  In 
very  high  degree  polynomials  along  with  the  attendant  logic  problems  of 
root  sorting  to  choose  the  correct  intersection.  The  computational 
complexity  of  the  Intersection  problem  coupled  with  the  fairly  large 
number  of  ways  a given  geometry  can  be  specified  via  various  combina- 
tions of  known  parameters  forces  the  geometry  program  to  accept  many  dif- 
ferent input  combinations  from  which  must  be  computed  the  same  basic  set 
of  QUICK  input  coordinates  and  slopes.  If  the  nose  code  did  not  accept 
a fairly  broad  range  of  input  combinations,  then  the  typical  user  (who 
didn't  have  quite  the  right  input  combination)  might  have  to  spend  a 
great  deal  of  time  deriving  and  debugging  the  intersection  equations. 
After  he  had  done  all  this,  the  computation  and  assembly  of  the  final 
QUICK  input  would  seem  a fairly  minor  task.  In  the  meantime,  he  would 
have  surely  forgotten  all  about  the  nose  code. 

Such  is  the  logic  upon  which  is  based  the  structure  of  the  nose 
geometry  code  KWIKNOSE,  the  subject  of  this  user's  manual.  KWIKNOSE 
should  be  viewed  as  a complement  to  PREKWIK:  whereas  the  latter  is 
Intended  to  model  the  overall  vehicle,  KWIKNOSE  is  specialized  for 
arbitrary  asymmetrical  nose  shapes.  With  these  two  programs,  then,  a 
flow  program  deriving  its  geometric  data  from  QUICK  could  handle,  as  an 
example,  the  symmetrically  ablated  bent  cone  of  Fig.  2.  A KWIKNOSE/ 
QUICK/BLUNT/STEIN  chain  could  carry  the  computation  to  some  point  on  the 
forecone,  and  a PREKWIK/QUICK/STEIN  chain  could  proceed  from  that  point 
through  the  bend  to  the  end  of  the  aft  cone.  (The  user  may  have  to 
construct  a fictitious  sphere-cone  forebody  on  the  PREKWIK  model  to 
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place  the  aft  body  at  the  same  position  relative  to  the  geometric  orgln 
of  the  KWIKNOSE  model.) 

The  asymmetrical  nose  geometries  which  KWIKNOSE  can  generate  are 
limited  to  cases  where  the  top,  bottom,  and  side  body  lines  can  be  , 
described  In  terms  of  sequential  series  of  straight  line  and  circle 
segments  (see  Figs.  3 and  4)  as  was  done  for  the  symmetrical  case.  The 
upper  and  lower  half  of  every  cross  section  must  be  descrlbable  in  terms 
of  an  ellipse,  though  the  ellipses  at  successive  cross  sections  need  not 
be  concentric  or  geometrically  similar  (may  have  different  centers  and 
different  major  and  minor  axes).  All  filleting,  rounding,  and  slicing  . 
capabilities  are  available  in  the  asymmetrical  case.  However,  if  the 
tip  of  the  asymmetric  nose  does  not  fall  on  the  centerline  of  the  after- 
body (or  the  first  conical  section  of  a stack) , the  resulting  geometry 
model  will  be  incompatible  with  the  PREKWIK  geometry  model  for  the 
afterbody.  However,  KWIKNOSE  does  not  distinguish  between  nose  and 
afterbody  sections,  and  a fairly  general  asymmetric  nose  atop  a stack  of 
cones  and  ogives  (with  bends)  can  be  moldeled  entirely  with  KWIOOSE 
without  reference  to  PREKWIK  (Fig.  5).  Some  additional  geometric  cal- 
culation may  have  to  be  done  by  the  user. 

The  remainder  of  this  report  describes  In  greater  detail  the  pro- 
grammed geometries  and  the  code  Itself. . A description  of  the  Input  and 
output  is  given,  and  a sample  case  Is  presented  to  Illustrate  use  of  the 
code.  The  necessary  processing  information  Is  listed,  and  certain- 
portions  of  the  programming  logic  are  elaborated  upon. 

2.0  PROGRAMMED  GEOMETRIES 

The  configurations  which  can  be  modeled  via  KWIKNOSE  fall  into  two 
general  categories:  axlsymmetrlc  and  nonaxlsymmetrlc . The  symmetric 
configuration  Is  an  arbitrary  sequence  of  conical  and  ogival  sections 
which  may  or  may  not  be  tangent  at  the  intersection  points.  The  ogival 
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sections  may  be  convex  or  concave  relative  to  the  centerline  of  the 
vehicle.  Nontangent  intersections  may,  at  the  user's  option,  be 
smoothed  with  circular  fillets  or  rounds.  The  fillets  (or  rounds)  may 
be  inserted  between  cone-cone  intersections,  cone-ogive  intersections, 
or  ogive-ogive  intersections.  The  fillets  need  not  be  very  small  in 
relation  to  the  Intersecting  curves  and  may  therefore  form  major  seg- 
ments of  the  body.  The  only  restriction  is  that  the  adjacent  segments 
being  filleted  Intersect  in  the  real  plane.  For  some  cone-ogive  and 
ogive-ogive  combinations,  there  may  exist  a maximum  filleting  radius 
beyond  which  a real  solution  does  not  exist.  In  general,  no  segment  may 
contain  a point  with  infinite  slope,  the  nosetip  being  the  one  exception. 
The  tip  of  the  nose  may  be  either  sharp  or  blunt  and  is  always  used  as 
the  geometric  origin  for  the  entire  geometry.  No  segment  may  wrap  back 
on  itself  so  as  to  require  a double-valued  radius  from  any  single  point 
on  the  centerline  axis.  The  program  dimensions  limit  the  geometry  model 
to  ten  or  fewer  segments,  and  each  fillet  occupies  one  segment  (see 
Table  1 , note  8) . Since  KWIKNOSE  does  not  distinguish  between  nose  and 
afterbody  portions  of  the  vehicle,  the  user  may  allocate  the  ten  seg- 
ments between  these  two  portions  as  needed. 

The  program  contains  logic  to  Insert  slicing  planes  into  the  basic 
axlsymmetric  vehicle  described  above.  Provisions  are  made  for  an  upper 
slice,  a side  slice,  and  a lower  slice,  each  being  composed  of  nine  or 
fewer  straight-line  segments.  In  the  side  view  the  top  and  bottom 
slicing  planes  are  seen  as  segmented  lines,  and  the  side  slicing  plane 
is  seen  as  a segmented  line  in  the  top  view.  Thus,  there  are  no  provi- 
sions for  skewed  slices.  The  slices  may  not  cross  the  centerline  of  the 
vehicle  and  may  not  intersect  each  other  inside  the  region  of  the  un- 
sliced  vehicle  at  a particular  station.  Corners  at  the  intersection  of 
adjacent  segments  may  be  rounded  or  filleted  at  the  user's  option. 
Although  the  slices  are  assumed  to  be  planar,  the  fillets  need  not  be 
small  and  may  therefore  comprise  major  portions  of  the  slicing  segments. 
It  should  be  understood  that  the  slicing  body  lines  are  input  to  QUICK 
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Independently  of  the  other  body  lines  and  thus  do  not  have  to  intersect 
the  unsllced  vehicle  at  all  points  or  any  points.  QUICK  automatically 
checks  whether  the  slicing  body  lines  intersect  the  unsllced  body  and 
chooses  the  appropriate  curve  for  the  outer  skin  of  the  vehicle.  The 
slices  therefore  may  enter  and  leave  the  body  as  needed.  There  are  no 
provisions  for  filleting  or  rounding  in  the  circumferential  direction 
(as  viewed  in  planes  normal  to  the  centerline  of  the  vehicle). 

The  asymmetrical  case  retains  all  combinations  of  cone  and  ogive 
intersection  options  (or  line-circle  intersections  in  a particular 
radial  plane),  the  rounding  and  filleting  of  the  outer  body  lines,  and 
slicing  with  rounds  and  fillets.  The  asymmetrical  vehicles  treatable  by 
KWIKNOSE  must  be  such  that  the  top,  bottom,  and  side  surface  body  lines 
may  be  represented  by  segments  of  lines  and  circles.  When  viewed  in 
cross  section,  the  vehicle  must  be  descrlbable  with  two  ellipses  (not 
necessarily  identical)  for  the  upper  and  lower  half-planes.  At  a given 
cross  section,  the  horizontal  line  dividing  the  two  half-planes  need 
not  intersect  the  centerline  (or  geometric  axis  through  the  nosetip), 
and  provisions  are  available  for  defining  this  deviation.  For  the 
asymmetrical  case  the  top  and  bottom  surface  body  lines  may  cross  the 
geometric  axis,  and  the  map  axis  (from  which  QUICK  swings  the  radius 
vector  which  must  remain  single  valued)  may  also  depart  from  the  geo- 
metric axis  if  needed.  Thus,  the  asymmetrical  option  allows,  for  example, 
the  modeling  of  an  ablated  nose  atop  a bent  cone. 

For  certain  asymmetric  configurations  there  is  a potential  inaccu- 
racy of  which  the  user  should  be  aware.  Consider  a geometry  such  as  is 
illustrated  in  Fig.  3 but  without  the  round  between  the  fore  and  aft 
sections.  This  geometric  model  was  constructed  with  the  geometric ' axis 
(X-axis)  through  the  nosetip  and  parallel  to  (but  not  coincident  with) 
the  aft  cone  centerline.  In  the  actual  hardware  (before  rounding),  the 
intersection  of  the  fore  and  aft  sections  when  viewed  from  the  side 
would  appear  as  a slanted  line,  Discontinuities  in  the  surface  slope 
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would  occur  at  this  line  in  moving  along  a body  line.  In  the  math 
model,  however,  the  slope  discontinuity  would  occur  only  on  the  upper, 
lower,  and  side  body  lines  and  would  be  smeared  everywhere  in  between. 
This  smearing  results  from  two  conditions:  (1)  QUICK  assembles  cross 
sections  in  planes  normal  to  the  geometric  axis;  (2)  as  modeled  in  Fig. 
3,  the  slanted  line  between  the  fore  and  aft  bodies  Is  oblique  to  the 
geometric  axis.  While  a carefully  constructed  manual  model  for  input  to 
QUICK  probably  could  represent  this  case  accurately,  programming  of  the 
general  case  for  arbitrary  Intersections  of  cones  and  ogives  was  consid- 
ered beyond  the  scope  of  KWIKNOSE. 

The  user  has  one  alternative  If  he  requires  an  accurate  model  from 
KWIKNOSE:  provide  the  input  data  relative  to  a geometric  axis  normal  to 
the  slanted  line.  Note,  however,  that  the  geometric  axis  might  pene- 
trate the  skin  and  require  bending  of  the  map  axis  to  keep  QUICK* s 
radius  vector  single-valued.  In  addition,  once  the  geometric  axis  has 
been  made  normal  to  a particular  Intersection  plane,  no  other  similar 
intersection  with  a different  slant  angle  can  be  modeled  with  equal 
accuracy.  For  two  conical  sections  of  different  half-angles  and  oblique 
centerlines,  the  intersection  "plane"  is  really  curved.  Thus,  reori- 
enting the  geomtric  axis  can  only  reduce  (but  not  eliminate)  the  inaccu- 
racy in  modeling. 


3.0  PROGRAM  DESCRIPTION 


3.1  INPUT 

The  input  to  KWIKNOSE  is  described  in  Tables  1 through  4,  and  a 
listing  of  the  KWIKNOSE/QUICK  Jobstream  with  the  input  for  the  sample 
case  of  Figs.  4 through  6 is  given  in  Table  5.  It  is  suggested  that  the 
user,  prior  to  assembling  the  input  deck,  make  a qualitatively  accurate 
sketch  of  the  upper,  lower,  and  side  body  lines  and  indicate  the  known 
geometric  data.  This  task  has  been  found  useful  in  choosing  appropriate 


11 


AEDC-TR-77-89 


Input  options  and  deciding  on  the  values  of  the  logical  Input  variables 
(DIR  and  SLOPE),  which  are  needed  by  the  code  for  selection  of  a segment 
Intersection  from  multiple  roots. 

The  first  Input  card  Is  for  Identification  and  may  contain  any  80 
characters,  though  QUICK  will  carry  only  60  of  them.  The  remaining 
input  is  via  namelist  /INPUT/.  For  a symmetrical  body  (with  or  without 
asymmetrical  slices),  only  one  namelist  Input  is  required.  For  an  , 
asymmetrical  body  three  separate  namelists  must  follow  the  single  Identi- 
fication card.  The  first  must  define  the  upper  body  line  (ZUCL,  see 
Figs.  7 and  8),  the  second  the  lower  body  line  (ZLCL),  and  the  last  the 
side  body  line  (YSID).  Note  that  the  Input  variable  Z Is  used  for  the 
radial  coordinate  regardless  of  whether  the  Input  data  represents  the  Y 
or  Z direction.  The  slicing  body  lines  (ZUSL,  ZLSL,  and  YSSL)  and  the  Z 
coordinates  of  the  side  and  map  axis  (ZSID  and  ZMAPAX)  may  be  Included 
in  any  of  the  three  namelists,  if  they  are  input  at  all.  The  upper, 
lower,  and  side  body  lines  are  processed  independently  of  each  other  and 
need  have  no  particular  correspondence  with  each  other  except  for  the 
same  final  X coordinate  at  the  end  of  the  last  segment.  Thus,  each  may 
be  an  entirely  different  combination  of  lines,  circles,  and  fillets. 

The  sign  notation  for  input  coordinates  and  angles  is  that  of  a cylin- 
drical coordinate  system.  Angles  (TH(I))  are  positive  away  from  the 
centerline  and  negative  pward  it  relative  to  a line  parallel  to  the 
centerline. 

The  input  rules  for  the  three  main  body  lines  are  somewhat  different 
chan  for  the  slicing  body  lines  and  will  be  discussed  separately.  For 
each  of  the  three  main  body  lines  (upper,  lower,  side)  to  be  defined, 
the  user  must  input  the  element  of  the  character  array  CURVE(I)  as 
'LINE'  or  'OGIV'  for  each  segment.  Based  on  this  choice  for  each  indi- 
vidual segment,  the  user  then  chooses  an  input  option  from  Table  2 for 
each  of  the  segments,  preferably  an  option  which  matches  the  combination 
of  geometric  data  he  has  for  the  segment.  If  any  redundant  data  is  to 
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be  Input,  case  17  or  18  (see  ICASE  values  In  Table  2)  must  be  chosen. 

If  the  redundant  data  is  not  self  consistent  to  within  ERRMAX,  execution 
will  be  halted.  All  data  items  indicated  by  X's  in  Table  2 for  a parti- 
cular ICASE  value  must  be  input  or  execution  will  be  halted.  Note  that 
ICASE (I)  is  not  an  input  variable.  The  data  elements  described  in  Table 
1 are  then  coded  according  to  the  standard  rules  for  namelist  input  as 
found  in  any  Fortran  manual.  The  user  should  note  that  X(I),  Z(I),  and 
TH(I)  are  the  coordinates  and  slope  at  the  end  of  the  Ith  segment 
(CURVE(I)).  Further,  certain  input  cases  may  not  be  used  following 
certain  other  cases.  For  example,  cases  5 and  9 constitute  inputting 
only  the  equation  of  the  segment  and  in  themselves  do  not  Include  suffi- 
cient information  to  locate  the  intersection  of  the  end  of  the  segment. 
Thus,  those  cases  which  assume  that  the  coordinates  of  the  end  of  the 
previous  segment  are  known  may  not  follow  cases  5 and  9.  Such  limita- 
tions, the  need  for  which  will  be  clarified  in  Section  3.4,  are  given  in 
the  footnotes  to  Tables  1 and  2. 

The  direction  indicator  variable  DIR(I)  requires  additional  clari- 
fication. Whenever  ogival  segments  are  involved,  the  possibility  of 
multiple  Intersections  with  adjacent  segments  exists.  In  addition,  when 
an  ogive  is  to  be  Inserted  between  two  given  points,  there  are  usually 
two  circle  equations  which  will  satisfy  all  the  numeric  data.  The  user 
therefore  must  indicate  to  KWIKNOSE  which  of  the  possibilities  is  to  be 
chosen  using  DIR(I).  Table  2 Indicates  which  input  cases  require  DIR(I) 
and  Table  3 defines  the  acceptable  values  of  DIR(I)  for  each  ICASE 
value.  The  nomenclature  for  these  acceptable  values  is  defined  in  Table 
4.  The  ’OPUP',  'OPDN',  'TOPS',  and  'BOTS'  values  apply  to  all  three 
main  body  lines  (upper,  lower,  and  side)  as  though  they  were  rotated 
about  the  centerline  for  viewing  as  oriented  in  Table  4.  That  is, 

DIR(I)  = 'OPUP'  indicates  that  the  centerline  sees  a convex  segment  and 
'OPDN'  a concave  segment.  Similarly,  'TOPS'  always  indicates  the  inter- 
section farther  from  the  centerline  than  the  'BOTS'  intersection. 


13 


AEDC-TR-77-89 


If  a nontangent  Intersection  Is  to  be  filleted  or  rounded,  this 
fact  Is  Indicated  by  giving  the  radius  of  the  round  or  fillet  RHO(I)  to 
be  Inserted  at  the  end  of  the  Ith  segment.  Fillets /rounds  are  not 
Included  when  numbering  segments  for  Input,,  and  no  data  other  than  the 
radius  are  input  for  the  fillet.  The  Intersections  and  center  are 
computed  Internally  by  KWIKNOSE.  There  are  several  restrictions  on  the 
filleting  (or  rounding)  process.  Segments  already  tangent  may  not  be 
filleted.  Only  adjacent  segments  which  actually  intersect  may  be 
filleted  because  the  location  of  this  Intersection  is  used  to  choose,' 
between  multiple  roots.  Further,  some  types  of  Intersections  may  have 
maximum  filleting  radii  beyond  which  real  solutions  do  not  exist.  If 
KWIKNOSE  detects  any  of  these  conditions,  usually  through  negative 
square  roots  or  division  by  zero,  execution  will  be  halted. 

For  as3nmnetrlcal  bodies,  such  as  a bent  cone,  the  user  may  opt  to 
define  the  Z-coordinate  of  the  side  (ZSID,  Figs.  6 and  8)  as  different 
from  zero,  in  which  case  he  must  Input  two  or  more  values  of  XSID  and 
ZSID  as  described  in  Table  1.  XSID(1)  and  ZSID(1)  must  both  be  zero; 
and  if  ZSID  is  to  be  Input  at  all,  at  least  two  values  of  each  (XSID(1) 
and  XSID(2),  ZSID(1)  and.ZSID(2)  must  be  given.  If  the  vehicle  is  bent 
so  far  that  the  geometric  axis  penetrates  the  skin,  the  user  must  bend 
the  map  axis  via  XMAPAX  and  ZMAPAX  so  that  a radius  vector  swung  about 
the  map  axis  at  any  station  along  the  vehicle  will  be  single-valued. 

Care  should  be  taken  in  using  models  generated  by  QUICK  in  which  the  map 
axis  is  bent.  The  stock  version  of  the  QUICKLOK  phase  of  QUICK  returns 
to  the  user  program  radius  vectors  relative  to  the  map  axis  rather  than 
the  geometric  axis.  Thus,  a bent  cone  whose  map  axis  follows  the  cen- 
terlines of  the  cones  may  appear  to  the  user  program  as  not  being  bent 
at  all.  Input  follows  the  same  rules  as  for  XSID  and  ZSID  above.  Note 
that. both  ZSID  and  ZMAPAX  may  be  negative.  These  variables  are  used 
only  if  TYPE  = ’ASYM' , for  which  case  ZSID  defaults  to  zero  if  not 
input,  and  ZMAPAX  defaults  to  ZSID  if  not  input.  These  input  variables 
may  be  included  in  any  one  of  the  three  namelists. 
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Input  for  the  slicing  planes  as  well  as  the  processing  of  the  data 
by  KWIKNOSE  is  completely  independent  of  the  main  body  lines.  Any 
combination  of  the  three  slices  may  be  used  regardless  of  the  symmetry 
aside  from  the  slices  (the  value  of  TYPE  has  nothing  to  do  with  symmetry 
of  the  slices).  The  minimum  slicing  plane  input,  if  present  at  all,  is 
the  coordinates  of  the  beginning  and  end  of  one  segment.  At  the  user's 
option,  different  combinations  of  X and  Z (or  Y)  coordinates  and  slope 
angle  may  be  input.  Note  that  TUSL(I)  is  the  slope  angle  of  the  segment 
ending  at  XUSL(I),  ZUSL(I),  for  I > 1.  Thus,  TUSL(1)  = TUSL(2),  and 
similarly  for  the  other  two  slices.  Redundant  data  will  be  checked  for 
consistency;  but  if  any  one  of  the  three  is  not  input  for  a given  seg- 
ment, it  will  be  computed.  The  slicing  planes  need  not  cover  the  entire 
length  of  the  vehicle,  but  the  first  and  last  segments  will  be  extrapo- 
lated to  the  extremities  if  not.  Nontangent  Intersections  of  adjacent 
segments  may  be  filleted  or  rounded  by  inputting  the  fillet/round  radius 
as  indicated  in  Table  1.  Note  that  RHOUS(I)  fillets  the  intersection  at 
XUSL(I),  which  is  the  end  of  the  I-lst  segment;  otherwise,  all  the  rules 
of  body  line  filleting/rounding  apply. 

The  remaining  input  variables  are  program  control  integers.  IPRINT 
set  different  than  zero  will  cause  KWIKNOSE  to  add  cards  to  the  QUICK 
input  deck  to  exercise  the  math  model  and  generate  printouts  of  coordi- 
nate and  slope  data  at  various  points  along  each  body  line  segment. 
IPRINT=0  (the  default)  suppresses  such  printout.  The  variable  IR  is  the 
unit  number  (Data  Set  Reference  Number)  of  the  facility  card  reader,  and 
IW  is  the  unit  number  of  the  printer.  IF  is  the  unit  number  of  the 
KWIKNOSE  output  file  of  card  images  to  be  passed  to  and  read  by  the 
QUICKDEF  phase  of  QUICK.  IP  is  the  unit  number  for  the  facility  card 
punch.  If  IF  and  IP  are  different,  the  data  from  unit  IR  is  copied  as 
card  images  with  an  80-character  array  to  IF,  from  which  the  input  cards 
are  then  printed  on  IW.  IF  is  then  rewound  and  read  as  namelist  data. 
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3.2  OUTPUT 

Table  6 is  the  KWIKNOSE  printout  for  the  sample  case.  The  printout 
is  composed  of  three  sections:  the  input  data,  the  results  of  the 
computations,  and  the  QUICK  Input  file.  The  first  section,  labeled  "INPUT 
CARDS"  is  a listing  of  the  data  cards  exactly  as  punched.  The  data 
between  the  first  two  lines  of  periods  is  the  input  data  after  the  first 
namelist  has  been  read.  Notice  that  geometric  variables  not  input 
through  the  namelist  default  to  1.E70.  This  defaulting  is  critical  to 
the  program's  logic  and  should  not  be  trampered  with  by  inputting  user- 
chosen  default  values  which  have  no  geometric  meaning.  The  printout 
from  this  point  to  the  third  line  of  periods  contains  the  results  of  the 
geometric  computations.  First,  the  number  of  segments  exclusive  of 
rounds  and  fillets  is  printed' followed  by  the  code’s  determination  as  to 
the  input  options  the  user  has  chosen  for  each  segment.  Next,  the 
results  of  the  computation  of  segment  intersections  and  equations  are 
printed.  In  the  printout  labeled  "SEGMENT  END  POINT  COORDINATES  AND 
SLOPES",  X and  Z are  the  coordinates  at  the  beginning  of  the  segment, 
except  for  the  last  point  (with  values  less  than  1.E70),  which  is  the 
end  of  the  last  segment.  TA  and  TF  are  the  slope  angles  in  degrees  at 
the  aft  and  fore  ends,  respectively,  of  the  segment.  The  segment  equa- 
tions in  the  next  section  of  the  printout  are  defined  in  terms  of  the 
center  coordinates  and  radius  for  ogival  sections,  and  slope  angle  and 
Z- intercept  for  straight-line  segments.  The  X value  printed  here  is  the 
coordinate  of  the  end  of  the  segment.  In  these  terms  the  line  and  ogive 
equations  are 

line;  Z = X tan  T + B 
ogive:  (X-H)^  + (Z-K)^  = R^ 

Next,  the  results  of  the  fillet  and  round  computation  is  printed.  The 
equation  of  the  fillet  is  defined  as 

fillet:  (X-KSI)^  + (Z-ETA)^  = RHO^ 
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and  straddles  the  printed  X value.  XI,  Z1,  T1 , and  X2,  Z2,  T2  are  the 
coordinates  and  slope  angle  (degrees)  at  the  fore  and  aft  intersections, 
respectively,  that  is,  the  points  of  tangency  with  the  intersected 
segments.  The  fillet/round  intersections  are  then  printed  as  additional 
segments  included  with  the  previous  unfilleted  body  line  under  the 
heading  "FINAL  QUICK  ABRAY".  Note  that  the  original  intersection  point 
has  been  removed  from  the  arrays.  This  completes  the  geometric  com- 
putation for  the  upper  body  line  (ZUCL).  When  TYPE  = 'ASYM'  in  the 
namelist  input,  this  computation  and  printout  are  repeated  first  for  the 
lower  body  line  (ZLCL)  and  then  for  the  side  (YSID).  The  third  printout 
includes  the  results  of  the  geometrical  calculation  for  slices.  The 
section  of  the  printout  headed  "SLICE  GEOMETRY"  gives  the  slope  angle 
for  each  segment  ending  at  the  indicated  coordinate  values.  These 
results  are  exclusive  of  rounds  or  fillets,  for  which  immediately 
follows  the  computational  results.  BI  and  BJ  are  the  Z- intercepts  of 
the  two  lines  formed  by  the  locus  of  the  fillet/round  center  while 
rolling  on  the  two  intersecting  segments.  A3  and  A4  are  the  Z and  X 
coordinates  of  the  fillet/round  center;  and  XI,  Z1,  T1 , etc.,  are  the 
tangency  points  as  described  for  the  main  body  line  fillet/round 
printout.  The  final  portion  of  the  computational  printout  is  the 
slicing  body  lines  with  the  rounds  and  fillets  Inserted.  The  KWIKNOSE 
printout  concludes  with  a listing  of  the  card  images  of  the  QUICK  input 
file.  This  listing  will  not  be  generated  if  IF  = IP.  The  contents  of 
this  output  deck  may  be  understood  from  Table  7 and  Refs.  1 through  3 and 
will  not  be  discussed  here.  The  user  should  understand  that  this  deck 
comprises  the  input  to  the  QUICKDEF  phase  of  QUICK  and  not  the  QUICKLOK 
phase,  which  is  the  portion  of  QUICK  employed  in  the  fluid  dynamics  pro- 
grams to  obtain  geometric  information. 

3.3  PROCESSING  INFORMATION 

KWIKNOSE  is  a Fortran  IV  program  containing  about  2,000  cards 
(Table  8).  The  code  was  developed  on  an  IBM  370/165  computer  system 
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under  the  H Level  21.7  compiler.  Compilation  of  the  entire  deck  re- 
quires 230  K bytes  of  core  and  40  sec  of  central  processing  unit  (CPU) 
time.  Execution  of  the  code  requires  80  K bytes  of  core,  and  all  geometries 
processed  to  date  have  required  less  than  2 sec  of  CPU  time.  Users 
should  note  that  because  of  the  character  manipulation  performed  in  the 
program,  it  cannot  be  compiled  under  the  G level  21  compiler.  Execution 
requires  one  peripheral  storage  device  for  unit  IF  (defaulted  to  unit 
1).  This  should  be  a sequential  file  formatted  for  80-character  records 
as  illustrated  in  Table  5. 

3.4  PROGRAMMING  CONSIDERATIONS 

The  overall  flow  of  logic  in  KWIKNOSE  is  basically  simple  and  can 
be  best  understood  via  the  series  of  CALL  statements  in  the  main  program 
(Table  8)  and  the  subroutine  description  of  Table  9.  Execution  begins 
with  the  initialization  of  certain  data  (subroutine  INIT),  after  which 
the  first  namelist  is  read  (INPT).  Next,  the  input  option  chosen  by  the 
user  is  determined  for  each  segment  (CASE)  in  terms  of  the  value  of 
ICASE(I).  This  variable  is  then  used  to  control  branching  in  deter- 
mining the  segment  equations  (SEGE)  and  locating  the  segment  intersec- 
tions (INTR).  This  completed,  the  results  of  the  computation  to  this 
point  are  copied  into  the  QUICK  arrays  XQ,  ZQ,  and  TQ  (QARY)  and  printed 
(PRNT).  Next,  the  round/fillet  computations  are  performed  and  added  to 
the  QUICK  arrays  (ROFL).  If  slices  are  present  (indicated  by  ISL  ^ 0), 
the  necessary  computations  are  performed  and  the  results  assembled  in 
the  QUICK  slicing  array  V (SLCE).  If  rounds  or  fillets  are  present  on  . 
the  slices,  these  computations  are  performed  and  the  results  Inserted 
into  the  original  arrays  (RFSL),  If  the  case  is  symmetrical,  the  QUICK 
input  card  image  file  is  generated  (QUCK)  and  printed  (LIST) . If  the 
case  is  asymmetrical,  the  completed  body  line  with  rounds  and  fillets  is 
copied  into  the  three-dimensional  QUICK  array  VB  (COPY).  When  all  three 
main  body  lines  are  complete,  the  slicing  with  optional  rounds /tlliets  is 
performed  and  copied  into  the  three-dimensional  arrays.  Next,  ZSID  and 
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ZMAPAX  are  defaulted  if  not  Input,  and  the  remaining  body  lines  are 
checked  for  validity  (BLCK).  Finally,  the  QUICK  input  file  for  the 
asymmetrical  case  is  generated  (QK3D)  and  listed  (LIST),  thus  completing 
the  processing. 

For  a more  detailed  examination  of  the  code,  the  variable  descrip- 
tion in  Table  10  may  be  useful.  Users  Interested  in  the  intersection 
computation  should  find  the  logic  as  coded  in  SEGE  and  INTR  fairly 
straightforward,  but  a few  clarifying  comments  about  the  rounding  and 
filleting  logic  are  necessary.  The  straightforward  approach  to  locat- 
ing the  center  and  tangency  points  of  a circle  filleting  two  other 
circles  is  to  use  analytic  geometry,  by  writing  four  circle  equations 
and  two  tangency  equations  in  terms  of  six  unknowns.  However,  this  is 
equivalent  to  a polynomial  of  degree  48  for  one  unknown,  which  is  an 
absurd  thing  to  try  to  solve.  The  approach  used  in  ROFL  is  to  write  the 
equation  of  the  locus  of  fillet  centers  as  the  fillet  rolls  on  the  two 
adjacent  segments,  where  the  intersection  of  these  two  loci  of  centers 
is  the  center  of  the  filleting  circle.  This  approach  reduces  the  problem 
to  second  degree  but  requires  additional  logic  to  decide  when  to  add  or 
subtract  the  fillet  radius  to  or  from  the  segment  radius  (ogive)  or  Z- 
intercept  (line).  Fortunately,  the  Information  necessary  to  make  this 
decision  can  be  extracted  from  the  slopes  of  the  adjacent  segments  at 
the  intersection  point.  For  the  line-ogive  intersection,  the  final 
second  degree  problem  can  then  be  solved  by  choosing  the  set  of  tangency 
points  which  straddle  the  intersection  point.  For  the  ogive-ogive 
Intersection,  the  filleting  circle  is  chosen  whose  center  is  closest  to 
the  original  Intersection  point.  The  advantages  of  using  the  loci-of- 
centers  approach  are  clear;  but  the  major  disadvantage  is  that  segments 
which  do  not  intersect  cannot  be  filleted,  whereas  the  general  48th 
degree  problem  includes  this  case. 

Because  of  the  complexity  of  the  above  logic  as  coded  and  to  the 
many  input  options  available,  KWIKNOSE  contains  many  small  blocks  of 
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coding  to  check  on  the  validity  of  the  input  and  the  running  computa- 
tion. When  a violation  of  the  rules  is  detected,  a condition  code  is 
generated  with  the  Fortran  statement  STOP  xxxx,  where  xxxx  is  the  four-  * 
digit  condition  code.  Execution  is  halted  and  an  error  message  Is  < 
usually  generated.  A program-generated  condition  code  without  an 
accompanying  error  message  would  Indicate  either  a programming  error  or 
computer  error  (unlikely) . The  condition  code  is  usually  the  same  as 
the  statement  number  of  the  format  that  generated  the  accompanying' 
message.  For  five-digit  statement  numbers,  the  left  most  digit  is 
truncated  to  obtain  the  condition  code. 

The  accuracy  of  the  code's  geometry  computation  has  been  checked  by 
comparison  with  manual  computation  on  an  electronic  calculator.  To  the 
best  of  the  author's  knowledge,  each  block  of  coding  has  been  checked 
against  a tailored  check  case.  However,  the  user  is  encouraged  to 
carefully  scrutinize  the  final  QUICK  math  model  for  accuracy.  The 
author  has,  to  date,  constructed  about  30  manual  card  input  decks  to 
QUICK  and  hundreds  of  PREKWIK  and  KWIKNOSE  decks,  and  no  errors  have 
been  traced  to  QUICK.  Thus,  errors  in  the  final  math  model  will  most 
likely  be  caused  by  either  problems  in  the  setup  codes  or  their  input 
decks.  Checking  of  the  QUICK  math  model  can  be  effectively  accomplished 
if  a plotting  device  is  available  with  stock  software.  With  very  little 
effort,  the  author  has  constructed  a cathode  ray  tube  (CRT)  plotting 
program  using  the  QUICKLOK  portion  of.  QUICK  and  has  found  it  indispens- 
able in  debugging  the  geometry  models.  This  approach  is  highly  recom- 
mended to  all  users. 
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4.0  SUMMARY 

A user's  manual  has  been  presented  for  a Fortran  computer  program 
KWIKNOSE  which,  for  certain  axisymmetric  and  asymmetric  nose  configura- 
tions, provides  simplified  geometric  input  to  the  Grumman  QUICK- geometry 
code,  which  in  turn  provides  a geometric  math  model  to  various  numeri- 
cal flow-field  programs.  A discussion  was  presented  as  to  how 
KWIKNOSE  may  easily  be  used  to  model  any  vehicle  whose  upper,  lower,  and 
side  body  lines  can  be  described  in  terms  of  sequential  segments  of 
circles  and  straight  lines  and  whose  upper  and  lower  half  cross  sections 
can  be  modeled  with  ellipses.  It  was  further  shown  how  the  body  may 
be  sliced  on  three  sides.  A description  was  given  of  the  input  and 
output,  and  a sample  case  was  presented  to  illustrate  use  of  the  code. 
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Figure  4.  Oblique  view  of  asymmetrical  nose  from  sample  case. 
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Figure  5,  Oblique  view  of  asymmetrically  ablated,  bent,  and  sliced  cone  (sample  case). 
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ZUCL,  ZUSL,  YSID,  etc.  are 
logical  names  used 
Internally  by  KWIKNOSE. 
Some  are  also  input 
variables  (Table  1). 


lie  10 


Lower  Centerline  Body 
Line  (ZLCL)  LB  - 2* 


Figure  7.  Body  line  notation. 
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Table  1.  Input  to  KWIKNOSE  via  Namelist/INPUT/ 


Namelist 

Variable 

Val-Uf! 

Default 

Value 

Heanlng 

TYPE^ 

X(I) 

(O,*) 

1 ,E70 

Centerline  distance  to  end  of  I th  segment  (I&IO) 

R*4 

Z(I) 

1.K70 

Radial  distance  (normal  to  contorllnc)  to  end  of  Ith  segment 

R*4 

TH(I)® 

(-90,90) 

1 .E70 

Slope  angle  (deg)  at  end  of  Ith  segment,  poeitive  away  from  ccntcriinr 

R*4 

CURVE( I ) 

'LINE' 

•bbbb'l 

Segment  I Is  a straight  line 

R*4 

•*CIV' 

Segment  I is  an  ogive 

B(I)^ 

(0,") 

1 .E70 

Radius  of  ogive  of  Ith  segment  (CURVK( I)-'dGlV' ) 

R*4 

H(I)’ 

(—CO^no) 

1.E70 

X coordinate  of  center  of  ogive  ((nJRVE(  I)**'4GIV' ) 

R*4 

K(I)® 

(.00^ «) 

I.E70 

Z coordinate  of  center  of  ogive  (C7URVE(I)-'4GIV' } 

R*4 

B(I)* 

(_«),  ‘90) 

1.E70 

Z intercept  of  line  segment  ( CURVE-' LINE ' ) 

R*4 

OIR(I)^ 

See  Table  3 

•#PI)N' 

Direction  indicator  for  CURVLC  I)-*^G1V'  when  oft  intersection  with  adjacent  segment  may  have  two  solutions 

R*4 

SUPEd) 

tipl 

•F' 

Indicates  ogive  is  to  be  made  tangent  to  aft  end  of  previous  segment  (do  not  iiso  for  segment  1-1) 

L*1 

ERRHAX 

>10"' 

S.li-9 

Error  criterion  for  checking  conslstrncy  of  redundant  Input  data 

R*4 

RH#(I)®’® 

(0,*) 

1 .E70 

Fillet/round  radius  at  end  of  Ith  segment 

R*4 

XUSL(I)® 

|0,«>) 

1 .E70 

'I 

R44 

ZUSL(I)® 

lO,") 

1.E70 

>X  and  Z coordinates  and  slope  (deg)  of  upper  slicing  plane.  ISIO 

R*4 

TUSL(I)® 

(-90', +90') 

1 .E70 

; 

K*4 

XLSLd)® 

lo.*) 

1.E70 

) 

R*4 

ZLSLd)® 

lO.o) 

1.E70 

>X  and  Z coordinates  and  slope  (drg)  nf  lower  slicing  plane.  I<I() 

R*4 

TLSLd)* 

(-90*  ,490') 

1.E70 

\ 

R*4 

XSSLd)® 

o.«) 

1.E70 

) 

R*4 

VSSLd)* 

lO.w) 

1.E70 

>X  and  Y coordinates  and  slope  (deg)  of  side  slicing  plane,  1£10 

R*4 

TSSLd)® 

(-90' .90') 

1 .E70 

R*4 
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Table  1.  Continued 


Namelist 

Variable 

Value 

Default 

Value 

Heanlog 

TYPE^ 

RH^USd)^ 

(0,»J) 

1.K70 

^Flllet/round  radius  on  the  upper,  lower,  and  side  slices, 

R*4 

RH^LSd)^ 

<0,’=) 

1.E70 

>re8DCctl«relv.  RH^ — (1)  fillets  or  rounds  the  bcRlnninR  of 

R*4 

RB^SSd)^ 

(0,») 

1.E70 

)segment  1 and  the  end  of  segment  1-1.  2<I<10. 

R*4 

XSIDd)^ 

(0,») 

X and  Z coordinates  of  side  surface  body  line.  Input  only  If  body  Is  asyumetrlcal  and  side  surface 

R*4 

ZSlDd) 

(.00,  «0) 

body  line  is  not  in  some  horizontal  plane  as  the  geometric  X axis.  Defaults  to  X axis  (see  Fig.  B). 

R«4 

XIMPAXd) 

X and  Z coordinates  of  QUICK  map  axis  (the  axis  relative  to  which  the  vehicle  must  be  descrlbable  In 

R*4 

ZMAPAXd) 

(-»,  *) 

terms  of  a single— valued  radius  vector).  Input  only  if  body  is  asymmetrical  and  single— valued  radius 

requirement  cannot  be  met  by  default.  Default  to  XSID,  ZSID. 

R*4 

TYPE 

•SYHH* 

•ASYM* 

'SYHH* 

Vehicle  is  axially  symmetric.  Provide  input  for  one  body  line. 

Vehicle  is  as^mnetrlc.  Provide  input  for  three  body  lines: 

DCL,  LCL,  and  SID  (see  Fig.  8).  Input  the  needed  combinations  of  CURVE.  X,  Z,  lU.  R.  B,  K,  B, 

DIR,  and  SL^PE  in  three  successive  ftINPUT  namelists.  Provide  slicing  data*  XSID,  ZSID,  XHAPAX • 

ZHAPAX  with  any  one  of  three  if  required.  First  of  three  namelists  is  upper  centerliney  second  is 

lower  centerline,  third  is  side. 

R*4 

11.991 

5 

Unit  number  for  namelist  input 

1*4 

iwio 

|l,99l 

6 

Unit  number  for  printed  output 

1*4 

11.991 

^ZP 

-IP 

1 

Unit  number  for  output  of  QUICK  input  data 

Unit  IF  will  be  rewound  and  printed  oa  unit  IW 

Unit  IF  will  not  be  rewound  or  printed.  Punch  QUICK  input  data 

1*4 

IP 

il.99| 

7 

Unit  number  for  card  punch  (need  not  be  dummied  if  not  used) 

1*4 

IPRINT 

=^0 

-0 

0 

Cards  will  be  appended  to  QUICK  input  (unit  IF)  to  exercise  QUICK  math  model  for  checking 

No  exercise  cards  generated 

1*4 
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Table  1.  Concluded 


Noteu: 


^'WU‘  - four  bl»nk  churacters. 

^TB(I)  and  B(I)  constitute  segment  equation  for  CURVECI)  -'LINE',  l.e.,  '/(I)  - X<I)  * TAN(TH(I))  + B(I). 

^R(I),  H(I),  and  K(I)  constitute  segment  equation  for  CimVE(l)-’dGIV',  l.e.,  (X(I)-H(I))**2+(Z(I)  -K(I))**2  - R(l)**2. 

^See  Tables  3 and  4 for  complete  explanation. 

^Do  not  Input  these  variables  unless  rounds,  fillets,  or  slices  are  present. 

^SID(l)  and  ZSlD(l)  nust  both  be  input  as  zero.  If  input  at  all,  at  least  two  values  of  each  must  be  input. 

7 

Variable  type:  L • logical,  I - integer,  R • real.  Number  following  « is  number  of  bytes. 

^When  defining  a series  of  segments  via  X,  Z,  TH,  R,  H,  K,  and  B,  do  not  skip  an  1 value  when  RH0(I)  is  given.  RHdCI)  will  fillet/round  the 
intersection  of  the  1th  and  I+lst  segments.  For  example,  to  fillet/round  two  line  segments,  input  X(l),  Z(l),  Z(2),  Z<2),  and  RH^d). 

Note  further  that  each  fillet/round  reduces  by  one  the  total  number  of  allowable  segments:  ten  segments  with  two  fillets  will  overrun 
QUICKCS  dimensions,  whereas  eight  segments  with  two  fillets  will  not.  Though  fillets/rounds  are  not  treated  as  segments  in  the  namelist 
input,  they  are  treated  as  separate  segments  in  the  output  after  their  intersections  with  adjacent  segments  have  been  computed  and  inserted 
into  the  XQ,  ZQ,  and  TQ  arrays. 

9 

May  be  negative  when  TYPE  -'ASYM*  and  if  map  axis  is  bent. 

^^If  the  unit  numbers  are  to  be  changed,  they  must  also  be  changed  In  subroutine  INXT  because  data  are  read  from  IR  and  bitten  on  IW  before 
the  namelist  data  arc  read  as  namelists. 

^^If  IF  were  set  equal  to  1P(07)  in  the  namelist,  the  input  data  would  still  be  read  from  unit  1R(05)  copied  to  unit  01. 

Note:  Character  means  capital  letter  O 
"0”  means  numeric  zero  0 

I ] means  end  points  arc  included;  ()  means  end  points  are  excluded. 
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Table  2.  KWIKNOSE  Input  Options'*  for  Three  Main  Body  Lines  (ZUCL,  ZLCL,  YSID) 


CURVE(Z) 

X(I) 

Z(I) 

H(I) 

K(I) 

K(I) 

SL«FE(I) 

TH(I) 

B(I) 

DIR(I)^ 

Hay  Follow 
Case  S.  9.  197 

lUSE(I)^ 

'♦GIV 

X 

X 

— 

— 

z 

— 

— 

— 

X 

No 

1 

582 

X 

X 

— 

— 

— 

X 

— 

— 

— 

No 

2® 

774 

X 

— 

— 

— 

X 

X 

— 

— 

X 

Mo 

3® 

^34 

X 

X 

— 

— 

— 

— 

X 

— 

— 

Ho 

4 

526 

— 

— 

X 

X 

z 

— 

— 

— 

X 

Tes 

5* 

624 

X 

X 

— 

— 

X 

— 

X 

— 

X 

Tes 

10^ 

590 

X 

— 

X 

X 

x 

— 

— 

— 

X 

Yes 

11® 

626 

— 

X 

X 

X 

z 

— 

— 

— 

X 

Yu 

12® 

628 

— 

— 

X 

X 

z 

— 

X 

— 

X 

Tea 

13® 

632 

X 

X 

X 

X 

z 

— 

X 

— 

— 

Tea 

18®*^ 

638 

— 

— 

— 

— 

X 

X 

— 

. — 

X 

No 

19® 

832 

X 

X 

— 

— 

X 

X 

— 

— 

X 

Yea 

21®'® 

838 

■LINE' 

X 

X 

— 

— 

— 

— 

— 

— 

— 

No 

6 

6 

B 

B 

B 

B 

No 

Ho 

8 

10 

12 

■9 

B 

H 

B 

Tea 

9* 

136 

D 

B 

B 

B 

Tes 

14® 

14 

D 

B 

^9 

B 

Yea 

15® 

138 

H 

X 

B 

H 

B 

H 

B 

Yes 

16® 

140 

B 

B 

B 

B 

B 

y» 

1,®.7 

142 

B 

B 

B 

B 

^Bi 

B 

B 

Bi 

No 

20® 

262 

512® 

1 

1 ^ 1 

1 1 

1 1 

1 1 

256 

1 ^ 1 

1 1 

^Onc  of  these  options  asist  be  chosen  for  each  segment.  Insufficient  or  redundant  data  will  cause  eaecution  to  be  halted. 
^CASE(l)  and  L are  uqt  Input  variables  but  are  used  for  program  control.  See  subroutines  CASEv  SEGE.  and  HVTR  In  Table  8. 

^IR(l)  Is  used  only  when  multiple  solutions  or  Intersections  can  occur.  When  multiple  roots  are  equals  DIR(I)  is  ignored. 
Site  Tables  3 and  4. 

^These  options  are  not  acceptable  for  last  segment  of  vehicle. 

4 5 

* Use  one  of  these  options  if  previous  seeiient  was  given  via  ICASB*5s9t  or  19. 

^These  options  are  not  acceptable  for  first  segaient  of  vehicle. 

^Theae  options  contain  redundant  data,  which  If  not  self  consistent  to  within  ERRMAX.  will  cause  execution  to  halted. 
^CUI(VE(I)-«GIV',  L>S12;  -'LIME;  L-0. 
q 

If  X(I)<1.E70.  IF  X(I)^1.B70  (default),  L*0.  Same  for  renalalng  variables  but  with  higher  powers  of  2. 
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Table  3.  Input  Description  for  Variable  DIR(I)'' 


ICASB(l) 

ICASE(I+1) 

ICASK(I-l) 

DIR(l) 

Meaning 

1,  3,  10,  11,  13 

— 

— 

'♦PUP'^ 

Ogive  opens  up 

**PDN*^ 

Ogive  opens  down  . 

5 

5 

— 

'P*RE' 

Use  fore  solution 

•aste* 

Use  aft  solution 

•tws’ 

Use  top  solution  Use  when  X coordinate  of  fore 

'B«S’ 

Use  bottom  solution  and  aft  solutions  are  equal 

5 

9 

1 

'F9RE’ 

Use  fore  solution 

9 

5 

_ f 

•APTE’ 

Use  aft  solution 

12 

— 

— 

'F*RZ’ 

Use  fore  solution 

'AFTE* 

Use  aft  solution 

20 

— 

5 

'T*PS’ 

Use  top  solution 

'B«S' 

Use  botcoo  solution 

21 

— 

’ i 

tjpUptZ 

Ogive  opens  up 

5 

21 

1 

'*PDN'^ 

Ogive  opens  down 

21 

— 

9 

— 

DIR(I)  not  needed  | 

19 

— 

— 

'OFF*' 

Ogive  opens  up,  use  fore  solution  ^ 

•UPAF’ 

Ogive  opens  up,  use  aft  solution 

1 

'irPT*‘ 

Ogive  opens  up,  use  top  solution 

1 For  DIR(l)  notation,  assume 

'UP'-'9PUP',  ’DH’-'^PDN', 

'UPB*' 

Ogive  opens  up,  use  bottom  solution  1 

'P*-'P*RE',  'AP'-'AFTE', 

•DKFt' 

Ogive  opens  down*  use  fore  solution  | 

'T**-*r*PS’,  and  'B*'-'B*TS' 

'DNAF' 

Ogive  opens  down,  use  aft  solution 

1 

•dot*' 

Ogive  opens  down,  use  top  solution  ^ 

and  use  Table  4. 

*Dl)B«' 

Ogive  opens  dom,  use  bottom  solution 

^For  the  ICASE(l)  value,  one  of  the  Indicated  DIR(I)  values  must  be  chosen.  See  Table  4 for  clarification. 
Except  for  the  noted  values,  DIR(I)  refers  to  the  aft  intersection  of  the  Ith  segment. 
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Table  4.  Nomenclature  for  Input  Variable  DIR(I) 


Table  5.  Jobstream  for  Sample  Case 


-o 


/•PRIORITY  5 

//VKF05333  JOB  (XXXt 

//  VRV06l46t01tPROJ-NeR)t *00000  USER  NAME*? 

//  MSOLEvEL«U»0)  t 
//  RE0l0N»160Kt 
//  CLASS»Xf 
//  TIME*(tl5) 

//  EXEC  FORTEPOSfPGMNOaVRVOOiee 

//GO.FTOlFOOl  00  UNlTaWORK»OISP« (NEW»PASS*0E1.ETEI tOSN»AAKWIKINt 
//  0CBa(RECFN«FBtLRECL»80tBLKSIZE*a00)tSPACE-(CYLtl) 

//GO.FTOSFOOl  OD  • 

hypothetical  asymmetrically  ABLATEOf  SLICEOt  BENT  CONE 
AINPUT 

TYPEa' ASYM* t IPRINT«1 » 

CURYC«2**06IV*(2**LlNE**DIR<2)«>*AFTC*t 

H(l)a.5tK(n«0.*R(l)«.5*TH(n«45.fH(2>«-.6490StK(2)al,14905»R(2)al.I25f 

X<3)«4,92*Z(3)«l,B4>TH(3)*I0.tRH0(2)».5*X(4)>9.22tZ(4)«2.92218fTHU)>l5.t 

AENO 

AINPUT 

CURVE>2*'00IV*f2**LlNE**DIR<2l«*AFTE*t 

H(l)B.37S»K(l)«0.tRCn«.375tTH(l)«30«tH(2)««U812S«K(2)>3,768B6tR(2)a4«t 

X<3)«S.0d242tZ(3)aI«a8fTH<3)«10*fRHO(2>««StX(4)«9*22tZ(4)a2.24i99tTH(4)«5«t 

AENO 

AINPUT 

CUR VEa2* • OG I V * 1 3*  * L 1 NE  * f 0 1 Ra « OPON  * * • UP AF  * , 
H(l)a,437S,K(naO.«R(na,4375tTHina36.325»SLOPE(2)a>T«* 
R(2)a2*5625rX(3)a4»99909t2(3)al.8253fRHO(2)a,5tX(4)a5,00333fZ(4)al,a253« 
X(5)s9«22tZ(S)a2.57l6S**TH(3)al0.* 

XSlOa0.»*eS$35t4,99909f 5»00333t9.22tZSIDa0* t*02666t*.02666«-. 06863* .SOOZA* 
XMAPAXaO**9«22fZMAPAXa0«»0.* 

XU$La6.e6t9.6S»ZUSLa3. 1 I t 1 .83* 

XLSLa3.79*7.97t9.94*ZLSL«1.7S*I.4*.7*RH0LS(2>«4.* 

XSSLa0**9.22*YSSLa2.3*2.3* 

AENO 

//  EXEC  FORTEPOS*PGMNOaVRV00O99(CONOa(8»LE) 

//GO.FTOSFOOl  00  0SNaL4KWIKlN*0ISPa<0LD»DELCTC) 

//GO.FTOTFoOI  do  UNlTB3330*VOLaSERa666666*OSNaUCA671l5* 

//  OCBa(RECENaU*BLKSIZEa80)tDISPa(OLO*KEEPI 
/* 
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Table  6.  KWI KNOSE  Printout  for  Sample  Case 


w 

00 


INPUT  CAROS 

IHYPOTHETICAL  ASYNHETRICALLT  ABLATEOt  SLlCEOf  BENT  CONE 
I lINPUT  • ' 

I TYPCa'ASTM'dPRINTalt 
I CURVC«2*'06tV' (Byline* iOIR(2)*>AFTC** 

I N(l)*.5tKtll>0»tR(l>>«S«THa>«4S.fNI2>»<A4«0S»KI2|Bl«l*905iR«2l«lal2S« 

I XO)44>92<Z(3|4t«a4>tHO>al0.iRHOI2l>aS>AI4)a9.22<ZI4)a2.9221B«TH<4)>lS.i 
• ( iENO  . - ■ 

I 4JNPUT  , . - 

I CURVCa2a<06IY'>2a'LINC'>0IR<2|a*AFTE«<  . 

» HU|a.3TStNll>ao,,R|Ua,3TS.THIl>a30.,N(2)a-l,8125.K(2|a3,TB886.R(2)a4., 

I X(3|aS.oa242fZ(3la|.e8tTHI3>a|0.iRHOI2la.5iAI4|a9.22tZ«4|a2.24t99<THI4ta9.« 

t MNO 
I 41NPUT 

I CURVEa2a'0«IY»i3a*LlNE**01Ra*0PI>N*»»UPAF*t 
I H(l)a,43T9«K(l)a0..R(l)a.43T5tTHma36,32S.Sl.OPE«2|a<T*t 

I RI2|a2.S42StX(3>a4.999a9«ZC3|al.e253iRHOI2)a.S«X<4»a5.00333tZI4)ai.a2S3« 

I XIS|a9.Z2tZI5|a2.STiaStiTH(3|al0.t  . ’ 

I XS10ao.i>a5535a4a99909»S.00333t9.22fZS10aO.«i02ea6«.<l2a64t-.OAa65tt3002ai 
I XNAPAXae«>9.22tZNAPAXa0<«0a> 

I XUSLa6aa6a9.ASaZU$La3atltlaa3a 

I |ll$Ca3aTBaT.9Ta9a94aZLSLal.T9tta4taTaRH0I.S(2|a4.t 
I XSSLa0a«9.22aYSS4a2.3t2a3t 
• 4EN0 


INPUT  data  I 1 ZUC4 


I 

X 

z 

TH 

H 

K 

1 

9.99999C 

69 

9a99999E 

69 

4,S0000C 

o” 

5aOOOOOE«01 

0*0 

2 

9a99999E 

69 

9a99999C 

69 

9.99999e 

69 

-6.490SOE-91 

1 .149052 

00 

3 

4.92000E 

00 

laSAOOQC 

00 

l.OOOOOE 

01 

9aU9999E 

69 

9.99999E 

69 

4 

9aZZOOOE 

00 

2.92218E 

00 

l.SOOOOE 

01 

9.99999C 

69 

9«99999e 

69 

9 

9a99999C 

69 

9a99999E 

69 

9.99999E 

69 

9a99999C 

69 

9*99999e 

6 

9a99999E 

69 

9a»9999e 

69 

9.99999C 

69 

9a9999*E 

69 

9.999992 

69 

7 

9a99999E 

69 

9a99999E 

69 

9.99999E 

69 

9a99999C 

69 

9.99999E 

69 

8 

9a99999E 

69 

9a99999C 

69 

9.99999E 

69 

9a99999C 

69 

9.99999C 

69 

8 

9a99999C 

69 

9a99999E 

69 

9.99999E 

69 

9a99999E 

69 

9.99999E 

69 

10 

9a99999C 

69 

9a99999E 

69 

9.99999E 

69 

9a99999C 

69 

9.99999E 

69 

aiNTUTt 

INa  ItUa  6aIFa  lalPa  TalPRIHTa 

4CN0 


ft 

e 

RHO 

CURVE  DIR  SLOPE 

5.00090€*01 

9,99999E 

69 

9.99999C 

69 

06tV  OROM 

I.I2500C 

00 

9,99999C 

69 

S.00900C*D1 

OOlV  6FTE 

9a99999C 

69 

9,99999C 

69 

9.99999e 

69 

LINE  OPON 

9.99999C 

69 

9,99999C 

69 

9.999992 

69 

LINE  OPON 

9a99999C 

69 

9.999992 

69 

9.999992 

69 

OPON 

9.99999E 

69 

9,94999E 

69 

9.999992 

69 

OPON 

9.99999E 

9.99999E 

69 

9.999992 

69 

OPON 

9a99999C 

69 

9.99999E 

69 

9. 999992 

69 - 

- OPON 

9a99999E 

69 

9.99999E 

69 

9.999992 

69 

OPON 

9a99999C 

69 

9.99999E 

69 

9.999992 

69 

OPON 

laERRNAXa  aa499999aaE«04 
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Table  6.  Continued 


MIMSER  OF  SE6MENT5  NSCO  ■ 4 

ICASEI  1)  ■ 13 
IC4$EC  2)  • $ 

ICASEI  3)  • 14 
ICASEI  4)  ■ 14 
ICASEI  5)  ■ 0 
ICASEI  6>  • 0 
ICASEI  T>  ■ 0 
ICASEI  61  ■ 0 
ICASEI  61  ■ 0 

ICASEIlOl  ■ 0 

SEQUENT  END  POINT  COOIIOINATES  ANO  SLOPES  I 

* I TA  TF 

1 0*0  0.0  4.4999oae  oi  s.ooosoof  oi 

2 1.46446TE-01  3.S35S3SE-01  S.S240T2E  01  4.SOOOOOE  01 

3 4.T20T13E>01  I.OSSTIOE  00  O.OOOOOTE  DO  9.0M40TE  00 

4 4.919999C  00  l.a39999E  00  1,41262SC  01  1.412a2SE  01 

5 9.2I9999E  00  2.9221T9E  00  9.999994C  69  9.999994E  69 

6 9.999994E  69  9,990994E  69  9.499994E  69  9.999994C  69 

T 9.999994E  69  9.999994E  69  9.999994E  69  9.999994C  69 

a 9.999994E  69  9.990994E  69  9.999994C  69  9.909994E  69 

9 9.999994E  69  9.999994C  69  9.999994E  69  9.999994E  69 

10  9,999994E  69  9.999994C  69  9.999994e  69  9,999994E  69 

11  9.999994E  69  9.999994E  69  9,999904E  69  9.999994E  69 

sequent  EQUATIONS! 


A H 

1 1.46446Te*01  S.OOOOOOE'Ol 

2 4.T20T13C>01  -a.AOOSOOE-Ol 

3 4.919999E  00  9.990994E  69 

4 9.219999E  00  9.99<994E  69 

5 9.999994E  69  9.999994E  69 

6 9.999994E  69  9,990994C  69 

7 9.999994E  69  9.990094C  69 

a 9.999994E  69  9,99<494C  69 

9 9.999994E  69  9.990994E  69 

10  9.994994E  69  9.9909942  69 


0*0 

9«OOOOOOPM»l 

9*9999946 

69 

4*5000006 

01 

9*9999946 

69 

)*|490S0€ 

00 

lel26000C 

00 

9*9999946 

69 

9*9999946 

9*9999946 

69 

69 

9, 9999946 

69 

I.TA326aE-01 

1*0000006 

01 

9.72AT1SE-0I 

9«999994e 

69 

9*9999946 

69 

9*9999946 

69 

1*5000006 

01 

4.Sia007C->OI 

9»999994e 

69 

9*9999946 

69 

9 *9999946 

69 

9*9999946 

69 

9*9999946 

69 

9.999994e 

69 

9*9999946 

69 

9*999994£ 

69 

9*9999946 

69 

9*9999946 

69 

9.999994C 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9*999994E 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9a999994E 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9«999994E 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

9*9999946 

69 

FILLET/ROUNO  EQUATIONS  ANO  INTERSECTIONS! 


NSI  ETA  RHO  XI  21 


1 1.4644TE>01  9.9999'<E  69  9,99994E  69  9.99999E  69  9.99999E  69  9.99999F  69 

2 4.T20TIE'-01  a.OT069E-OI  6,213I3E-0I  S.OOOOOE-01  4.|4969E>0I  T.a3694E>01 

3 4.92000C  00  4.99999C  69  9.99499E  69  9.99999e  69  9.99999C  69  9.94999E  69 

4 9.22000E  00  9.99949E  69  9.99999E  69  9.99999C  69  9.99999C  69  9.999992  69 

5 9.99999E  60  9.99999E  69  9.99999E  69  9.99999C  69  9.99999C  69  9.94949E  69 

6 9.99999E  64  9.99999C  69  9,99999E  69  9.94999C  69  9,99999E  69  9.99999C  69 

7 9.99999E  69  9.94999E  69  9.99999E  69  9.99999E  69  9.99999C  69  9.99999E  69 

a 9.99999C  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999C  69 

9 9.99999C  69  9.99999E  69  9.99999C  69  9.99999E  69  9.99999C  69  9.99999E  69 

10  9.99999E  69  9.99999C  69  9.99999E  69  9.99999C  69  9.99999E  69  9.99999E  69 


ri  K2  22  T2 


9*999996 

69 

9*999996 

69 

9*999996 

69 

9*999996 

69 

7*104666 

01 

O.OIOASE'OI 

1*113726 

00 

1*000006 

01 

m 

9*999996 

69 

0.909V4E 

69 

9*999996 

69 

9*999996 

69 

Q 

9*999996 

69 

0.99999E 

69 

9*999996 

69 

9*999996 

69 

9*999996 

69 

9.99999C 

49 

9*999996 

69 

9*999996 

69 

9*999996 

69 

9.99999E 

69 

9*999996 

69 

9*999996 

69 

H 

9*999996 

69 

9.99999E 

69 

9*999996 

69 

9*999996 

69 

X 

9*999996 

69 

9.99999E 

69 

9,999996 

69 

9*999996 

69 

9*999996 

69 

9.99999E 

69 

9*999996 

69 

9*999996 

69 

Nj 

9«99999E 

69 

0.99999C 

69 

9 999995 

69 

9,99999c 

69 

do 

<D 

FINAL  QUICK  ARKAVt 

X 2 

TA 

Table  6.  Continued 

TF 

1 

t 

1 

1 

0.0 

0.0 

4«499998e 

01 

6e999699E 

01 

2 

1.464467E*01 

3.S35S3SC>01 

7.1048S4C 

01 

4.600000E 

3 

4.149692e-»9) 

7.634936C>01 

l»99000oe 

01 

T«io«6a4e 

91 

4 

0.0I0449E-01 

l.ll3717e  00 

9.999997C 

00 

I.OOOOOOE 

01 

s 

4.9I9999C  00 

1.S39999C  00 

1.41262SC 

01 

1.41262SC 

01 

4 

9.219499C  00 

2.922179C  00 

9«999994e 

69 

9«999994C 

69 

7 

9.999994C  69 

9.999994E  69 

9.999994E 

69 

9.994994C 

69 

S 

9.999994C  69 

9»Q99994C  69 

9.999994C 

69 

9.999994C 

9 

9.999994E  69 

9.999994C  69 

9.999994C 

69 

9.499994E 

10 

9.999994E  69 

9.999994C  69 

9*999994E 

69 

9*999994e 

69 

INPUT  data  I 2 ZLCk 


1 

X 

Z 

TH 

H 

K 

R 

8 

RHO 

CURVE  OIR  OLORE 

1 

9.99999C 

49 

9.99999C 

69 

3.00000E 

01 

Z.TSOOOC'.Ol 

0*0 

3.T5000POI 

6.99999C 

69 

9*99999C 

69 

061V  OPON 

2 

9.99999e 

49 

9.99999C 

69 

9.99999C 

69 

•I.81290C 

00 

3*7oae6e 

00 

4*00000E 

00 

9,99999C 

69 

S.OOOOOC-OI 

06IV  APTC 

3 

5e06e42C 

00 

1.08000E 

00 

1 .OOOOOC 

01 

9.99999E 

69 

9«99999C 

69 

9*99999E 

69 

9.99944C 

69 

9.99999C 

69 

LINE  ORON 

4 

9e22000E 

00 

2.24I99E 

00 

5.00000C 

00 

ft  ftftftft^^ 

69 

9«99999E 

9«99999E 

69 

9.99999C 

69 

9.99999C 

69 

LINE  ORON 

S 

9*999996 

49 

9.99999C 

69 

9.99999C 

69 

9.99999E 

69 

9*99999E 

69 

9 .999996 

69 

9.99999C 

69 

9.99999C 

69 

OROR 

6 

9*99999E 

49 

9.94999C 

69 

9.99999E 

69 

9.99999C 

69 

9*999996 

69 

9.99999E 

69 

9.999998 

69 

9.99999E 

69 

' OROR 

T 

9*99999E 

49 

9.999998 

69 

9.99999C 

49 

9.99999C 

69 

9*99999E 

69 

9*999996 

69 

9.99999C 

69 

9.99999C 

69 

OROR 

8 

9*99999E.  49 

9.999998 

69 

9.99999C 

49 

9.99999C 

69 

9*99999E 

69 

9 *999996 

69 

9.99999C 

69 

9.999998 

69 

OROR 

9 

9 *999996 

49 

9.99999E 

69 

9.99999C 

69 

9.99999E 

69 

9*99999E 

69 

9 *999996 

69 

9.99999C 

69 

9.99999C 

69 

ORDN 

10 

9*99999C 

49 

9.99999C 

69 

9.99999C 

69 

9.99999E 

69 

9*99999E 

69 

9*999996 

69 

9.99999E 

69 

9.99999C 

69 

OROR 

AINPUTl 

IRa 

l.lva 

S.ipa 

l.lPa 

T.IPRINTa 

l.CmWAXa  0.49999988C«04 

ACNO 


NUHBCP  or  seoMCNTs  Nscf  ■ a 

ICAMI  1)  ■ 13 
ICASCI  21  ■ S 
KAMI  31  > 1* 

ICASCI  *i  • 1* 

ICASCI  SI  ■ A 
ICASCI  Al  a 0 
ICASCI  71  a 0 
ICASCI  SI  a 0 
ICASCI  «l  a 0 
ICASCIIOI  a 0 


SeOMCNT  CNO  POINT  COOROINATCS  ANO  SLOPCSI 

X Z TA  TF 


1 

2 

3 

4 

s 

A 

T 


0.4 

l.8TS001E*OI 
1.23062TC  00 
S.OS241«e  00 
4.210990E  00 
0.999094E  60 
4.99S994C  SO 


0.0 

3.24TS9SC-01 
1.20223SE  00 
1,ST«999E  00 
2.2419S9C  00 
9.94«994C  69 
9.999994E  69 


2:99999SC  01 
4.9T100SE  01 
9.99999TE  00 
4.9999S0C  00 
9.994994E  69 
9.994994C  69 
9.999994E  69 


e.999S99F  01 
3.00000ZE  01 
9.999997C  00 
4.9909S0E  00 
9.994994E  69 
9.999994E  69 
9.994994E  69 
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Table  6.  Continued 


0 

69 

9.999 994C 

69 

9*999994E 

69 

9*999994E 

69 

9 

«.9««9«*E 

69 

9«999994C 

69 

9*999994E 

69 

9*999994C 

69 

10 

Q.9999MC 

69 

9«999994C 

69 

9*999994E 

69 

9*9999962 

69 

11 

9.99Q9MC 

69 

9*999 994E 

69 

9«999996E 

69 

9*9999942 

69 

SCONENT  EQUATIONS: 

% 

N 

K 

R 

M 

T 

e 

1 

l.0TSa«lE-01 

a.TsooeoE-oi 

0*0 

3.T50000F-01 

9*9999942 

69 

3.000000E 

01 

9.9999942 

69 

2 

l.ZSQQZTE 

oo 

•l,ai2SOOE 

00 

J*76e659€ 

00 

4*ooooooe 

00 

9.909994E 

69 

9*9999«4E 

69 

9.999994E 

69 

3 

S.0a2419C 

00 

9.999994E 

69 

9*999994e 

69 

9.9999942 

69 

l.T63Z»Se-01 

UOOOOOOE 

01 

9.O3a324C-0l 

4 

9.Z19999E 

00 

9.999994E 

69 

9 *9999942 

69 

9.999994C 

69 

9*999994E 

69 

s.ooooooe 

00 

1 .4393452 

OO 

9 

9.999994E 

69 

9, 9999942 

69 

9*9999942 

69 

9.999994C 

69 

9*999994£ 

69 

9*999994E 

69 

9.9999942 

69 

6 

9.9999ME 

69 

9.999994e 

69 

9*999994E 

69 

9*999994E 

69 

9*999994E 

69 

9*9999946 

69 

9.999994E 

69 

7 

9.99999Ae 

69 

9*999994E 

69 

9*999994E 

69 

9 *9999942 

69 

9*999994E 

69 

9*9999942 

69 

9 .9999946 

69 

0 

9.99999AE 

69 

9«999994E 

69 

9*999994E 

69 

9*9999942 

69 

9 *9999942 

69 

9 *9999946 

69 

9 .9999946 

69 

9 

9<999994E 

69 

9«999994E 

69 

9*999994E 

69 

9*9999941 

69 

9.999994E 

69 

9 *9999946 

69 

9.999994E 

69 

10 

9a99999*E 

69 

9q999994E 

69 

9*9999942 

69 

9*9999942 

69 

9 *9999942 

69 

9*9999942 

69 

9.999994C 

69 

FILLET/ROUW  EguATIOMS  AND  IMTERSECTIOMSI 


**•  ET*  »HO  XI  II  n XZ  Zt  TZ 


l*07SOOE*01 

9«94999E 

69 

9.99999e 

69 

9.99999E 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

U23063E 

OO 

1*497392 

00 

T.*eu9E-oi 

5*00000E< 

■01 

1.1Z96ZE 

00 

1*070902 

00 

A.T35Z1E 

01 

I.AIOSTE 

00 

1.Z3ZSSC 

00 

1*000002 

01 

5*002422 

00 

9*999992 

69 

9*999992 

69 

9.99999E 

69 

9.99999E 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

09 

9*22000E 

00 

9*9999yE 

69 

9*999992 

69 

9.99999E 

69 

9*999992 

69 

9*999992 

69 

9.999992 

09 

9*999992 

09 

9*999992 

69 

9.999992 

09 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9.Q9999E 

69 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9.9WQE 

69 

9*999992 

09 

9*999992 

09 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

69 

9*999992 

69 

9. 99  WE 

69 

9.Q9999E 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9.99WE 

69 

9*999992 

09 

9*999992 

09 

9*999992 

69 

9*94«9«£ 

69 

9.99999E 

69 

9*999992 

09 

9*999992 

69 

9« 999992 

69 

9*999992 

69 

9*999992 

69 

9*999992 

09 

9*999992 

09 

9*999992 

69 

9*999992 

69 

9.999992 

69 

9.99999C 

49 

9>99999E 

69 

9*999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

09 

9*999992 

09 

einxl  Quick  array: 


i 


T»  Tf 


1 QeO 

2 Ue7»00IE-01 

3 l«l29624e  00 

4 U410S66E  00 

0 S*0O24l9e  00 

6 9«219999C  00 

7 9.999994E  69 

e 9.999994E  69 
9 9«999994E  69 

10  9*999994E  69 


0.0 

3.247$96e-01 
l«076a96€  00 
1*232&S3E  00 
l,aT9999e  00 
2«24l909e  00 
9«999994E  69 
9«999994E  69 
9.999994C  69 
9«999994E  69 


2*999998E  01 
4.T35210E  01 
leOOOOOOE  01 
9,999997E  00 
4 *9999002  00 
9*999994€  69 
9*9999942  69 
9*999994e  69 
9 *9999942  69 
9.999994E  69 


0*999699^:  01 
3*000002E  01 
4«73S21QC  01 
1 *0000006  01 
4*9999002  00 
9*9999942  69 
9.999994E  69 
9*999994E  69 
9*9999946  69 
9*999994E  69 


INPUT  06TAI  3 YSIO 


I*  2TH  H K R ORHO 


CU9VC  OIR  SLOPE 


9*999992 

69 

9*999992 

09 

3*632502 

01 

A.aTSOOE-DI 

0*0 

4.37500«^-01 

9 •999996 

69 

9.999992 

69 

9*999992 

09 

9*999992 

69 

9*999992 

09 

9*999992 

09 

9*999992 

09 

2*502502  00 

9.999992 

69 

5.000002*01 

4*999092 

00 

1*825302 

00 

1*000002 

01 

9.999992 

09 

9*999992 

09 

9.999992  69 

9*999992 

69 

9*99999C 

69 

5*003332 

oo 

1*825302 

00 

9*999992 

69 

9.999992 

69 

9*999992 

09 

9*999992  69 

9.999992 

69 

9.999992 

69 

9*220002 

00 

2*571052 

00 

9.999992 

09 

9.999992 

09 

9*999992 

09 

9*999992  69 

9.999992 

69 

9*99999C 

69 

OGIV  OPDN 
OGIV  UPAF 
LINE  OPDN 
LINE  OPDN 
LINE  OPDN 
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6 9.99999E  69  0<0  9.99999E  69  9i99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  OPON 

7 9a99999E  69  9«99999E  69  9.99999E  69  9a99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  OPON 

■ 9.99999E  69  9«99999E  69  9,99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  OPON 

9 9<99999E  69  9«99999E  69  9«99999E  69  9.99999E  69  9a99999E  69  9a99999E  69  9.99999E  69  9.99999E  69  OPON 

10  9.99999E  69  9.99909E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  9.99999E  69  OPON 


lUSL  X-SL  Z’SI.  T>SL  RHO>S 


1 

6.S6000E 

00 

SallOOOE 

00 

9.99999E 

69 

9.99999E 

69 

2 

9.65000E 

00 

t.SSOOOE 

00 

9.99999E 

69 

9.99999E 

69 

3 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

4 

9<99999E 

69 

9>99999E 

69 

9>99999E 

69 

9a99999E 

69 

9 

9.99999E 

69 

9«99999E 

69 

9.99999E 

69 

9.99999E 

69 

6 

9.99999E 

69 

9«99999E 

69 

9,99999E 

69 

9.99999E 

69 

7 

9.99999E 

69 

9«99999E 

69 

9.99999E 

69 

9.99999E 

69 

S 

9.99999E 

69 

9>99999E 

69 

9.99999E 

69 

9.99999E 

69 

• 

9.99999E 

69 

9«99999E 

69 

9.99999E 

69 

9.99999E 

69 

10 

9.99999E 

69 

9«99999E 

69 

9.99999E 

69 

9.99999E 

69 

ZLSL  X-M. 

Z-SL 

T-SL 

RHO-5 

1 

3.T90««E 

00 

leTSOOOE 

00 

•.90999E 

69 

9e99999E 

69 

z 

7.9T0(0E 

00 

le40000E 

00 

9.99999E 

69 

4«OOOOOE 

00 

3 

9.94000E 

00 

7eOO0OOE-Ol 

9,99999E 

69 

9e99999E 

69 

* 

9.99999E 

69 

9e99999E 

69 

9.99999E 

69 

9«99999E 

69 

5 

9.99999E 

69 

9*99999E 

69 

9.99999E 

69 

9t99999E 

69 

to 

6 

9.99999E 

69 

9*99999E 

69 

9.99999E 

69 

9.99999E 

69 

7 

9.99999E 

69 

9 •999992 

69 

9.99999E 

69 

9e99999E 

69 

8 

9«99999C 

69 

9.99099E 

69 

9.99999E 

69 

9.99999E 

69 

9 

««99999C 

69 

9.9999 9E 

69 

9.999998 

69 

9.99999E 

69 

10 

9.99999C 

69 

9.99999E 

69 

9.999998 

69 

9.99999E 

69 

VSSL  X*SL 

Z-SL 

T-SL 

RHO>S 

1 

0.0 

2.3O000E 

00 

9.999998 

69 

9.99999E 

69 

2 

9.22000E 

00 

2.30000E 

00 

9.99999E 

69 

9.99999C 

69 

3 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

4 

9.99999E 

69 

9e99999E 

69 

9.99999E 

69 

9.99999E 

69 

s 

9.99999E 

69 

9.99999E 

69 

9 .999998 

69 

9.99999E 

69 

6 

9 #999998 

69 

9.99999E 

69 

9.99999E 

69 

9.99999e 

69 

7 

9 #999998 

69 

9 .999990 

69 

9.99999E 

69 

9.99999C 

69 

8 

9*999998 

69 

9.9999 9E 

69 

9.99999E 

69 

9.99999E 

69 

9 

9#999998 

69 

9e99999E 

69 

9.99999E 

69 

9.99999E 

69 

10 

9.999998 

69 

9e99999E 

69 

9.99999E 

69 

9.99999E 

69 

ZNAPAX 

XSID 

ZSID 

1 

0.0 

0.0 

0.0 

0.0 

2 

9.22000E 

00 

0.0 

8.55350E-01 

2.e6600E-02 

3 

9.999998 

69 

9.99999E 

69 

4.99909E 

00 

2.066000-02 

4 

0.09990E 

69 

9.999998 

69 

5.00333E 

oo 

-6.e6*99E-02 

8 

9 #999998 

69 

9 .999998 

69 

9.22000E 

00 

3.00260E-01 

6 

9 #999998 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

7 

9.99999E 

09 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

• 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

9 

9 #999998 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 

10 

9.999998 

69 

9.99999E 

69 

9.99999E 

69 

9.99999E 

69 
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4INPUT1 

IP*  IflWa  btIFa  lilPa  TtlPRINTa  liERRNAXa  0.49M998SE-04 

UNO 


NUMBER  OF  SEONENTS  NSEO  a 5 


ICASEI  II  a 13 
ICASEI  21  a 19 
ICASEC  31  a 14 
ICASEI  41  a 6 
ICASEI  51  a A 
ICASEI  61  a 0 
ICASEI  Tl  a 9 
ICASEI  Bl  a 0 
ICASEI  91  a 0 
ICASEI 101  a 0 


SE6MENT  END  POINT  COORDINATES  AND  SLOPCSI 


X 

2 

TA 

TF 

1 

OeO 

0.0 

3.632498E 

01 

8.999899s 

01 

2 

U7B3406E-01 

3.S24806E-01 

S.a93269E 

01 

3.63249BE 

01 

3 

6.5S347SE*01 

1.094646C 

00 

9t999994E 

00 

9.999994E 

00 

4 

4.999069E 

00 

1.B25299E 

00 

OiO 

0.0 

S 

5.003329E 

00 

1.82S299E 

00 

1 .0037398 

01 

1.003739E 

01 

6 

9.219999E 

00 

2.S7I6S0C 

00 

9.999994E 

69 

9 .9999946 

69 

7 

9.999994E 

69 

9,999994E 

69 

9.9999946 

69 

9.999994E 

69 

8 

9.999994E 

69 

9.999994E 

69 

9.999994E 

69 

9.999994E 

69 

9 

9.999994E 

69 

9.999994E 

69 

9.999994E 

69 

9.9999946 

69 

10 

9.999994E 

69 

9.999V94E 

69 

9.999994E 

69 

9.999994E 

69 

11 

9.999994E 

69 

9.999994E 

69 

9.999994E 

69 

9.999994E 

69 

SEOHENT  EQUATIONS! 

X 

H 

K 

R 

H 

T 

B 

1 

|.T83406E*01 

A.3TS000E-01 

0.0 

4.3750OQa-01 

9.9999946 

69 

3.6325006 

01 

9.9999946 

69 

2 

8.5S34T5E-01 

*1 .3398936 

00^.  2.41T009E 

00 

2.562S00C 

00 

69 

9.9999946 

69 

69 

3 

00 

9.999994E 

69 

9.9999946 

69 

69 

l.T63268E*01 

1.0000006 

01 

9.4382586-01 

4 

S.003329E 

00 

9 .9999946 

69 

9 .9999946 

69 

9.Q9999AE 

69 

0.0 

0.0 

1.82S209E 

00 

5 

9.219999E 

00 

9.9999946 

69 

9.9999946 

69 

69 

1.7699996*01 

1.0037396 

01 

9.3971066-01 

6 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.9999946 

69 

7 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.999994C 

69 

9.9999946 

69 

9.9999946 

69 

9.9999946 

69 

8 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.9999946 

69 

9 

9.999994E 

69 

9.999994E 

69 

9.9999946 

69 

9.Q99994E 

69 

9.9999946 

69 

9.9999946 

69 

9.9999946 

69 

10 

9.999994E 

69 

9.9999946 

69 

9.9999946 

69 

9.999994C 

69 

9.9999946 

69 

9.9999946 

69 

9.9999946 

69 

FILLET/ROUNO  equations  and  INTERSECTIONS! 


KSI 


ETA 


RHO 


X| 


21 


Tl 


X2 


l.Ta341C-01 
8.5S34TE-01 
4.99909C  00 
S.00333E  00 
9.22000E  00 
9»99999E  69 
9.99999E  69 


9*99999E  69 
1>15439E  00 
9.99999E  69 
9a9999«E  69 
9>99999E  69 
9>99999E  69 
9i99999E  69 


9.99949E  69 
6.39662C-al 
9.99999L  69 
9.99999E  69 
9.99999E  69 
9.99999E  69 
9.99999E  69 


9.99999E  69 
S.OOOOOE-OI 
9t99999E  69 
9.99999E  69 
9.99999C  69 
9 •999991  69 
9.99999E  69 


9.99999E  69 
T.4T208E-01 
V.99999E  69 
9.99999E  69 
9.99949E  69 
9<99999E  69 
9.99999E  69 


9.99999E  69 
9.29842E-0I 
9.99999E  69 
9.99999E  69 
9.99999E  69 
9.99999E  69 
9.99999C  69 


9.99999E  69 
5.45243E  01 
9.99999E  69 
9.99999E  69 
9.99999E  69 
9i99999E  69 
9.99999E  69 


9.99999E  69 
1.06T56E  60 
9.999998  69 
9.99999E  69 
9.99999C  69 
9.99999C  69 
9.99999E  69 


22 

9.99999E  69 
1.1320TE  00 
9,99999E  69 
9.99999E  69 
9.99999E  69 
9.99999C  69 
9.99999E  69 


T2 

9.99999E  69 
l.OOOOOE  01 
9.99999C  69 
9.99999C  69 
9.99999E  69 
9.99999E  69 
9.99999C  69 
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• 

9e00990e  09 

9.999996 

09 

9,999996  69 

9.999996  69 

9.999996  69 

9.999996 

09 

9 

9«99999e  00 

9.499996 

09 

9.999996  69 

9.999996  69 

9.999996  69 

9.999496 

09 

10 

9e99999e  09 

9.999996 

09 

9,999996  69 

9.99*946  69 

9.999996  69 

9.999496 

09 

9.9999W  M 
9.9«999e  M 
*.9«4«4C  M 


9.9999«e  69  9.99999e 
9.99999C  69  9,99999E 
9.999996  69  9.99999C 


9.99999C  69 
9.999996  69 
9.999996  69 


rilUl.  MICK  MWAYI 


X 

2 

T9 

TF 

1 

OeO 

0.0 

3.6324006 

01 

0.9990996  01 

i 

1.7a3406E>01 

3.S24006C>01 

. 5.4524206 

01 

3,6324906  01 

3 

7.472076E>ai 

9.29a42oe>oi 

1 .0000006 

01 

5.4524206  01 

4 

1.0675656  00 

1.1320666 

00 

9.99999*6 

00 

1.0000006  01 

S 

4.9990096  00 

1.0292996 

00 

OtO 

OeO 

0 

5.0033296  00 

1.0252996 

00 

1.0037396 

01 

1.0037396  01 

7 

9,2199996  00 

2.5716506 

00 

9.9999946 

00 

9.9499946  69 

0 

9.9499946  69 

9.9999446 

99 

9.9999946 

60 

9.9999946  69 

9 

4.9999946  69 

9.9999946 

69 

9.9999946 

69 

9.4949946  69 

10 

9.999994C  09 

0.090904C 

99 

0.P0M94E 

00 

0*090004E  09 

SLICE  OEONEYRYI 

UPPER 

LOWER 

X 

z 

T 

X 

2 

1 

OeO 

6.2572306 

00 

•2.4644706 

01 

OeO 

2,0673446 

00 

2 

9.6500006  00 

i.aaeoooE 

00 

•2,4644706  01 

7.9699996  00 

1,400000E 

00 

3 

9.9999946  69 

9.99099*6 

99 

9.9999946 

00 

9e940000E  00 

7.0OOOOOE*OJ 

« 

9.9999946  69 

69 

9.9999946 

00 

9.9999946  69 

9e400994E 

09 

S 

9.9999946  69 

9,999994C 

99 

9.9999946 

00 

9*999994E  09 

««400094E 

09 

0 

9.9999946  69 

9,9999946 

69 

9.9999946 

09 

9.99999*6  69 

9.999994E 

09 

7 

9.9999946  69 

«e999994E 

99 

9.9994946 

09 

9.9999946  69 

9.«999946 

09 

a 

9.9999946  69 

9.990994C 

99 

9.9999946 

09 

9.9999946  69 

9.9999946 

09 

9 

9.9999946  69 

9«999994E 

99 

9.9999946 

09 

9e999994C 

09 

10 

9.9999946  69 

9.9999946 

99 

9.99999*6  69 

9.9999946  69 

9.4999946 

09 

St.IC6  R0UH0/FILL6T  CCNTCRS  AND  INTCRSeCTKMSI 
Bt  aj  AS 


•4.7663316 

•4.7663316 

•1.9M16S6 

9.9999946 

9.9949946 

9.9999946 

9.9999946 

9.9999946 

9.9999946 

9.9999946 


XI 


61 


Sloe 

X » T 


00 

0.0 

2.299999E 

00 

OeO 

00 

9.2199996 

00 

2.2999996 

00 

OeO 

01 

9e999994E 

00 

9.9999946 

09 

9e999994C 

69 

09 

9.99999*6 

09 

9.99999*6 

09 

9.9999946 

09 

09 

9e999994E 

00 

9.9999946 

09 

9.9999946 

09 

69 

9»999994C 

00 

9^000094C 

09 

9.9999946 

09 

09 

9.9999946 

00 

Oe099994E 

09 

9.9999946 

09 

09 

4.9999946 

69 

9«099994C 

09 

9.9999946 

09 

09 

9e999994C 

09 

9.9999946 

09 

9.9999946 

09 

09 

9.9999946 

09 

Oe090004E 

09 

9.99999*6 

09 

71 

KE 

22 

72 

66  *1.946656  96  -1.39367C*02  •2.946766  09  7,119426  99  7.453166  60  1.443276  00  >4.766336  00  6,456706  00  1.626356  00  •1.956166  61 


FINAL  MICK  ARRAYS  FOR  R0UNDeD/FILLE7760  SLICESI 


UPPCR 

X 2 7 


X 


LOWER  0106 

2 7 X Y 7 


OeO 

6.2572366 

00 

•2.4644766 

01 

0.0 

9.6500006 

00 

la630000E 

00 

-Se4044T«E 

01 

7.4531746 

00 

9.9999946 

00 

9,9999946 

09 

9.9994946 

09 

S.450690E 

00 

9.99999*6 

00 

9.9999946 

09 

9.9999946 

69 

9.9400006 

00 

9.9999946 

09 

9.9999946 

09 

9.9999946 

09 

9.9999946 

00 

9.9999946  69 

9.999994E 

69 

4.9999946 

09 

9.9999946 

69 

9.9999946 

09 

9.999990C 

09 

9e940994E 

69 

9e999994e 

69 

9.9999946 

09 

9,99999*6 

09 

9e999994E 

09 

9,99999*6 

69 

9.99999*6 

09 

9.9999946 

69 

9.9999946 

69 

9e999994e 

09 

9.99999*6 

09 

0e099094E 

09 

9.9999946 

09 

9.9999946 

69 

2.0673446  00  •4.7663316  DO 
1.4436746  00  •4.7063316  00 
1.9561656  01 
1.9961656  01 


1.2663516  00 
7.900000e<01 

9.9999946  69 

9.9999946  69 

9.9999946  69 
9.9994946  69 

9.9999946  69 

9.9999946  69 


9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  64 


0.0 

9.2199996  00 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

4.9999946  69 

9.9999946  69 


2.2999996  00 
2.2999946  00 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 


0.0 

6.0 

9.9999946  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 
9.9999446  69 

9.9999946  69 

9.9999946  69 

9.9999946  69 
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t/i 


IHVPOTNETICAL  *STHMCTnlC*U.V  Aai.ATeO«  SLICCDt  SENT  CONE 
t I 

I I 6 nomcircular  cross  section 

ILCLL  ICLLl  PlEC  BLCL  SIO  LSCP 

lUELk  2CLLI  PJEC  SID  TUCL  USCP 

ILBSL  SLINE  PIEC  BLSL  LSkCP 

ILSSL  AklHE  PIEC  LSLCP  $$(.  (.ELL 

lUSSL  SLINE  PIEC  SSL  U5LCP 

lUTSL  SLINE  PIEC  USLCP  TUSL  UCLL 

( I MAP 

ill  0.0  S.tZOOO 

• VUCL 

I 1 LINE  PIEC  KV5 

I 0.0  0,0  0.22000  0.0  A 0.0  A 0.0 

1-1 

IZUCL 

I 1 ELLA  PIEC  KVO 

I 0.0  0.0  O.IASAS  0.3S3SSA  S9.99S90A  AA.0999S 

I 2 ELLA  PIEC  KVO 

I O.USAS  0.353S5  0.41A9T  0.TB3S9A  tS.OOOOOA  Tl. 0*0(10 

I 3 ELLA  PIEC  KVO 

I 0.41*97  0.703*9  0.00104  1.113T2A  T1.04SS4*  10.00000 

I * LINE  PIEC  KVS 

I 0.00104  1.11372  4.92000  1.04000*  10.00000*  10.00000 

I S line  PIEC  KVS 

I 4.92000  1.04000  9.22000  2.92218*  14.12025*  14.12025 

1-1 

I2LCL 

I I ELLA  PIEC  KVO 

I 0.0  0.0  0.18750  -0.32470A-00.99B99A  -29.99998 

I 2 ELLA  PIEC  KVO 

I 0.1S7S0  -0.32*70  1.12902  •1.0Ta90*-30. 00002*  -47.35210 

I 3 ELLA  PIEC  KVO 

I 1.12902  -1.07*90  1.41057  -1.2325SA-47. 35210*  -10.00000 

I 4 LINE  PIEC  KVS 

I 1.41057  -1.23255  5.082*2  -I.BBOOOA-IO. 00000*  -10.00000 

I 5 LINE  PlEC  KVS 

I 5.08242  -1.88000  9.22000  -2.24199*  -4.99998*  -4.99998 

1-1 

lYSIO 

I I ELLA  PIEC  KVO 

I 0.0  0.0  0.17834  0.35248*  89.99899*  36.32498 

I 2 ELLA  PIEC  KVO 

I 0.17834  0.35248  0.74721  0.92984*  36.32*98*  54.52428 

I 3 ELLA  PIEC  KVO 

> 0.74721  0.92984  1.0*750  1.13207*  54.52*28*  10.00000 

I 4 LINE  PIEC  KVS 

I 1.00750  I. 13207  4.99909  1.82530*  10.00000*  9.99999 

I 5 LINE  PIEC  KVS 

I 4.99909  1.82530  5.00333  1.82530*  0.0  * 0.0 

I 0 LINE  PIEC  KVS 

I 5.00333  1.82530  9.22000  2iS716S*  10.03739*  10.03739 

1-1 

IZUSL 

I 1 line  PIEC  KVS 

I 0.0  0.25724  9,65000  1.83a00*-24.0*470*  -24.64478 

1-1 

IZLSL 

I 1 line  PIEC  KVS 

I 0.0  -2.06734  7.45310  -1.44327*  4.78633*  4.78633 


LELL 

UCLL 
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Table  6.  Concluded 


1 2 CLLX 

PIEC  KVO 

1 T.*s3ia  •i.aisaT 

e *43070  *1*226350 

4.766334 

19.56169 

1 3 LINE 

PIEC  KV5 

1 

1 StiSaTa  <>1.22635 

9.94000  -0.700004 

19.561654 

19.56165 

1-1 

irssL 

P 1 LINE 

’ PI6C  KVS 

1 0.0 

2*30000 

9.22000  2.300004 

0.0  4 

0*0 

1-1 

IZSIO 

P 1 LINE 

PIEC  KV5 

1 0.0 

0.0 

0.05535  0,020664 

1.919084 

1.91906 

P 2 LINE 

PIEC  KVS  ' 

t 0.05535  0.02066 

4.99909  0.026664 

6,0  6 

0*0 

' 1 

t 3 LINE 

PIEC  KV$ 

1 6.90909  0.02066 

9.00333  -0.066654-57.505024 

-67.90902 

1 4 line 

PIEC  KVS 

1 

* 5« 00)33  -0«060«s 

9*22«00  ••300266 

4,90OO«A 

4«99999 

t-i 

tZPUP 

- 

1 

11  line 

PIEC  KVS 

1 

I 0.0 

0.0 

9.22000  0.0  4 

0.0  - 4 

0.0 

1 

irusL 

YUCL 

• V ■ . 1 

irLSL 

YOCL  . 

PTLCL 

YUCL 

« 

PYN6P 

YUCL 

*■  ■ 

ITUSCP 

YSID 

trtscp 

VSIO 

trusLCP 

YSSL 

lYLSLCP 

YSSL 

PZUSCP 

ZUCL 

PZLSCP 

ZLCL 

PZUSLCP 

ZliSt 

1 

PZLSLCP 

ZLSL 

I 

IZSSL 

ZSIO  , 

1 1 2 

1 0 

o.ooioo 

9.22100  0.09210 
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Table  7.  Definitions  of  Logical  Names  for  QUICK 
Cross-Section  and  Body  Line  Models 


Cross-Section  Models 


Logical 

Name 

Definition 

LELL 

Lower 

ellipse 

UELL 

Upper 

ellipse 

UTSL 

Upper 

top  slicing  plane 

USSL 

Upper 

side  slicing  plane 

LSSL 

Lower 

side  slicing  plane 

LBSL 

Lower 

bottom  slicing  plane 

Body  Line  Models 

LCL 

Lower 

centerline 

LSCP 

Lower 

slope  control  point 

MAPAX 

Map  axis 

SID 

Side 

UCL 

Upper 

centerline 

USCP 

Upper 

slope  control  point 

LSLCP  Control  point  for  lower  horizontal  and  vertical 

slicing  planes 

LSL  Centerline  on  lower  slice 

SSL  Side  on  side  slice 

USL  Centerline  on  upper  slice 

USLCP  Control  point  for  upper  horizontal  and  vertical 

slicing  planes 

Note:  See  Fig.  8 also. 
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Tables.  KWIKNOSE  Source  Lbting 

C KWIKNOSE  > PROGRAM  TO  SET  UP  QUICK  INPUT  FOR  AN  AAISNMETRIC  STACK  OF  100 

C CONES  AND  OGIVES/  NOT  NECESSARILY  TANGENT  200 

C ORIGINAL  PROGRAM  VERSION  2/IB/7T  300 

c— ......J... AOO 

c INPUT  Instructions  (SEE  AEDC-TR-77-89  FOR  omiLS)  soo 

C first  INPUT  CARO  IS  10  data  (SOAliQUlCK  READS  ONLY  60A1}  GOO 

C REMAINING  INPUT  IS  VIA  NAMELIST  /INPUT/  700 

C SET  IPRINT  > 1 TO  MAKE  QUICK  EXERCISE  HATH  MODEL  SOO 

C provide  do  CARO  FOR  UNIT  IF  (DEFAULT  w II  WITH  DCBw (RECFMwFBt  900 

C LRCCLwSOI.  1 CYL  LONG  SEQUENTIAL  FILE  1000 

C IF  IP»1F  QUICK  INPUT  DECK  WILL  BE  WRITTEN  ON  IP  WITHOUT  REWINDING.  USE  UOO 

C TO  PUNCH.  1200 

C>— — — . 1300 

COMMON  /RA/XllOI.2llO).TH(10)«H{10>.K(10I.R(lVI.B<10)tMtlOI.  lAOO 

1 CURVE(10).OIRI10).XlllUI.X2UO).Zl(10I.Z2(10ltTl<10>iT2<10).  1500 

2 KSI(10>«CTA(10I.RHO(10I«XQI|1I.ZQI11).TQI11»2I.KRV(3I.CL<2I.  1600 

3 KV(2).VI10t3.3l tSLI3l .HHOSL(10.3ltCRVSL(10t3>  1700 

A.VB(10.3.2.SI  .CRVB  (lOtSI  .XHAPAXllOl  .ZMAPAX  ll»|  .XSIOdOl  .ZSIOdOl  . IBOO 

S BLBIB)  1900 

COMMON  /lA/ICASEIlO) .JCASEI30) tNCASEI7) .NPI3I •NS(3> «NB(ei  2000 

COMMON  /RS/K!.M|tKU«MJ.MIJ.L|NE<HLNKtRlNFiERRMAX.OPONtCLLXt  2100 

1 YUCL.ZERO.FORE.AFTE.TOPStBOTStOPUP.OGIV.PIOlSO.XJtZJiXl.ZI.HIt  2200 

2 RIiBI.THI.BJ.Al.A2.A3.bIJiTHU.S.THIJ.SIJ.AI.t1JrRJ.A*.0X  2300 

3iTYPEtASTM.UPFO.UPAF.UPro.UPBO.ONFO.DNAF.DNTO.ONBO  2*00 

COMMON  /IS/IRtlw.IF.IP.lPHINTtNSEG.NSEGPl.lSTOP. JSTOP.KSTOP.  2500 

4*.  1 l.J.L.lPl.LMItIl«NPTS.12.ISL.LB.N.NM1.13  2600 

OO  COMMON  /LA/SLOPE(10>.IDEN{SO).CARDI80)  2700 

COMMON  /LS/TI.T'  2800 

REAL  KtM.KSl.KltMI.KJ.HJ.MIJ.LlNEtKV.KRV  2900 

LOGICAL*!  SLOPE. lOEN. CARO. Tl.T  3000 

LBwl  3100 

CALL  INIT  3200 

CALL  IHPT  3300 

IFITYPE.EO.  A5YMIG0  TO  100  3*00 

CALL  CASE  3500 

CALL  SEGE  3600 

CALL  INTfi  3700 

CALL  OARY  3800 

CALL  PfiNT  3900 

CALL  ROFL  *000 

CALL  SLCE  *100 

CALL  RFSL  *200 

CALL  QUCK  *300 

CALL  LIST  **00 

STOP  ... 

100  00  120  LB«t.3  *600 

IFILB.EQ.IIGO  to  110  *700 

CALL  INIT  *S00 

CALL  INPTl  *900 

110  CALL  CASE  5000 

CALL  SEGE  5100 

CALL  INTR  S200 

CALL  OARY  5300 

CALL  PRNT  5*00 

CALL  ROFL  5500 

CALL  COPY  5600 
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Table  8.  Continued 


■1^ 

VO 


120  CONTINUE 
CilLL  SLCE 
CitLL  RPSC 
00  130  LfiMtO 
CALL  COPY 
130  CONTINUE 
CALL  BLCK 
CALL  OK30 
CALL  LIST 
STOP 
END 

SUBROUTINE  INIT 


C INITIALIZE  variables 


CONNON  /RA/X(10I  iZ IIO)  tTHI  tOI  •HdOltAtlO)  (R IIOI  .BaO)  tN(10>  • 

1 CURVEI10)<OIR(10liXl  liO>«X2llO>.2l(10>tZ2UOI  fTl  1101  *T2 (101 1 

2 RSI(IO)«ETA(10>.RNO(|OI»XO(lll.Z«llll,TO(llf2l tKRV(3ltCL(2l« 

3 KV(2| iV(10f3«3l •SL(3) tRHOSLIIOtSl tCRVSLI10i3l 

«iVB(10t3>2<SI«CRVB  (lOtSl  (XMAPAX  (101  tZMAPAXdOl  *XS10<10>  t ZStOdOl  • 
5 BLBiai 


CONNON  /I  A/ICASE  (10)  t JCASCdOl  (NCASC  (Tl  •NP(3I  *NS  (31  .NB(a) 

CONNON  /RS/KI«NItKJtNJ>NlUtLlNE*BLNX«RlNP»ERRNAX<OPDNtCLLX» 

1 VUCL*ZeRO«FORCtArTEfTOPSteOTS>OPUP«OetV>PIC180«XU»ZJtXliZlfHl« 

2 RltBUTHI<BU«Al.A2tA3tBlJ«THJ($fTHtJ«SlU»Al>NJtRUtAAtOX 
3tTTPEiASVN<UPFO.UPAF.UPTO<UPBOtONFO.ONAP«ONTO.ONBO 

CONNON  /IS/IR<lH*iFtIPdPRINT*NSEO>NSCOPI«lSTOP«slSTOP«KSTOP» 

1 ltJtLtlPl«LNltIl«NPTSil2*lSL«LB<NtNNl«I3 
CONNON  /LA/SLOPE(IO)»IO£N(80ltCARO(SO> 

CONNON  /LS/TlfT 

real  K«N>KSIfKl»NI.KJ.MJ.NIJ.LINE>KV>Knw 
LOeiCAL*!  SLOPEiIOENtCARDtTIiT 
1000  CONTINUE 

lF(LB.eT.l)00  TO  1100 

1PR1NT«0 

IRbS 

1HB6 

IFbI 

IPbT 

ERRNAXaS.E-S 
1100  NSEBbIO 

PI01S0B3.)«t5926/ia0. 

ISTOPbO 

JSTOPaO 

KSTOPbO 

RETURN 

end 

CB*BBBBB*BB**B»B***#B*B*BBBBB*BBB#BBBBB***BBBB**BB**BBBB**BBBBBBB*#«BB#B 

subroutine  INPT 

C— 

c obtain  input  data 


CONNON  /RA/X(10)  >2(10)  <THdOI  tHdO)  tKIlO)  *R(10I  >BdO>tN(IOI» 

1 CURVE(10)t0IR(10)IXl(l0)tX2(10)>Zld0l  tZ2(10)«TI(10)»T2dO)« 

2 KSI(l0)«ETAd0)tHNO(10)tXa(lll»ZQ(lll  »Te(llt2l«KRV(3)«CL<2)t 


5700 

BBOO 

5900 

0000 

0100 

0200 

6300 

6«00 

0500 

0600 

OTOO 

oaoo 

0900 

7000 

7100 

7200 

7300 

7400 

7500 

7000 

7700 

7000 

7900 

6000 

oieo 

0200 

0300 

0400 

BSOO 

0000 

6700 

0000 

B900 

9000 

9100 

9200 

9300 

9400 

9500 

9000 

9700 

9000 

9900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

loooe 

10900 

11000 

11100 

11200 
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Table  8.  Continued 


o 


3  KVI2I • V (10t3t3) (SL (31 tRHOSL I t0«3t (CRYSL (10t3) 

*«VB(I0(3>2«e)»CRVB  (t0>«>  (XMAPAX (10>  <2HAPMI  10>  iXS10(10>  »Z$ID(IO) • 
S BLBI8> 

COMMON  /IA/lCASEII0lf>fCASC(30>  tNCASe(T)<NP(3ltNS(3>«NB(e> 

COHHON  /flS/KItMI<KJ«MJ»NlJ«LINE»BCNK«RINFtERRMAXfOPON«eLLXt 
I YUCLtZEROtFOREtAFTEfTOPStBOTStOPUPiOGIVtPlOlSOtXJtZJtXltZItHl* 

Z RItBI»THl»ejfAt<A2«A3tBIJtTHJi5>TN|JtSIJ»Al«HJ>RJ<AA<0X 
3*TVPEtASYM«UPFOiUPAFfUPTOtUPBOiONFOtONAF<DN10«ONBO 
- COMMON  /I5/lRtlWf IFtlP.lPRINTrNSEG.NSEOPltlSTOPtJSTOPtXSTOPt 
I t«J,LtIPl«I.MttntNPTS>12>lSL>I.BtN>NM)tI3 
COMMON  /LA/SLOPCIIOt<10eN(80)  iCAROIBOl 
COMMON  /LS/TItT 

REAL  K«MtKSl*Kt«HI tKJ<MJ«N|J,LINEfKV«KRV 
LOGICAL*!  SLOPEtlDEN<CARO»Tt»T 
DIMENSION  - XUSL(10l«ZUSL(10ltTUSL<l0>t 

1 XLSL(10»*ZLSL(10)iTLSL<t0)» 

2 ' XS$L(10)tVSSL(10)tTSSL(l0) 

EGUIVALENCC 

1 (V(l»lt!)tXUSL(l>tt(V<l«2tl)«ZUSL(l))f  (VlU3illtTU$L(ln 
2f (V(lflf2)«XLSL(l)><(Vll«2t2)iZLSL<l)if (V(lt3>2ltTLSL(l)> 
3«(V(I>lt3)«XSSL(l>l t(V<!»2t3>trSSL(l)) t(Vllt3«3>«TSSLIl)! 

DIMENSION  RMOUS(10liRHOLSI!0)«RHOSSllGI 
EQUIVALENCE  (RHOSL (l*ll iRHOUS (1 t ) 

1 «(RH0SL(lt2l«RH0LS(ll> 

2 « (RHOSL (ItSliRHOSS (111 
NAHELIST/INPUT/X»ZtTHfCURVE«RtHiKf  Bt IRtIHt IF »ERRMAX»OtR« IPRINTt 

1 IP,RHOfSLOPE«XUSLtZUSL*TUSL<XLSL<ZLSL<TLSL«XSSLtYSSLtTSSL 
2 t RHOUS t RHOLS • RHOSS t type • XM AP AX • ZMAP AX • XS I 0 t ZS I 0 
NANELlST/INPUTI/IR(lH*lFfIP*lPRINT»ERRMAX 

2000  CONTINUE 
IF(IF.EO.IP)00  TO  2009 
REMIND  IF 

2002  READ(IRfOOO»ltENO*200A)CARO 
HRITE ( IF t 0000 I • ENO>2003 > C ARD 
GO  TO  2002 

2003  STOP  2003 
200A  IR*IF 

REMIND  IF 
2005  CONTINUE 
2010  FORMAT (80AI) 

HRITE(IM<200II 

2001  FORNAT(n««30X.23('<*>MF31Xi23(><*>M/ 

- 3IXt»<*>**43X**«*>«/ 

31Xt«<*>  program  KMIKNOSE 

1 TO  SET  UP  QUICK  INPUT  FOR  NOSE  GEOMETRIES  «•>•/ 

2 31X.»<*>».63Xf»**»*/ 

3 31Xt'<*>  VON  KARMAN  GAS  DYNAMICS  FACILITY  (VKF/ADPI  ••13X* 

4 •<•>•/ 

5 3IXt»«*>t, *3*, #<•>•/ 

6 3IX««<*>  ARNOLD  ENGINEERING  DEVELOPMENT  CENTEPt  TENNESSEE  3T389t 
T.  TX. •«•»•/ 

8 31X»‘«*>**6SX. •<•>*/ 

9 31X #23 I *<*»»1/3|X. 23 (*«•>» >////» 

IFIIF.EO.IPIOO  TO  2013 

REMIND  IF 
HRITE (IH#200«I 


11300 

llAOO 

11500 

11600 

IITOO 

11800 

11900 

12000 

12100 

12200 

12300 

12AOO 

12500 

12600 

12T0O 

12800 

12900 

13000 

13100 

13200 

13300 

13A00 

13500 

13600 

I3TOO 

13800 

13900 

lAOOO 

1*100 

14200 

14300 

14400 

14500 

14600 

14T00 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

15600 

15T00 

15800 

15900' 

16000 

16100 

16200 

16300 

16400 

16500 

16600 

16T0D 

16800 
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Table  8.  Continued 


200ft  F0f)H*T<>  INPUT  CMOS*)  IftOOO 

mlTE<l«i200r)  ITOOO 

200T  FORMOTC  ••S2<'-*))  ITIOO 

2008  RE«D(IF««OOBt>ENO>2011ICARO  IT200 

MR|TE(tM«2O09ICARO  IT300 

2009  FOR)UT<«  l••80•l••l•)  ITftOO 

00  TO  2008  ITSOO 

2011  WR|TE4tW«2O0T>  ITbOO 

2013  CONTINUE  ITTOO 

IFIIF.NE.IPIREWINO  IF  17800 

REAO(lR«2010tENO>2O20)IOEN  17900 

ReAO(|R>lNPUTtENOa20ftO>  18000 

2019  HRITEdMtOOSOI  18100 

WRITE<lNt21IllL8>BLBILB>  • 18200 

2111  FOflHAK*  INPUT  DATA*  •«  I2<2A«A4I  18300 

HR1TE<IW«2112)  18400 

2112  FORMAT«>0  6X«>  A'tllXi*  Z'tllXf •TH>«11X><  HifllXt*  K«tllX>  lOSOO 

1 • R'tllXf*  B'tllXi'RNci*)  18000 

HR1TE<I«i2113I  18700 

2113  FORMAT<«**t)lTXf*CURVE  U|R  SLOPE*)  18000 

HRITC<IH.2114)  18900 

2114  FORMAT)*  ••132(*-*))  19000 

WRITE (IW>21lSl  lI«X(tl«2(t)*TH(ll  «H(I)«K(l)«Rttl  iBIDtRHOdl*  19100 

lCURVC<l>tOIRIt>«$LOPC(ll*  19200 

2 IxltlO)  19300 

2115  FORMATdOl*  • tUt  lP0E13.5«llXf  2 1 lX<  A4)  «4X<A1/) ) 19400 

tF4TYPE«Ee«ASVM.ANOiLeiLT«3IOO  TO  2023  19900 

DO  2118  L<1«3  19000 

WR|Te(IW«21|ft)SL<L)  19700 

2110  FORHAT(*0<tA4»  3Xt*X«SL*t  9X<*Z-SL*«  SXi'T'SL**  8Xt*RHO-S*/  19800 

I • i«54<*«*>)  19900 

write (IWt21 17) II<<V(|tU*L)«J>Ii3>>RH0SLIIiL) fl«l>10)  20000 

2117  FORMATIlOf*  • • 12* 1P4E13*S/) > 20100 

2118  CONTINUE  20200 

WRITCilWt2021)  20300 

2021  FORMAT  I •0*«6X<*XMAPAX*«7X|I2HAPAX*.8X.*XSIO*>9X«*ZSID*  20400 

1/*  •t54(*-*|)  20500 

WRITCIlw«2022)  il.XMAPAXtl)  tZMAPAX  ID  tXSIDd)  *251011)  •l>l«10l  20000 

2022  FORHATIIOI*  •.I2.1P4C13*5/) > 20700 

2023  WRITEIIW.INPUTI)  20800 

WR|TCIIW*90SO)  20900 

GO  TO  3000  : 21000 

2020  WR|TEIlMi2030)  21100 

2030  FORMAT  I*  COMO  CODE  2020 1 MISSING  ID  CARO*)  21200 

NRlTE<Iw*2035)  21300 

2035  format  I tOEXECUTION  HALTED*)  21400 

STOP  2020  21500 

2040  HRITE<1W«2050) IR  21600 

2050  FORMAT  I*  COHO  CODE  20401  EOF  ON  UNIT  IRa'tlZ)  21700 

WRITE! 1W« 2035)  21800 

STOP  2040  21900 

3000  RETURN  22000 

ENTRY  INPTl  22100 

lFIIF.NE.IP)IRaIF  22200 

REAOIIR*INPUT*ENO«2040>  22300 

GO  TO  2019  22400 
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Table  8.  Continued 


9050  FORHkTC*  ••132I*.MI  22500 

90001  FORMAT (aOAl I 22000 

90002  format (•!' I 22TOO 

END  22B00 

C*”********************************************************************  22*«® 

SUBROUTINE  CASE  23000 

C— — — — . 23100 

C DETERMINE  tCASE  . 23200 

C— — — — 23300 

COMMON  /RA/X(10)i2llO)tTH(10)<H(10)tK(10)tR(10l«a(10)<M|10>«  23A00 

1 CURVCI10)tO]R(10)tXl(10l«X2<|0)«Zl(10UZ2llO).Tl(10>«T2llOlt  23500 

2 KSI(10ltETA(10ltRHOI10lfXO(lliiZO(llltTO<llt2)«KRV<3)«CL<2>*  23600 

3 KV<2)«VllO«3i3).SLI3i»RMOSL<10«3ltCRVSLaOt3l  23700 

A«VB(10f3t2ta)<CRVe  (10ia)tXMAPAXa0)«ZNAPAXfl0)«XSID<10lfZSlOllOI>  23B00 

5 BLBtai  ; . . 23900 

COMMON  /IA/ICASE(10>iJCASE(30l«NCASE<Tt«NP(3)iNS(3liN8<e>  2A000 

common' /RS/KttNI«KJtMJiMtJtLINE<BLNK»RINF*CRRMAXiOP0N«ELLXi  26100 

1 YUCLiZCRO»FORE«AFTEtTOPS«BOTS>OPUPtOOIV»P10iaOfXJiZU<XI»2t<Hl«  26200 

2 Rl.BlfTHl«aj<AliA2tA3tUlJtTHJtStTHIJf$IJ«AIfHJ,R3,A6«0X  26300 

StTYPCiASVNiUPFOiUPAFtUPTOfUPBOtONFOtONAFfONTOtONBO  26600 

COMMON  /IS/lRtlH.lFtlPtlPRiNTtNSEOfNSeOPltlSTOPtJSTOPiKSTOPi  26500 

1 IfJtLtlPlfEMltlltNPTStlZ^ISLiLBfN^NMlfia  26«00 

common  /LA/SLOPEaOI«10eN(80>>CARO<eOI  26700 

COMMON  /US/TltT  26000 

REAL  KtMiKSItKI>HIfKJfHJ«MlJ>LINE«KVtKRV  26900. 

Cft  LOOICAL*!  $LOPEitOENiCARO«Tl»r  25000 

DIMENSION  E(l0t7»  25100 

EQUIVALENCE  lElltUtX  (ll>  25200 

1 tlEIU2>«Z  (ll>  25300 

2 tlEI|p3)*TH(ll>  25600 

3 . tIEIlfAlpH  (ll>  25500 

6 .IE(lt5)*K  (lit  25000 

5 flEd.OltR  (in  25700 

6 tlE<l«7t»e  till  25800 

3000  continue  25900 

C DETERMINE  ND.  SEOMENTS  NSEO  20000 

00  3010  J61.10  20100 

I>10>J*1  20200 

ICA$EI1I>0  - 20300 

IFICURVEIlIpNE.BLNKIGO  TO  3010  20600 

NSEO-I'l  20500 

3010  CONTINUE  20600 

. MRlTE<IWt3015INSEG  20700 

3015  FORMAT t'ONUMBER  OF  SEOMENTS  NSEO  -•■13)  20800 

C make  sure  CURVEb'LINE*  OR  'OOtV  SPECIFIED  FOR  EACH  SCO  20900 

00  3060  I61,MSEG  27000 

DO  3020  3*1 f 2 27100 

IFICURVEIll.CO.KRVIJIlOO  TO  3060  27200 

3020  continue  27300 

HRITEIlM«303IIIIpCURVCII)tKRVIIItKRV(2l  27600 

3030  FORMAT!*  COND  CODE  30301  CURVEI»» 12*' I »A6i • NOT  ■ TO  OR  • 27500 

ltA6>  27600 

ISTOPal  27700. 

3060  CONTINUE  27800 

IFCISTOP.EQ.IISTOP  3030  27900 

DO  3080  lalpNSEO  28000 
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L«0 

DO  3050  JaItT 

IF  IE  1 1 1 .LT.R|NFII.aL*2aaj 
3050  CONTINUE 

IF (SLOPE  1 1 1 .eo.T I LaL*Z56 
IF(CURVEIII.Ea.KRV(2)ILaL*Sl2 
DO  3050  ja|,2l 
IF<L«NE.JC*SEIJII00  to  3060 
ic*sE<naj 
00  TO  3000 
3060  CONTINUE 
ISTOPal 

MRITEdMtSOTO)  I 

3070  F0RM6TI*  COND  CODE  30T0I  INPUT  OAT*  FOR  SEONENT  ••IZi*  DOES  NOT  N* 
ITCH  ANY  OF  PROGRAHMEO  OPTIONS') 

3080  CONTINUE 

MRITE(IW<3D90)  (I.ICASEIII  dal.ioi 
3090  FORMATIiOl/'  ICASE I 't IZt « I Bi«t3)> 

IFlISTOP.EQ.OlOO  TO  AOOO 
WRITE (lMt203S) 

STOP  3070 
4000  RETURN 

2035  FORMAT! 'OEXCCUTION  HALTED') 

END 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa.**.*. 

SUBROUTINE  SE6E 

C......— 

~C  DETERININE  SEGMENT  EQUATION 

C................ 

COMMON  /R*/X(10)»ZCIO)fTHI)0) tH< 10) «K (10) «R(10 I iB) 10) •HIIOI • 

1 CURVEI10)«01RI10)  tRiaO)fX2llO)«21llO)  (ZZIIO)  iTl  ( 10)  tTZIlO) « 

2 KSI(10)<ETA(10l>RHO(lOI<XOIlll.2aill).TOai>2).KRV|3),CL(2>> 

3 KVt2) •V(10«3«3) (SL 13) tHH0$Lll0i3) tCRVSL I10<3) 

AtVeil0i3>2f«) tCRVB  (lOtO) (XHAPAXIIO) iZMAPAX I 10) «X$ID I 10) tZSIO 110) • 
5 BLBIO) 

common  /1*/1C*SE<10).JC*SE<30)«NC*SCIT) tNPI 3) tNS (3) .NOIB) 

COMMON  /RS/KI.HI.KJ«MJ*HIJ.LlNE>BLNK«RINFiERRMAX.OPON.ELLX. 

1 VUCLtZERO*FORC«*FTC.TOPSiSOT$*OPUPt061V<P101BO.XJtZJ>Xl«ZItH|. 

2 RI>BI>THItBJ»*l«*2.*3>UIUiTMJtStTHIU<Sl«><*liHJ,RJ«6AiOX 
3<TVPE>ASVM<UPFO<UPAF.UPTO>UPeO«DNFOtON*FtONTOiONBO 

COMMON  VlS/IR.IW«lF«IPilPRINTiNSCG.NSEGPltlSTOPtJSTOPtKSTOP« 

1 l(J<LtIPl<LMI>ll«NPTS«12iI$L>LSiN«NMltI3 
COMMON  /LA/SLOPE ( 10) i IDENISOI fC*R0(80) 
common  /LS/TltT 

REAL  KtMiKSltKItNItKJiMJtMIJiLlNEtKVtKRV 
LOGICAL*!  SLOPE* IDEN.CAHOtTI.T 
4000  CONTINUE 
XJaO. 

ZJaO* 

DO  4990  laUNSEG 
IlalCASEII) 

IFIt.EO.DGQ  TO  4007 

IF(tC*5Elt>l).NE.5.AM0.IC*SEII-ll.NE.9)00  TO  4007 
00  4006  Jal,T 

IFIll.NE.NCASEIJDOO  TO  4006 
HRITC(Iw*400S)ICASEII)*1CAS£II>I) 


28100 

28200 

28300 

28400 

28500 

28600 

28T00 

28000 

28900 

29000 

29100 

29200 

29300 

29400 

29500 

29600 

29700 

29800 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31200 

31300 

31400 

31500 

31600 

31700 

31800 

31900 

32000 

32100 

32200 

32300 

32400 

32500 

32600 

32700 

32800 

32900 

33000 

33100 

33200 

33300 

33400 

33500 

33600 
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4005  format ( 7 .COND^COOE  4005J  CASE  ‘tlEi*  MAY  NOT  bC  OSEO  AFTER  CASE  •• 
1I2»  • . ; , ■ 

ISTOP-1  -r' 

00  TO  '4007  ■ “ • 

4000  CONTINUE  -7  . - 

4007  CONTINUE  -V 
XI«X<I) 

zi«z<n 

Ml4M(I>  ' -ly 

K1»KUI 

Ri4Ra> 

BI4B(I>  ' 

TIaSLOFEdI  - 
TMIaTMdl  vi 

C.  3 4 5 0 7 b 9 10 

C . -lIlTi'lZ  13  14  15  1 0 1 7 1 8 1 9 20  21 

00.  TI>.(4U7O>4090t4115«4l2U>4l5D>4OlO«4020.40JU<41SO«4lOO« 

--1  ' . 4170«41BOt4l90442UOt4300«4400«4410<4420t4430<4400(4490 
2)  til 
C ICASEaO 

4010  Mlt)«l21-2U>/(XI«XJ) 

THIllaATANIHM)  )/PIOiaO 
e(i)4Zi-Ma>4xt 
60  TO  4040 
C tCASe-T 

4020  H(l>«TAN|TH(t|4Pt01SOI 
- ZIIlaMU>a|Xltl«XJI*2J 
BII)4Z4I>-NII|4XI 
GO  TO  4040  % . . ' 

C ICASE-B 

4030  M4I)>TAN(TH<n4P|01Ba) 
xa)>xj*izi-zj>/M<ii 
B<l)<zt-H(l>*x<l) 

60  TO  4040 

C CHECK  REDUNUANT  INPUT  DATA 
4040  JSTOPaO 

IFIXl.LT.RINF.ANO.ABSIXt-XdH^GT.ERRMAXIjSTUHal 
IFIJSTOP.CQ.llMRlTEdb440Sa>  JSIOPd.CRRMAX 
IF(2t.Lr.RINF.AND.ABSlZl>ZdH.GT.ERRMAXI  J$T0Pa2 
IF  IJSTOR.E0.2IHRITEdM<40SO>JSTOP>I,ERRMAX 
IFiai.LT .RtNF.ANO^ABSIbl-bdl  ) aGT .ERRMAX)  JSTOPaS 
IFI3STOP.E0.3)WRITEdw«40S01  JSrOPf  ItERRMAX 
IFITMl4LT.RINF,AND«A8S(THi«THd) ) •GT.ERRMAXI3ST0P44 
IF <JST0P.EQ.4) WRITE (IM. 40501 JSTOPt I lERRMAX 
4050  FORMAT!'  COND  CODE  4050-'«lltd  REDUNDANT  INPUT  DATA  FOR  SEGMENT  • 
1«I2«'  NOT  self  CONSISTENT  TO  vlTHIN  'tlPEIS.TI 
IF IJSTOP.NE.O) WRITE (IH. 4055) 

4055  FORMAT !•  SUGGEST  USER  CHECK  IS  INPUT  DATA') 

GO  TO  4150 
C ICASEal  : 

4070  HUa-|XI-XJ)/<ZI-2JI 
B Ja2| *2 J>Mja  < X I *X J ) 

BJa0.54bJ 

Ala|XlaNJ4(2|>BJ)  >/d.*HJ442) 

A2a  I X I *Nja  ( 2 1'BJ)  I >■2'*  1 1 • *HJ4a2 ) • I X 1442*  (Z  I>b J)  aaj.R  1 *42 ) 

A2aS0R T ( A2) / ! 1 . »M J442) 


33700 

33800 

33900 

34000 

34100 

34200 

34300 

34400 

34500 

34600 

34700 

34800 

34900 

35000 

3S100 

35200 

35300 

35400 

35500 

35600 

35700 

35800 

35900 

36000 

36100 

36200 

36300 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

3/200 

37300 

37400 

37S00 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

38300 

38400 

38500 

38600 

38700 

38800 

38900 

39000 

39100 

39200 
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Hin«*l»A2  39300 

39400 

IF(0l9(li<Ee.0m>N.*N0.K<i>.LT.ZII60  TO  4000  39S00 

H(I|»«1-A2  39600 

K(IlsHJ4HI|l*ej  39TOO 

4000  H(I|a>IX(ll«H(l) l/(Z(I)>K(II>  39000 

TH(I>aATAN(HCII)/PI0100  39900 

00  TO  4130  40000 

C ICASEa2  40100 

4090  IFII.6T.II00  TO  4||0  40200 

lffllTE<l«f4100)  40300 

4100  ro«HATI>  COHO  CODE  *1001  FIOST  SEOMENV  MAY  NOT  BE  SPECIFICO  TANGCN  40400 
IT  TO  PREVlOVS  SeOHENT'l  40500 

WRITE IIH|2015I  40600 

STOP  4100  40700 

4110  MJ«-<Xt>XJ)/(ZI«ZJ>  40000 

e<«*zi*z«i-Hj*ixi«xji  40900 

BJ*0*S*BJ  41000 

M|«>l./TANIPI0180*TH(t>l))  41100 

BtaZJ~Rl*AJ  41200 

Htll  — <BJ«Bn/(HJ-NI>  41300 

K<ll-MJ>H(I)*aJ  41400 

Rn)aSORTCai>Hlin«*2*IZI-K(t)  >••21  41500 

THI1|mATAN||XI>H(III/<Z1-K(III  >/PI01S0  41600 

GO  TO  4130  41700 

C/i  C ICASEO  41S00 

4115  9>la  41900 

1F(DIRIII.NE>OPON>S»-1.  42000 

HI1|pXJ«R|4SIN(PI01S04THJI«S  42100 

Kill aZ J-RI4COS (PIOl BO^THJI •$  42200 

Z 1 1 ) BK 1 1 ) *S*SeRT IRI ••2- ( X 1-H  C I ) 1 ••Z)  42300 

TH(II«>ATAN< (Xl-H(III/(Z<I)oK(t> I >/P|OieO  42400 

60  TO  4130  42500 

C tCASEM  42600 

4120  NJ«-1,/TANIPI0180^TH{I>>  42700 

BJaZI-MJ^Xt  42800 

HIJ>-<X1-XJ>/<2I-ZJ>  42900 

6IJ«ZI*ZJ>H|J^IXI*XJ)  43000 

BIJaO«54BIJ  43100 

H(t)a>IBJ-BlJI/<MJ>NIJI  43200 

K(t)«MJ^H(l|*BJ  43300 

R«I)aSQRT  < UI-N< I) > ••2« (ZI»K( t ) >••21  43400 

GO  TO  4130  43500 

C CHECK  redundant  data  43600 

*130  KSTOPbO  43700 

IFCXI  .LT.RINF.ANO.ABSIXI  -X  (II I .OT.ERRH*X)KSTOP>I  43800 

lF(KSTOP.Ea.llHRITE(tWt4l40)KSTOPit<CRRHAX  43900 

IFIZI  .LT.RINF.ANO.ABSIZI  >2  (1) I .OT.CRRMAXIKSTOP-2  4*000 

lF(K5T0P.EO.2>WRITEIlHf4140>KSTOPiItERRN*X  4*100 

lF(TH|.LT.RlNF.AN0.Aa$(THl>TH(l)I.GT.CRRNAXIKST0Pa3  4*200 

tF(KSTOP<EO.3>HRlTE(|Hi4l40>KSTOP«ltERRMAX  4*300 

tF(HI  .LT.RtNF.4NO.AeS(Hl  •H  ( I ) ) .OT.ERRHXXI KSTOPa*  4*400 

IF<KSTOP.EO«*)WRITEIIMi*140)KSTOPtliERRM*X  44500 

1F(K1  .LT.RINF.ANO.ABSIM  •K  III ) .OT.ERRMAXIKSTOPaS  4*600 

IF(K$TOP.CQ.SlHRITE<lHf*140IKSTOPttteRRH*X  4*700 

IF  IRI  .LT.RINF.ANO.ABSIHI  aR  1 1 1 1 •GT.ERRNAXlKSTOPaO  4*800 
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IF(KST0P.Ee.6)HRITE(IM««U0)KST0P>t<ERRMAX  ««9a0 

4l«0  F0RMAT(«  cord  code  REDUNDANT  INPUT  DATA  FOR  SE6NENT  • 45000 

i>12«*  NOT  SELF  CONSISTENT  TO  WITHIN  •>1PC1S.T|  45100 

tF(K$TOP.NEiOIWRITC<IW*405S>  45200 

60  TO  4150  45300 

4150  CONTINUE  ’ - 45400 

XJaXII)  45500 

ZJ>EtI)'  45600 

THJ>TH(1»  45T00 

60  TO  4990  45600 

C ICASEalO  45900 

4160  Sal.  ' 46000 

IFIDIRIU.EO.OPUPISa-l.  46100 

H(l|aXI*SaRiaSIN(PI018>i4THI)  46200 

KII)aZIaS*RlaCOS(PIOiaoaTHI)  46300 

60  TO  4130  46400 

C ICASEal 1 . . 46500 

4170  sal,  46600 

tFIOIR(I>.e«.OPUPlSa.|.  46TaO 

Z<l|aK|*SaSaRT<Rtaa2>IXI>HI|aa2|  46600 

THIt)a-ATAN(<XI-Ht>/(2lt)*KII l/Pt016C  46900 

60  TO  4130  „ 47000 

C ICASEalZ  47100 

4160  5a|,  47200 

IFIOIRdl.EO.FORE.OR.OIRlD.EO.AFTElOO  TO  4166  47300 

L/i  WRITE (16,4105)  I 47400 

ON  4105  FORHATC  COND  CODE  4I05I  DIR  MUST  OE  EITHER  "FORE*'  OR  "AFTE"  F 47500 

. . lORE  SEOMENT'dS)  47600 

WfllTE<IW,2035)  47TOO 

STOP  4165  47000 

4166  CONTINUE-  . : ■ 47900 

IF(OIR(l).EO.FOREISa-l.  46000 

X(I)aHI«SaSeRT(Riaa2>(2I>K|)*a2)  46100 

TH(I)a.ATAN((X(II>HI>/)Zl«Kl))/PIOieO  40200 

60  70  4130  46300 

C ICASEaU  46400 

4190  S-l.  46500 

IF(OIRIt) .EO.OPUP>Sa-l<  46600 

X(I)aHl-SaRiaStN)Pt01S'faTHn  46700 

ZIIIaK|*SafliaCOSIPIO160aTHI)  46800 

60  TO  4130  40900 

C ICASEal4  49000 

4200  B(I|aZI-XIaTAN<PI01S0aTHll  49100 

60  TO  4040  ’ 49200 

C ICASEalS  49300 

*300  Z(I)aXlaTAN(PI0180*THII*dI  49400 

60  TO  4040  49500 

C ICASEal6  49600 

4400  X(t)a(ZI-BII/TAN(PI01BaaTHI)  49700 

60  TO  4040  49800 

C ICASEal?  49900 

4*10  IF(ABS(2II>>XII>aTAN(PIO180aTH(I))-B(I))>LT.ERRMAX>60  TO  4150  .50000 

USTI)Pa5  50100 

WRITE(IW,4050)USTOP,I,EHRH*X  S0200 

60  TO  4150  50300 

C ICASE-IB  50400 
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44Z0  IFUBSI<X<|>-NII)I**2*I2(I>>KII|>**2*R(I)**2>«LT.ERRH*X|60  TO  4150  50500 

KST0P«T  50400 

WRlTE(lHt4140>KSTOPtliErtRRAX  50TOO 

00  TO  4150  50800 

C ICASEbIO  50000 

4430  S«0.  51000 

IFIDIR<l>.Ea.UPFOI5-*l.  51100 

IFIDIR<l).Ca.UPAFISa*l.  51200 

IFIDIR<1I.EQ.UPT0IS»1»  51300 

lFI01R<I).Ea.UPBOISa»l.  51400 

IFIOIR<1>.EO.ONFOI5«*1«  51500 

IFIDlRIll.EO.DNAFISa*!.  51400 

IFIDIR(l>.Ea.ONTO)Sa-l>  51T00 

tFIOlRIll.EO.ONBOISa’l.  51800 

IFIS.NE.0.160  TO  4440  51000 

WRITE (IM«4431) 1 52000 

4431  FORHATC  COND  CODE  44311  DIRIl)  EQUAL  TO  ILLEOAL  VALUE  FOR  CASE  10  52100 

I FOR  SESHENT'tlS)  52200 

WRITE (1842035)  52300 

STOP  4431  52400 

4440  IFI|.N£*NSEGIC0  TO  4450  52500 

WRITE (IWt444l) t 52400 

4441  format (•  CONO  CODE  44411  ICASEalO  MAT  NOT  BE  USED  FOR  LAST  SEGNENT  52700 

l«tI3l  52800 

WRITE (IW42035I  52000 

STOP  4441  53000 

4450  H(l)aX3>SaRl  aSlNlPIOlWOaTMJI  53100 

KII)aZJ4SaRI  •COS<PI0180aTHJI  53200 

IF (DIR ( 1 ) .EO.UPFOt  DIR ( I ) aFORE  53300 

IF(DIR(I).EO.UPAFIDIR(I)aAFTE  53400 

lF(DIR(II.Ea.UPT0)01R(l)«T0P$  53500 

1F(DIR(1I  .EO.UPBOlDIRdlaOOTS  53400 

IF (DIR ( 1 1 .EO.DNFOIOIR ( 1 > aFORE  53700 

IF (DIR ( I ) .EG.ONAFIOIR ( 1 ) aAFTE  53800 

IF (DIR ( 1 ) .EG.DNTOIDIR ( 1 ) aTOPS  53000 

lF(DIR(l>4EQ.0NB0)DlR(I>aB0TS  54000 

ICA$E(I>a5  541 00 

60  TO  4130  54200 

C ICASEaZO  54300 

4440  INlal-l  54400 

AlaS0RT((Xl.H(lMl)|a»2.(ZI.K(IHl>>aa2t  54500 

A2aARCOS(R(|Hl>/Al>/PI0180  54400 

A3aATAN( (Z1<K(1M1>I/(X1-H(1M)>I l/PIOlOO  54700 

SaO.  54800 

1F(DIR(I>.EO.TOPSIS-»1.  54000 

IF(DlR(l>.Ea.BOTSISa-l.  SSOOO 

IF(S.NEiO«)ao  TO  44TD  55100 

WRITE(lWt4441>I  55200 

4461  FORNATC  COND  CODE  44611  DIR  MUST  BE  EITHER  TOPS  OR  BOTS  FOR  SE6ME  55300 

INT'iIS)  55400 

WRITE(iW.2035> , 55500 

STOP  4441  55400 

4470  IFd.NE.nOO  TO  4480  55700 

WR1TE(1H»44T1>  55800 

447 1 FORMAT (•  COND  CODE  44711  lCAS£a20  MAT  NOT  BE  USED  FOR  FIRST  SE6MEN  55900 

ITM  54000 
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URITE(IM«2035I 
STOP  4«7) 

««ao  THIII*A3-S«l90.-*2) 

M ( 1 > aTAN IP |01SO*TH 1 1 1 1 
B<naZI-MII|*XI 
00  TO  40«0 
C ICASEaZl 

4490  IFll.NE.llOO  TO  4500 
HR]TEIIW»449I) 

4491  FORMAT  I • COMO  COOE  44911  ICASE-21  MAT  NOT  UE  USED  WITH  FIRST  SESEM 
IENT*> 

HRITEIlNtZOaSI 
STOP  4491 

. 4S00  IF(CURVC<I<I)<EO«OOEVII>0  TO  4540 
C LINE-OGIVE 

IMlal«l 

HJ  ■TANIPIOlOOMHIIHll) 

AlBNJ4XII>*aiIMll 

$■•1. 

. 1F(A1.0T.ZI1)IS«-1. 

OJbBI INI t *S4R 1 1 lycos (P1U1004TH I IHl 1 1 
AI«1.*HJ442 

azb2.4 ( le J>z I I 1) •hjbx I 1 > I 

A3>X 1 1 1 ••2* iBJ>Z 1 1 1 1 ••2-R 1 1 1 442 

A4aA2/2./Al 

■ ' AS4A3/A1 

°°  AS4A4442.XS 

IFIAS.OE. 0.160  TO  4515 
IFIABStASl.GE.ERRMAXiGO  TO  4510 
IIRtTE<l»t450lll 

4501  FORMAT!*  COMO  COOE  45011  SEOMENT'.IS.*  LOCUS  OF  CIRCLE  CENTERS  00 
lES  NOT  INTERSECT*! 

WRITE<tR*2035l 
STOP  4501 
4510  A540. 

4515  ASaSORTIAS) 

A14.X44A5 

A24.A4-.AS 

A3aA14MJ*BJ 

A44A24MU*ej 

IFIAI.BT.AZIGO  to  4520 

HIM-Al 

KII)4A3  ■ ■ 

00  to  4530 
4520  H<I|4A2 
KlIlaAA 

4530  TNI1|4-ATANI1XI1).H(II1/IZ(I).K(IIII/Pt0180 
60  TO  4130 
C OGIVE’OOIVE 

4540  |Nl4|.l 

IFIOIRIIMll.Ea.OPON.OR.UlRIlHlI.EO.OPUPieo  TO  4550 

4541  HRITEIIH.4S42IIM1 

4542  FORHATI*  CONO  COOE  45421  SCCNENT*«I3f • OIR  NOT  EQUAL  TO  OPUP  OR  OP 
10N*I 

WR!TCIIWt2035l 
STOP  4542 


50100 
50200 
56300 
5O40O 
50500 
50000 
S6T00 
56000 
56900 
57000 
57100 
57200 
57300 
57400 
57500 
57000 
57700 
57000 
57900 
50000 
50100 
50200 
50300 
50400 
50500 
50600 
50700 
50800 
50900 
59000 
59100  ' 

59200 
59300 
59400 
59500 
59600 
59700 
59000 
59900 
60000 
60100 
60200 
60300 
60400 
60500 
60600 
60700 
60000 
60900 
61000 
6110O 
61200 
61300 
61400 
61500 
61600 
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lyi 

VO 


*SS0  ir<0IR<l) .EO.OPUP.OR.DIMII).EO.OPDN>60  to  *S60 

IH1»I 

60  TO  4541 

4560  S>*1> 

tFIO[RIIMlUEO.OIR<I)>Si»l. 

IF(ZIII.E0.K(IM1m00  to  4630 

8J«X 1 1 ) ••2»Z ( 1 1 •*Z»H CIHl ) ••Z-K < ! Ht » ••Z»R 41Ml|*(R|lHl) *Z.*S*R (II) 
ejaBJ*.S/(Z(l)-KllM|l| 

A2B2.4(iej-Z<II|4H,>-X(l)l 

*3ax [ I ) aaz*  <BJ>Z (III aa2-R  < 1 1 a*2 

*4a.5a*2/«| 

ASaxl/Al 

ASaA4aa2>AS 

IFIAS.GE.O.ieO  to  4500 

IF(AOS(AS),GE.ERRMAX)60  to  4570 

HR|TE<lPt4S61)I 

4561  FORMAT (•  COND  CODE  4561 « SEGMENT • t I 3 f • LOCUS  OF  CIRCLE  CENTERS  00 
lES  NOT  INTERSECT*) 

HRITE(|H«2035) 

STOP  4561 
4570  ASaO, 

4580  ASaSORTIAS) 

Ala.A4*AS 

A2a.A4-AS 

A3aN0aAl*BJ 

A4aMJaA2*BJ 

4585  IF(DIRIIHl>.EO.OtRII)ICO  TO  4622 

IF(D1R<IMI>.EO.OPDN.AND.DtRII).E0.OPUP)OO  TO  4590 
tF(0]R(lNll  .EOiOPUP.ANOiOIRID.EOiOPONIGO  TO  4620 
STOP  4581 

4590  IFIA4.6T»A3)60  TO  4619 
4600  MIDaAl 
K I I ) aA3 
GO  TO  4660 
4610  NIIIaAZ 
K<ll>A4 
GO  TO  4660 

4620  1F(A4«ST.A3IG0  TO  4600 
GO  TO  4610 

4630  A1  a|XII|aa2*Z(I)a*2-IHI]Ml|a*2*K|IMl)a«2|«R|lHl)a|R(|Hl)*2.a 
1 SaR(lll|a.S2(X(tl-HIIMl)> 

A2aR(l)aa2.(«l  -Xll)l*a2 
IFIAZ.GT.O.IGO  TO  4650 
IFIABSIAZI.GE.ERRMAXIGO  to  4640 
HR|TE(U«4561)I 
STOP  4562 
4640  A2a0, 

4650  A2a50RTIA2) 

A3aZII)*A2 

A4a2<|).A2 

A2aAl 

GO  TO  4585 

4622  IFIOIRIlMD.EO.OPDN.ANDtOIRID.EO.OPDNIGO  TO  4620 
1F(OIR(IMI).EO.OPUP.ANO.OIR(1I.EO.OPUPIGO  to  4590 


61700 

61800 

61900 

62000 

62100 

62200 

62300 

62400 

62500 

62600 

62700 

62800 

62900 

63900 

63100 

63200 

63300 

63400 

63500 

63600 

63700 

63800 

63900 

64000 

64100 

64200 

64300 

64400 

64500 

64600 

64700 

64800 

64900 

65900 

65100 

65200 

65300 

6S400 

65500 

65600 

65700 

65800 

65900 

66000 

66100 

66200 

66300 

66400 

66500 

66600 

66700 

66800 

66900 

67000 

67100 

67200 
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' M C • 

STOP  *623  67300 

*660  TH(II*-AT*N<<XtI)-H<in/IZIl)>K<Iin/PIOiaO  67400 

00  TO  *130  67S00 

*990  CONTINUE  67600 

IFIJSTOP.NC.0IWRITE<lW«2035)  67700 

IFIKSTOP.Ne.0IWRITEaWt2035l  67800 

IFI1STOP.NE.OIMRITE<1M<2035I  67900 

IFUSTOPtHE.OISTOP  *050  . 68000 

IFIKSTOP.NE.OISTOP  *1*0  68100 

IFIISTOP.NE.OISTOP  *005  68200 

RETURN  68300 

2035  F0RH8T(*0EXECUTI0N  HALTED')  68*00 

END  68500 

68600 

SUBROUTINE  INTR  68700 

C—————— 68800 

C CONPUTE  segment  INTERSECTIONS  FOR  IC«SE>5.9  68900 

C»— 69000 

COMMON  /RA/X(lO>tZ<IO)«TH<lO)iH<l0)tK(lO)fRllO).B<lO)>M(lO>>  69100 

1 CURVE<lO)fO|R(I0)fXl(10)»X2(l0>«21C10)«Z2(10>»Tlll0>>T2<10)«  69200 

2 KSII10l.ET*(10).RHO(10l>XO(ll>.ZU(ll).TQ<Il«2>,KRVl3I.CL<2>t  69300 

3 KV<2l>V(10t3t3lfSLI3)»R<0SL(|0i3>tCRVSL<l0f3)  69*00 

*»VB<IOt3.2»8)«CRVB  1 10«8)  iXMAPAX  1 10)  tZMAPAX  1 10)  ,XS10(  10)  (ZSIOdO)  t 69500 
5 BLBIB)  69600 

COMMON  /l*/lCASEI10)>JC*$El30)<NCASE(7)«NPI3)tNS(3l*NBIB)  69700 

Ov  COMMON  /RS/KI.MI«KJ<MJ.MIJtLlNC,BLNKtRINF<ERRM*XtOPON»ELLXt  69800 

O 1 TUCLfZEROtFOREf  AFTCtTOPS«eoTS«OPUP.06IVtP10180tXJ>ZU>XliZIfH).  69900 

2 RI<BItTHI<BJ<Att*2fA3«UIJtTHJtSfTHIU<SlJ»AI«MUtRJ«A*fOX  .70000 

3«TYPE**SYM,UPFO«UPAFiUPTOiUPBO<ONFOtONAF«ONTOiONSO  70100 

COMMON  /IS/IR>lH<IFt|PtlPRtNT.NSEG»NSEGPltlSTOP»JSTOPiKSTOPt  70200 

1 I«J«L«IP1«LMI«JI.NPT5>12>I5L«L8«N«NN1.I3  70300 

COMMON  /L*/SLOPEI10)ilUEN(a0)«C*RDI80)  70*00 

COMMON  /LS/TI.T  70500 

REAL  KtMtKSltKItMI.KJiMUiMlJtLlNEtKViKRV  70600 

LOGICAL*!  slope* IOEN«CAHO«TltT  70700 

5000  CONTINUE  . TOBOO 

DO  5320  I*ltNSEO  70900 

IPUI'l  71000 

IFIICASEII) .NE.5.AN0.ICASEII) •NE*9)G0  TO  5320  TIlOO 

. 1F<I.NE.NSEG)G0  TO  5020  71200 

NR1TCIIW»5010)  71300 

SOlO  FORMATC  COMO  CODE  50101  EQUATION  IS  INSUFFICIENT  FOR  LAST  SEGMENT  71*00 
IM  • , ■ - 71500 

MRITE(1H«2035)  71600 

STOP  5010  71700 

5020  CONTINUE  71800 

IFIICASEII) ,EO.S.ANO.CURYEItPl).EO.OGIV)60  TO  50*0  719Q0 

lF<IC*SE(I),Ee.5.*ND.CURVE<IPl).Ea.LINE)60  TO  5250  72000 

IF<|CASE(I),Ee.9.AND.CURVE<lPl).E(I.OGlV)GO  TO  5260  72100 

IF(1C*S£I1).E0.9.AND>CURVEUP1).EO.LINE)GO  to  52*0  72200 

WRITEdWtSOiS)  72300 

S025  FORMATC  COMD  CODE  50301  OOPSI  PROGRAMMER  ERROR.  •)  72*00 

NRITE<IW«2Q35)  72500 

5030  STOP  5030  - 72600 

C OGIVE'OGlVE  INTERSECTION  72700 

50*0  CONTINUE  72800 
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THIJaATAN2IKI(Pl>>K<I)(H(IPl)-HCl) I/PIOISO  72900 

siJ>saRT(iH(iPt>-H<ni*<2*iK<iPii-K<in>*2  > raooo 

Ala  (R(n*»2*SU**2->RIIPl>««2l/l2.*RCt)*SIJI  73100 

tF(ABSUI)  .9T.t.)AI*SI6Nll..AI>  73200 

AlaARCOSCAll/PIOlOO  73300 

X|aH(II*R(ll*COSI ITHIJ*AI1*PI0180>  73A00 

X<I|aRll>*RII)*COS<(7HlJ-An*P|OIOO)  73S00 

tFIAeS(X(II>XII  .LE.ERRHAXIGO  70  5090  73000 

lFII.LT.NSE(}.ANO.SLOPE(l«ll.Ea.7>GO  70  5070  73700 

1F<DIR<II.CO.FOREIOO  70  5060  73800 

IFIOIRilUEQ.AFTElOO  70  5080  73900 

WR17E(1M«50»0)1  7A000 

5050  F0RNA7(I  COND  CODE  50501  DIR  MUST  BE  EI7HER  "FOREX  OR  ••AFTC»i  F 70100 
lOR  SEGMEN7  •tt2>  70200 

WRITE (IHf2035>  70300 

STOP  5050  70000 

5060  IFIXl.OE.XdllOO  TO  5070  70500 

5065  X(I>>XI  70600 

Z<[)«KIII*RI1>*SIN(PI0180*{THIJ*AI>>  70700 

GO  TO  5230  70800 

5070  Z(I}aK(I)*RII>*5INIPI0180«|THIJ-Al>>  70900 

GO  TO  5230  TSOOO 

5080  IF(XI.GT.XII>>00  TO  5065  7SI00 

GO  TO  5070  75200 

5090  ZI«K<I>*R<I)OSlNCPI0100»(THia*Atn  75300 

Z(I»-KIII»RI|>*$INIPI0100*(TH1J«A1I|  75000 

IF<ABS(2<I>>ZI>.LE.ERRMAX>00  TO  5230  75500 

IFCOIRdl.EO.TOPSIGO  TO  5200  75600 

IF(DIRCII.E0.B0TS)OO  TO  5220  75700 

WRITEdWfSlOOl  I 75800 

5100  FORMATC  CONO  code  SlOOl  OIR  HUST  BE  EITHER  "TOPS"  DR  "BOTSit  F 75900 
lOR  segment  'tlZI  76000 

WRITE  dWtZOlS)  76100 

STOP  5100  76200 

5200  IFlZI.GT.Zdl  100  TO  5210  76300 

5205  Xd|aHdt*Rdl*CO$<PIOiaO<ITHlJ>AIH  76400 

GO  TO '5230  76500 

5210  ZdOZII'-  76600 

Xd|aHdl*R(I|*C0S<PI0180«<THU*Al»)  76700 

GO  TO  5230  76800 

5220  IFlZI.GT.Zdl  100  TO  5205  76900 

GO  TO  5210  77000 

5230  THdl»ATANI(Xdl-H<ll)/IZd)-Kdlt)/PI0100  TTIOO 

GO  TO  5320  77200 

C LINE-LINE  INTERSECTION  77300 

5240  Md  laTANIPIOIOOaTHd  ))  77400 

MdPtlaTANIPIOleO»THdPl>)  T7500 

Xdl  — (edPll-Bd)>/IHdPll-N<tl>  77600 

Zd)><Hdl4edPl>-M<lPl|4Bdl  l/IHdl-HdPl)>  77700 

GO  TO  5320  77800 

C OGIVE-LINE  INTERSECTION  T7900 

5250  CONTINUE  78000 

H dPl  I ■TAN<PIO1804TH  ( IPl  1 1 78100 

HlaMdPlI  78200 

BlaBdPlI  78300 

HJaHd)  78400 
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K JaK  < 1 > T6S00 

RJaR<Il  ' ' 76600 

00  TO  S2T0  70700 

C LINE'OOIVE  INTERSECTION  70000 

S260  CONTINUE  70000 

MII|aT«N(RI0100*TH(I|l  79000 

MlaHIII  79100 

Blaoai  79200 

HJaHlIPlI  79300 

KJaK(IPl)  79400 

RJaRlIPlI  79500 

00  TO  5270  79600 

5270  Alal.4Hla*2  79700 

A2aNI*IBl’>KJ)-HU  79B00 

A3aHja42*IBlaKJ)*a2-RJ*42  79900 

A4aA2**2-Al*A3  BOOOO 

IFCA4.0T.O.IOO  TO  5273  80100 

IFIABS<A4).CT.eRRMAX»0O  TO  5272  60200 

HRITEItH»52TltI  80300 

5271  FORNATO  CONO  CODE  52711  SORTKOI  FOR  SE0NENT<i|3«'  NO  INTERSECTIO  80400 

IN.'l  80500 

NRITE1lllt203St  ■ 80600 

STOP  5271  80700 

5272  A4a0.  80800 

5273  A4a$0RTCA4l  80900 

' XIa(-A2*A4|/Al  . 81000 

X<l)a(-A2>AAI/Al  81100 

IFIAB5<X(II>XI).LE.ERRHAXIC0  to  5300  81200 

IFlI.CT.NSEO.ANO.SLOPEdall.EO.TlOO  TO  5300  81300 

lFIDtR<I).e0.F0RE>00  TO  5290  81400 

IFIOIRIII.EO.AFTEIOO  TO  5310  81500 

WRITE* tHt52aO It  81600 

5200  format**  CONB  CODE  52801  DIR  MUST  BE  EITHER  ••FORE"  OR  ••AFTE**  F 01700 

lOR  SE6HENT*«t2l  OIBOO 

WRITE*tM«2015)  61900 

STOP  5200  02000 

5290  IFIXI.ST.XdllOO  TO. 9300  82100 

5295  XlllaXt  82200 

5300  Z(llaMI4X*|)*BI  82300 

00  TO  5320  02400 

5310  IFIX1.GT.X*1I|00  TO  5299  62500 

GO  TO  5300  82600 

5320  CONTINUE  82700 

RETURN  02000 

203S  FORMAT<*OEXECUT10N  HALTED*)  82900 

END  83000 

C*aaaaaaaaaaa>*aa**»a*a»a*«aaa*aaaaa*aaaa*aa*aaaaaa«*aa«**aa*aaM****aa»  03100 
SUBROUTINE  QART  83200 

........ .....  03300 

C SET  UP  FINAL  ARRAYS  FOR  OUICK  MODEL  83400 

C.. ..........................  03900 

COMMON  /RA/X(10)>Z<10>>TH<10>»H(10ltK(t0l>R*10>iB(10)4M*10>.  63600 

1 CURVE*10)  »0IRil0l<X|*10l  tX2*)0)  >21  (10) .22*  101  tTl  1101  (TZOO).  63700 

2 KSI(IO)  tETA(10).RHO*lOI  .XOIl  1 1 .20*11) .TO* 1 1.21  iKRVOI .CL*2I  t 83000 

3 KV(2).V(10t3.3).SLI3).R  0SL*10.3).CRVSL*10.3)  83900 

4.VB*10.3.2.8I.CRVB  I 10.01 .XHAPAX(IO) .ZHAPAX 1101 .XSIO(IO) .ZS10( 10) • 84000 


AEDC-TR-77-89 


Table  8.  Continued 


On 


S BLBCS)  Ml  00 

COMMON  /I*/tC«se<IOltJC«»E(3SltNC*SE(T)tNPO)«NS<3t>N9<8>  M200 

COMMON  /RS/KI«MI<K3tMJ<HIJiLlNEtBLNKfRINF>ERRHAX«OPDN(CLLK(  M300 

1  VUCLtZCROtFORE>RFTEtTOPS«eorS<OPOP«001ViPIOieOtXJ«ZJ>XI.ZI,HIt  a«400 

Z RIiBItTHItBJ«*l<*2>A3>eli>(TMJ>SiTHt3«SIJt*ltH3.RJ«A*>OX  MSOO 

3 1 T YPE I ASYMiUPFOiUPAF  • UPTOfUPBO  t ONFO  tONAF  t ONTO  • DNSO  M600 

COMMON  /I$/IR«IMtIF<IPtlPRINrtNSE6iNSE6PI>lST0PtJST0P«KST0Pi  MTOQ 

1 ltJiLt|Pltl.Ml,Il,NPTS»I2«ISL>LBtN>NN|.I3  OASOO 

COMMON  /LA/SLOPE(10) tIOEN(80l «CAROI8t>  8A900 

COMMON  /LS/TIiT  85000 

real  KfMiKSItKIiMtiKJiMJiNlJ>LINC«KV.KRV  85100 

LOOICAL*!  SLOPCtIOEN«CAHDtTttT  85200 

C TOlItll  IS  AT  BEGINNING  OF  SEGMENT  1 85300 

C T0<|«2I  IS  AT  END  SEGMENT  I 85A00 

C XO<I)  AND  ZO(I>  ARC  AT  BEGINNING  OF  SEGMENT  I 85S00 

6000  CONTINUE  85600 

XOIIIMO.  85T00 

ZO(1>bO>  85800 

00  6010  iBltNSEG  85900 

IF (CURVE ( I I .EG .KRV (2) ) CURVE ( I I bELLX  86000 

X0(I*1>«X(II  86100 

ZO(I*nBZ(II  86200 

6010  CONTINUE  86300 

DO  6030  IbIiNSEG  86400 

IP1b1*I  86500 

IF(CURVE(1I,NE.LINEIG0  TO  6020  86600 

to(I«iibatani izo(iPi>-zeili )/(xo(iPi>-xa(i>i )/Pioiao  86T00 

T0II»2|bT0(1»I>  86800 

GO  TO  6030  86900 

6020  IF(ZeiII.E0.K(l)>00  TO  6025  87000 

TO(ltllB-ATAN(  (XQIII>H(l>)/IZO(lfK(|H)/PI01BO  87100 

GO  TO  6024  87200 

6025  T0(ttll490.  87300 

6024  TO(l«2|a>ATANI(XQIIPl)>H(I||/IZG(IPll-K(|)|l/PI01BO  87400 

6030  CONTINUE  87500 

NSC8PI«NSEG«I  87600 

C SET  ALL  RIGHT  ANGLES  (90  OEGt  TO  ''IBM  RIGHT  ANGLES"  (89.999  OEG>  87700 

DO  6040  IbIiNSEGPI  87800 

DO  6040  JbI.Z  87900 

iF(TO(ltJI.6T.R|NF)GO  TO  6040  88000 

1F(  ABS(TO( If Jl I .LE. 89.999)60  TO  6040  88100 

TOdf  J)BSIGN(89,999tT0<lf  Jl ) 88200 

6040  CONTINUE  88300 

RETURN  88400 

END  88500 

C44444444f4444*4*«««44***4»444444»f«4«4f4*«*««*««»***f»»<**f*»»«*»««.,«,  88600 

SUBROUTINE  PRNT  88700 

C-..— 88800 
C PRINT  RESULTS  88900 

C..............  89000 

COMMON  /RA/X(10).Z(10><TH(10)fHllU)fK(10)«RllO)fB(10)tH(10)f  89100 

1 CURVE(10)fDIR(10)fXl(10)fX2(10).ZII10)fZ2()OI f T1 (10) •T2(10l . 89200 

2 KSKIOI  tETAdOlfRHOdO)  *X0 1 1 1 1 t ZQ(  1 1)  • TO  1 11  f2>  tXRVO)  tCL  (2)  f 89300 

3 KVI2>fV(10f3f3>fSLO>tRHOSL(10f3).CRVSLIIOf3)  89400 

4fVB[  10. 3.2.81  .CRVa  (10.8)  .XMAPAXdO)  .ZMAPAX  ( 10)  .XSIDdO)  .ZSIOdOl  . 89500 

5 BLB(8I  89600 
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COMMON  /IA/ICASE(10><JC*SCI30liNC*SCITI.NP<3>tNS(3t.NB(ai  80700 

COMMON  /RS/KI*MI«KJtNJ»MlJiLINEtBLNK<RINF»ERRHAA«OPONiEULXt  80800 

1 YUCLfZERO.FOREtAFTEiTOPSiBOTSiOPUPiOOIVtPIOUOtXJfZJtXIfZltHlt  80900 

2 RIt81<THl.B3iAliA2tA3t8IJ«THJ<SiTHIvl>SlJtAliHJ,RJ«A4«DX  90000 

3«TYPEtASYHtUPFOiUPAFiUPTOiUPBOtDNFO«DNAF«ONTOfONBO  90100 

' COMMON  /IS/IPa]U.lF>IP<lPR[NT.N9EG.NSEOPltIST0P<JST0P>KST0P«'  90200 

1 ltJrL»IPl>LH|<IliNPT$tt2ilSL>LBaNfNMltI3  90300 

COMMON  /LA/SLOPEMO)  iIOEN(80>«CAROI80)  90400 

COMMON  /LS/TIaT'  90500 

REAL  KfH.KSIiMlaMItKJtMJ.MIJ.LlNEaKVaKRV  90600 

LOOICAL*!  SLOPEalOEN.CAKOaTlaT  90700 

TOCO  CONTINUE  90800 

HRITE<lHt70a5>  90900 

7005  format COSEGMENT  END  POINT  COORDINATES  AND  SLOPES! • I 91000 

MRITE(IWtT010>  . 91100  ' 

TOlO  F0RMAT(*0*i9  Xt>X*al3Xt'2!<12X«tTA>«12X>«TF*/<  91200' 

00  7020  t>lall 91300 

- HRITE(IWt9000a)I<XO<ntZO<I>aTQ(lt2)«TO(ltl>  91400 

7020  CONTINUE  - 91500 

WRITE  (ItltT025>  91600 

7025  FORMAT  I 'OSEeMENT  EOUATlONSIt)  91700 

NRITC<IWt7030>  91800 

7030  FORMAT!  'O't  9X» 'X*?  13Xt 'H*  aUXt 'K' tlSXf 'RlaiaXi  tHt  ,|3X»  >7  • « 91900 

2 13X.>B<./'  •«100l»-*)l  ' 92000 

00  7040  IMI.IO  ' ■ • 92100 

ON  WRITEIIW>900aBIIiXlIlaH(tl«KII>aR(I»«H(tl>THIl),B(I)  . 92200 

7040  CONTINUE  - 1'  92300  ' 

RETURN  - 92400 

90008  FORMAT!* '••I2tlP9E14,6>  92500 

END  92600 

C***********************************************************************  92700 

SUBROUTINE  ROFL  92800 

— 92900 

C INSERT  ROUNDS  AND  FILLETS  93000 

C— — — — 93100 

COMMON  /RA/XllO><Z!10t aTHIlOlaHIlOt •KIIOI •RI10|.B<10).M«10>.  93200 

1 CURVEIlOl.OIRilOl tXI (10).X2ll0|fZl llO>aZ2!10l<Tl!IO).T2!10l t 93300 

2 KSI!lO)«ETAI10I.RHO(10)«XO|ll),ZOIll)«TOIIl«2)tKRVi3>iCLI2l«  93400 

3 KVI2>tV!10<3<3l<SLI3l«RHOSLI10t3)tCRVSL!l0t3)  93500 

4«VBI10«3a2a8l aCRVB  !10t8) f XMAPAX !10| aZMAPAX I 10) aXSIDllO) tZSID! 10) t 93600 
5 BLB(8>  93700 

COMMON  /lA/ICASEIlO) IJCASEI30) •NCASE!7> tNP!3) iNS!3) iNBIB)  93800 

COMMON  /RS/KI<MI«KU«MJtMtJlLINE»BLNKaR|NFfERRHAXiOPONiELLXt.  . ‘.93900 

l-YUCLiZEROaFORCaAFTEtTOPSteoTStOPUPiOOlVaPIOlOOaXJaZJiXIaZlaHIf  94000 

2 R|<BlaTHI«BJ<Al«A2.A3<BIJ.THJtS«THIJ,StJ>AlaHJ.RJtA4.0X  94100 

3<TVPEiASYH*UPFO>UPAF<UPTOaUPBOaONFO.ONAF«ONTO«ONBO  94200 

COMMON  /IS/IRaIW«lFalP<IPRlNTaNSE&*NSEGPlt|5TOPtJSTOPaKSTOPi  94300 

1 I<J>L«tPlfLMlillaNPTS>I2tISL>LBtNaNN)tI3  94400 

COMMON  /LA/SL0Pe!l0>«I0EN(80)aCARD!a0)  94500 

COMMON  /LS/TIaT  94600 

REAL  NaN.KSItKI.MIiKJiMJ.HIJ.LlNEaKV.KRV  04700 

logical*!  SLOPE«tOEN,CAHD<TIfT  94800 

6100  CONTINUE  94900 

--  - ISKIP*0  ..  95000 

DO  6800  iKltO  95100-. 

IPlalal  95200 
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tFIRH0(II.0T.«>INF>60  TO  6000  95300 

ISKtF>t  95*00 

tFCCURVEdI  .EO.LINE.ANO.CURVE(IPl)  •CQ.LINEtOO  TO  6110  95500 

IFICURVEII)  .EQ.LINC.ANO.CURVEdPll  .EO.ELLX>SO  TO  6200  95600 

IF(CURVE<t>.EQ.ELLX.AN0.CURvEdPl>.Ee.UINE)O0  TO  6300  9ST00 

IFCCURVECn.EO.ELLX.ANO.CURVEdPlt.EO.ELLXtOO  TO  6600  95600 

MtTEdM«6ieil  I 95900 

6101  FORMATd  CONO  COOE  6101S  PROGRAM  OR  USER  ERROR*  SEGMENT<tI3>  96000 

WRITE  dN*203SI  96100 

STOP  6101  96200 

C LINE'LINE  FlLLETTING/ROUNOINO  96300 

6110  CONTINUE  96*00 

lFIABS(T0d.2)-TQdP|il>).GT.ERRHAX>G0  TO  6120  96500 

WRITEdW.6111)I*|Pl  96600 

6111  FORMAT  I • CONO  COOE  61111  FILLETS/ROUNOS  HAT  NOT  BE  INSERTED  BETMEE  96TO0 

INLINE  SEGMENTS  WITH  SAME  SLOPEiSEGMENTS* *216)  96600 

WRITE  dH. 2035)  96900 

STOP  6111  9TO00 

6120  Sal,  97100 

IFITOdPltll  .6T.Tad,2l  >Sa.|.  97200 

B I aB  1 1 1 aSaRHO  d I /COS  (PI  0 1 60*70  (I  • 2 M 97300 

BJaB(|Pll-S*RHOdl/COS(PI0160*TadPltl>»  97*00 

Md>aTAN(P10180aTMd))  97500 

MdP|)aTAN(P10ie0*THdPlM  97600 

ETAd)B(HdPl|aBl-Md>*aj)/(HdPl)-Mdl>  97700 

OS  KSIdla*(BJ>BI)/INdPll-Nd)>  97800 

^ Ald)aKSId)-S*RHOd>*SlN«PI0160aTHd  ))  97900 

X2  d I >KSt  ( I ) -SaRHO  ( t > •SIN(P10160*THdPl>  1 96000 

21  d)*ETAd)*SaRHOd>*C05(P10180*THd  >)  96100 

22 (I ) bET A d > aSaRHO ( I ) aCOS (PI0I80*TH dPl } 1 96200 

TldiaTHd)  96300 

T2d)BTHdP|)  96*00 

GO  TO  6600  96500 

C LlNEaOGIVE  FILLETTING/ROUNDING  96600 

6200  continue  96700 

C OGlVE-LlNE  FILLETTINO/ROUNDINO  98600 

6300  CONTINUE  96900 

6201  SIbO*  99000 

S2b0.  99100 

53*0.  99200 

tl"!  99300 

IF (CURVE ( IPl > .ER.ELLXI llalPI  99*00 

I2«l  99500 

lF(U.C0.ni2aIPl  99600 

IF(T0(Ii2)-TadPl. 1116210*6205*6230  99700 

6205  WRITE ( lw*6206) 1* IPl  99800 

6206  FORMAT (>  CONO  CODE  62061  F IllETS/ROUNOS  MAY  NOT  BE  INSERTED  BETWEE  99900 

IN  segments  with  same  slope*  SEGMENTS* *213)  100000 

WRITCdW*2035)  lOOlOO 

STOP  6206  100200 

C TH(IPl>>TNd)  100306 

6210  IF(Kdl>-20(IPll  16220*6215,6225  100*00 

6215  WRITEdW*6216)I«IPl  100500 

6216  FORMAT (•  CONO  CODE  62161  K a 2 NOT  A PROGRAMMED  OPTION*  SEGMENTS**  100600 

I 213)  100700 

WRITE  dw*2035l  lOOBOO 
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STOP  6216  100900 

6220  SlBl.  101000 

S2al.  101100 

S3-1*  101200 

00  TO  6250  101300 

6225  Sl«l.  101400 

S2«-l.  101500 

$3al.  101600 

60  TO  6250  lOlTOO 

6230  1FIK<I1I-ZO<IP1)I6240<6235|6245  101000 

C THCIPIKTHIII  101900 

6240  SlB-1.  102000 

S2B-1.  102100 

S3B-li  102200 

GO  TO  6250  102300 

6235  HRITE(IWt6236)I«IPl  102400 

6236  FORHATC  CONO  CODE  62361  K ■ Z NOT  A PROGRANHEO  OPTION^  SEOHENTS'f  102500 

1 2131  102600 

HRlTEClHiZOSSI  1O2TO0 

STOP  6236  102000 

6245  SlB-1.  102900 

S2a}.  103000 

S3a-1.  103100 

GO  TO  6250  103200 

6250  CONTINUE  103300 

C^  BI«BI12)*S1*RHO(II/COS(PI01B04THI12>)  103400 

A1b1.*M(121**2  103500 

A2b2.*(IBI-K(I1)I*M(12)-H<I1))  103600 

A3aH  (111  MZ*  IBI-K  ( 1 1 1 1 ••2-  (R 1 1 1 1 *S2>RH0  ( 1 M ••2  1O3T00 

A4bO.S*A2/A1  103800 

A5BA3/A1  103900 

ASbA4**2-A5  104000 

IF(AS.OE.O.|GO  to  6260  104100 

IF(ABS(A5I.GT.ERRHAAIG0  TO  6254  104200 

A5a0>  104300 

GO  TO  6260  104400 

6254  WRITC(lMt62S5IItlPl  104500 

6255  FORMAT  I • CONO  CODE  62551  REAL  SOLUTION  TO  FILLET/ROUND  PROBLEM  00  104600 

lES  NOT  EXIST,  possible  PROGRAM  ERROR.  SEGMENTS  '.2121  104T00 

WRITE (IMt2035l  I04BOO 

STOP  6255  104900 

6260  ASbSORTIAOI  105000 

AlB«A4»A5  105100 

A2aN(t2l*Al*BI  105200 

A6bHII1I*R(11I»|A1-H(I1)I/S0RTIIA1-N(I1|I««2*|A2-KII1|)*«2I  105300 

ATb|M(I2I*(A2«BII2II*A1I/I1.*N(I2)*«2I  105400 

IF(X<1I.0E.AMIN1(A6.ATI.AN0.X(1|.LE.ANAX1IA6.AT|>G0  to  6265  105500 

Ala>A4*A5  105600 

A2bM(12|4A1*BI  lOSTOO 

A6aH(ll)*Rlll)*(Al>H(Ill  I/SORTKAI-HIIII  I••^*(AZ-K(I1>)••^)  lOSBOO 

ATa|M(I2)*(A2-B(12ll*All/ll.*N(I2l*«2l  105900 

1F(X(II.GE.ANIN1(A6.a7).AN0.XI1|.LE.ANAX1(A6.AT|I00  TO  6265  106000 

WRITE(IW.6261II  106100 

6261  FORMAT  I • CONO  CODE  62611  PROGRAM  UNABLE  TO  CHOOSE  BETWEEN  TWO  FILL  106200 

ICT/ROUND  SOLUTIONS.  SEGMENT* .131  106300 

WRITE! IW.20351  106400 
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STOP  6261 
6265  KSIII)-*1 
ET«I1)>A2 

IF(CURVE<I>.E0.ELLXI60  TO  6270 

X1(I)-AT 

X2I|»bA6 

Z2n>aK(ll)*R(Ill*IETAIll-K(IlM/SaRT(Ul-H(Il)|«*2*IA2-Klllll**2) 

Zl(llBM(I2)*Xl(l>*e(I2l 

TllItBTHIl)  . 

T2inB-ATANI(X2(I)-HIIl>>/<Z2<I)-Kllll))/PI0ia0 
GO  TO  6800 
6270  XIIIIBA6 
X2II>bA7 

Z2(IlaM<t2)*X2(II*B(I2) 

Zlll|BK||l»*RIIl)*CETAin-K(Il))/S0RTI(Al-HIIlM*>2*U2>KIIIt)»*2) 

Tl(I|a>ATANI(XIII>-H(IlM/<Zl(I)-Kinii)/P10160 

T2(IIbTH(IP1I 
GO  TO  6800 

6600  CONTINUE 

C OGIVE-OGIVE  FILLETTING  ROUNDING 
AlBTOIIPlfll-TO(I>2l 
A2«ZiI>-K(I> 

A3bZ(I>-K(IP1) 

IFIA1.GT.ZERO.ANO<A2>GT>ZERO.AND.A3,GT«ZERO)00  to  6610 
lFIAl.LTiZERO.ANO.A2,OT.ZEROtANO«A3.0T.ZCRO>00  TO  6615 
1FU1*LT>ZERO«ANO.A2.LT.ZERO.ANO>A3.LT.ZEROIGO  to  6610 
IFIA1.0T«ZER0.AND.A2.LT.ZER0.AN0.A3.LT.ZER0IG0  to  6615 
IFIAULT. ZERO. AND. A2.6T«ZERO.ANO«A3.LT>ZERO)00  TO  6620 
IFIAl.OT. ZERO. AN0.A2.BT. ZERO. AND. A3. LT.ZEROIOO  TO  6625 
IFIAI. 6T.ZER0.AND.A2.lt. ZERO. AND. A3. GT.ZEROIGO  TO  6620 
IFIA1.lt. ZER0.AND.A2.lt. ZERO. AND. A3.GT.ZERO>60  TO  6625 
WRITEIlW.6601II»IPl.AltA2.A3 

6601  FORMATC  COND  CODE  66011  SEGMENTS' t2I3t • t Al(A2.A3  HAT  NOT  BE  ZERO 

1.  A1«A2iA3b'«1P3E1S.7> 

HRlTEIIHtZ035t 
STOP  6601 
6610  SlBl. 

S2b1. 

GO  TO  6630 
6615  SlB-1. 

S2b-1. 

GO  TO  6630 
6620  Sl—1. 

S2b1. 

GO  TO  6630 
6625  Sl-1. 

S2a-1. 

GO  TO  6630 
6630  CONTINUE 

RIaR(ll*SlaRHO(II 

RJaR(lPl).S2aRH0II> 

IFIKlII.EO.KdPlDGO  TO  6665 
NIa.<H|IPl).HII))/IK(IPll-K(I)) 

Bla  |H(  IPl  I aaZ-HI  I ) aa2*K  ( IPl  I aa2.K  ( 1 1 aB2.Rj»2.Rj».2| /2  / ,K  ( ipi 
1 KID) 

Alal,.Hlaa2 


106500 

106600 

106700 

106800 

106900 

107000 

107100 

107200 

107300 

107600 

107500 

107600 

107700 

107800 

107900 

108000 

108100 

108200 

108300 

108600 

108500 

100600 

108700 

108800 

108900 

109000 

109100 

109200 

109300 

109600 

109500 

109600 

109700 

109800 

109900 

110000 

110100 

110200 

110300 

110600 

110600 

110600 

110700 

110800 

110900 

111000 

111100 

111200 

111300 

111600 

lllSOO 

111600 

111700 

111800 

111900 

112000 
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A2a2.»IIBI-K<I)l*Ml-Hll>)  112100 

«3>H( I ) IBI'K ( I I ) ••Z-RlaaZ  1 12Z00 

A*a.S*AZ/*l  112300 

A5aA3/Al  112A00 

A5aAA»Z-AS  112500 

IP IAS.6E. 0.160  TO  6635  112600 

IFIABSIA5).LE.ERRHAXI00  TO  6636  112700 

HRITEIIW(6631)I«IP1  112000 

6631  FORMAT!*  CONO  CODE  66311  REAL  SOLUTION  TO  FILLET/ROUNO  PROBLEM  OOE  112900 

IS  NOT  EXIST.  POSSIBLE  PROGRAM  ERROR.  SEGMENTS* .213)  113000 

MR1TE(1W.2035)  113100 

STOP  6631  113200 

6636  ASaO.  113300 

6635  A5aS0RTIA5l  113600 

«la-A6*«S  113500 

A2a-A6»AS  113600 

A3aMiaAl*Bl  113700 

A6aMiaA2*BI  113800 

6636  ASaSQRT((X(I)-Al>aa2*(Z(l)oA3)aa2)  113900 

A6aS0RT ( (X ( 1 ) -A2I aaZ* (Z I I I -A6) aaZ)  116000 

IF(A6.LT.AS)60  TO  6660  116100 

KSlinaAl  116200 

ETAII)aA3  116300 

GO  TO  6650  116600 

6660  KSI(I)aA2  116500 

OS  ETAII|aA6  1166O0 

00  GO  TO  6650  116700 

6665  Ala(H(IPl)aa2.H(I)aa2.Rjaa2*Riaa2)/2./(H(IPll>H(t))  116B00 

AZaAl  116900 

A6aRlaa2>IAl>H(II )aaZ  IISOOO 

IF(A6.6E.O.)00  TO  6668  IISIOO 

tF(A8S(A6) .GT.ERRHAXIGO  TO  6666  115200 

A6a0.  115300 

6668  A6a50RT(A61  115600 

A3aK(I>-A6  115500 

A6aK(l).A6  115600 

GO  TO  6636  115700 

6666  WRITE(lHt6667)ItlPl  115800 

6667  FORMATI*  CONO  CODE  66671  FOR  SCGNENTS*t2I3.*  K(I)  a KIIPIII  REAL  S 115900 

lOLUTION  TO  FILLET/ROUNO  PROBLEM  DOES  NOT  EXIST.  POSSIBLE  PROGRAM  E 116000 
ERROR.*)  116100 

MRITE(IH.2035)  116200 

STOP  6667  116300 

6650  DO  6660  U-UIPl  116600 

AlaSORT!  <KSI(I)-H(J)  )*a2*'<ETAII)-K(J)  )aa2)  116500 

A2aH(J)*R(J)a|KSllI)-H(J))/Al  116600 

A3aK<J)*R(3)a(ETA(I)-K(J))/Al  116700 

A6a.ATAN ( ( A2-H ( J) ) / I A3-K { J) ) ) /PlOl 80  1 16800 

IFIJ.EO.lPDGO  TO  6655  116900 

Xl<I)aA2  117000 

Zl«I)aA3  117100 

Tl(I)aA6  IITZOO 

GO  TO  6660  117300 

6655  X2(I)aA2  117600 

Z2IIlaA3  117500 

T2(I)aA6  117600 
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ON 

NO 


6660  CONTINUE 
6600  CONTINUE 

IF  I tSKlP.EO. 01  RETURN 

C PRINT  RESULTS  OF  FILLCT/ROUND  COMPUTATION 
HRITE(IUt6810l 

6810  FORMAT! 'OFILLET/ROUND  EQUATIONS  AND  INTERSECTIONS! •/ 

I'O'i  TXt«  X 'tlOXf 'KSIMlOXt'ETAtflOXf tRHO'dOXf • Xl««10Xii  ZPf 
2 lOXi'  Tl'tlOXt'  XZ'tlOXf*  ZZ'tlOXt*  T2'/ 

3'  •tl3ZI'*MI 

MRITEIlHtOaiS)  <I«XIIitKSIIIItETA(I)tRHOnitXl«II«Zl(IltTllItt 
I X2II»<ZZ(I).T2IIIilBltlO> 

6815  FORNATdOO  •tI2ilPI0E13.5/)) 

C INSERT  FILLET/ROUND  DATA  INTO  OUlCK  ARRAYS 

6900  CONTINUE 
JaO 

DO  6920  I-ltlO 
JBJ*I 

IFIRHOID.OT.RINFIOO  TO  6920 
JaJ*l 

NSEGaNSEO*! 

NSEGPlBNSEO*! 

IFINSC6.LE. 10)60  TO  6905 
WRlTE(l«t6901)I 

6901  FORMATI«  CONO  CODE  69011  INSERTION  OF  ROUND'«13t'  EXCEEDS  DIHEN5IO 
INS  110)  OF  QUICK') 

HRITE(IHt203S) 

STOP  6901 
6905  Lall 
6910  LaL«l 
LNIaL’l 

lFIL>Le.J)GO  TO  6915 
XOIDaXQILNl) 

ZOIDaZOILHl) 

TOIL'DaTOILMltl) 

Ta(Li2)aT0ILMlt2) 

CURVE  ID aCURVE ILMl ) 

GO  TO  6910 
6915  XQIJ)aXl|I) 

ZQIJ)aZl|I) 

TQIJtl)aTl(I) 

TOIJt2)aT2(ll 

TQ(J>li2)aTllI) 

CURVE  III*  1 > aCURVE  lil) 

XaiJ*l)aX2III 

ZaiJ*l)aZ2III 

TOIJ*lfl)aT2II) 

CURVE  III)  aELLX 
6920  CONTINUE 
C PRINT  RESULTS 

MRITEIIHt6930) 

6930  FORMATI'OFINAL  QUICK  ARRAY!') 

MRlTEIlHtTOlO) 

00  6940  laliio 

6940  HRlTEIIWi9000B)l*XOII)fZOIl)tTOII«2)«TQIlil) 

RETURN 

2035  FORMATt'OEXECUTION  HALTED') 


IITTOO 

117800 

117900 

118000 

118100 

118200 

116300 

118400 

118500 

110600 

118700 

118800 

118900 

119000 

119100 

119200 

119300 

119400 

119500 

119600 

119700 

119800 

119900 

120000 

120100 

120200 

120300 

120400 

120S00 

120600 

120700 

120800 

120900 

121000 

121100 

121200 

121300 

121400 

121500 

121600 

121700 

121800 

121900 

122000 

122100 

122200 

122300 

122400 

122500 

122600 

122700 

122000 

122900 

123000 

123100 

123200 
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7010  FORMOTCO'aO  Xi 'X*  • 13Xi 'Z' ilZXi 'TA*  ilZXt 'TF'/'  •fSBC-')) 
9000a  FORMAT!'  'tlZflPOElA.Ol 
END 


SUBROUTINE  OUCK 

C.— 

C SET  UP  QUICK  MODEL 

COMMON  /RA/XC10>tZllOMTMI10>iH(10l«K(10)<R<10l«ailOI  fMdOl* 

1 CURVEI10)iDIR(10ltXlilO>iXZ(10)<ZlllOI«ZZ(10)iTl(IO)>TZI10li 
Z KSI<I0><ETA(10>«RHOIlOI>XQ!ll>>ZOI|l>iTaill>Z)iKRVI3l>CL<Zli 
3 KVIZ)*Vll0t3t3ifSLC3)tRHOSLI10i3l«CRVSLI10t3l 

*fVB(IO«3«Z>a) .CRVB  (lOtB) tXMAPAXllOl iZMAPAX(lO) tXSIDIlO) tZSlDClOl • 
S BLBiai 

COMMON  /lA/lCASEIIO)  tUCASEOOl  tNCASECTI  tNPOl  iNS!3)  iNBIBI 
COMMON  /RS/KItNItKJ«NJtHlJiLlNEiBLNKtRINF«ERRNAX*OPON>ELLXt 
I YUCLfZEROiFOREiAFTEtTOPStBOTStOPUPiOOlViPIOlBOiXJtZJfXItZNHIt 
Z RItai*THI«ejtAl«AZfA3«UI>ltTHJtStTHIJtSIJfAl»HJ,RJfAA*0X 
3 I T VPE • AS VM • UPFO • UP AF f UP T 0 • UPBO • DNFO f ONAF f ONTO t ONBO 
COMMON  /IS/lRtlH*IF.|Pf IPRINTfNSEOfNSEGPltlSTOPf JSTOP>KSTOP< 

1 WJtLtIPliLMlf  IltNPTS>IZiISL<LB>N>NMltI3 
COMMON  /LA/SLOPEIIOI iIDENCBOl iCaRDIBO) 

COMMON  /LS/TltT 

REAL  KtMtKSI (Kl iNI iKJtMJfNIJiLlNEiKVtKRV 
LOOICAL*!  SL0PE>I0EN«CARD'T1*T 
BOOO  CONTINUE 

IFdF.NE.IPIREHiND  IF 
C CROSS  SECTION  MODELS 

MRITEdF.OOOOdlOEN 
lFdSL*NE.OIGO  TO  0001 
HRlTEdF«9000Z) 

MRITE(IFi90003l 
GO  TO  BOOZ 
BOOl  MRlTEdFtZlTOO) 
aOOZ  CONTINUE 
llBl 

HRITE (IF t9000A) t 1 • I 1 tXO d I <X0 INSE6PI I 
C BODY  LINE  MODELS 

MRITElIFiOOOOSIVUCL 

HRITEdFt90006>IlfKRVd)iKVd) 

WRITE dFf9000T)ZEROtZCRO>XOINSE6Pl)tZERO(ZERO>ZCRO 
WRITE  dFi90009) 

DO  B030  U>I«Z 

WRITEdF<9000SiCLIJi 

DO  aOlO  laltNSEG 

IPIB|*1 

I|wl 

IF (CURVE ( I ) .EO.ELLX) I laZ 
WRITE ( IF>90006I I (CURVE ( I I (KV (I 1 I 
60  TO  (BOOBfBOOSItll 
BOOS  CONTINUE 

NRITEdFfOOOOTIXOdltZadliXOdPDtZOdPdtTOdtDtTOdfZI 
GO  TO  BOlO 

8006  WRITEdFtOOOUIXOdlfZUdltXOdPlIiZOdPlI 
BOlO  CONTINUE 

WRITE (IF (900091 


1Z3300 

1Z3600 

1Z3SOO 

123600 

123700 

123S00 

123900 

124000 

124100 

124Z00 

124300 

124400 

124S00 

124600 

124700 

124B00 

124900 

IZSOOO 

IZSIOO 

125200 

125300 

125400 

125S00 

125600 

125700 

125000 

125900 

126000 

126100 

126200 

126300 

126400 

126500 

126600 

126700 

126B00 

126900 

127000 

127100 

127200 

127300 

127400 

127500 

127600 

127700 

127B00 

127900 

12B000 

120100 

12B200 

120300 

12B400 

12B500 

12B600 

12B700 

12BB00 
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DO  aOZO  laltNSEOPl 

ZOCIIa-ZOlII 

TOCItlla«Tg(I«l) 

TOdiZIa-TDIItZI 
BOZO  CONTINUE 
0030  CONTINUE 

HRITEIlFtOOOlO) 

1FIISL<EO.O)60  TO  B033 
DO  B032  La|.3 

MRITEdFfOOOOBISLIL) 

IlaNStLI 

DO  8031  Ial,ll 

IZal 

IF ICRVSL  d (LI .EO.ELLXI IZaZ 

WRITE ( IFtOOOOB) I tCRVSL  d tLI tKV ( IZ) 

WRITE  dF»9000T>VdiltLltVd«ZtL)«Vd*ltliLltVd*l«ZiL>»Vd«3>LI> 
I yd*lt3tLI 
8031  CONTINUE 

WRITE  dFf  900091 
803Z  CONTINUE 

WRITE  dFtZ1900l 
8033  CONTINUE 

WRITE  dF«900ll> 

IFdPRlNT.EO.OlOO  TO  8080 

OXa.OOOOl 

NPTSaS 

C SET  UP  CAROS  TO  EXERCISE  MOUEL 
Hal 
IZa2 

DO  80«0  IbIiNSEG 
AlaXOdlaOX 
AZaXOd*ll«OX 
' A3a|A2-All/(NPTS*ll 
AZaXOdal)*OX 

WRITE dFi9001ZIIl*IZfAltAZfA3 
BOAO  CONTINUE 

IFdSL.EO.OeO  TO  0050 

IlaZ 

IZaZ 

lAaNPIl) 

A3aZ. 

A«a-90. 

A5a90. 

AOaS, 

00  aOAA  Lai ,3 
ISaNPILI 
00  8043  Ia2,I3 
1FIL.EO.1IGO  TO  8042 
IFd.EO.JlOO  TO  8043 
DO  8041  Jal.IA 

IF(ABS(VIJtl.ll-Vd.l.LM.LE.ERRHAX)GO  TO  8043 

8041  CONTINUE 

8042  AlaVd.l.LI 
A2aAl*l. 

WRITEdF.90012IIl.I2.Al.A2.A3.A4.AS.A6 

8043  CONTINUE 


128900 

129000 

129100 

129200 

129300 

129400 

129500 

129600 

129T00 

129800 

129900 

130000 

130100 

130200 

130300 

130400 

130500 

130600 

130T00 

130800 

130900 

131000 

131100 

131200 

131300 

131400 

131500 

131600 

131700 

131800 

131900 

132000 

132100 

132200 

132300 

132400 

132500 

132600 

132700 

132800 

132900 

133000 

133100 

133200 

133300 

133400 

133500 

133600 

133700 

133800 

133900 

134000 

134100 

134200 

134300 

134400 
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80*4  CONTINUE 
SOSO  WRITEdFtOOOia; 

RETURN 

C.— 

C QUICK  INPUT  FORMATS 

c.— 

90001  FORMAT  IBOAU 

90002  FORMAT  I*  !•/•  1 2') 

90003  FORMAT (ILELL  lELLI  PlEC  BLCL  SIO  LSCP*/ 

1 'UELL  2ELL1  PlEC  SID  TUa  U5CPV 

2 • 1 MAPAX') 

9000*  FORMAT <2l2t6Xi2F10<5) 

90005  FORMAT  I A*l 

90006  FORMATI12<lXiA*t'  PlEC  i.A*) 

90007  FORMATI3FI0.S>2(F4.5><AMfF|0.5) 

9000B  FORMAT!'  ••I2|1P9E1*.6> 

90009  FORMAT<'-l'l 


90010  FORMAT !'VLCL 

YUCL'/ 

1 

'VMAPAX 

YUCL'/ 

2 

'ZMAPAX 

YUCL'/ 

3 

'ZSID 

YUCL'/ 

* 

'TUSCP 

ZUCL'/ 

5 

'ZUSCP 

ZUCL'/ 

6 

'YLSCP 

ZUCL'/ 

7 

'YSIO 

ZUCL'/ 

8 

'ZLSCP 

ZLCL') 

90011  FORMAT  I • •) 

90012  FORMATII2tl3t5X.6F10.5) 

90013  FORMAT!'  O'l 
9001*  format !3F10.S«F9.5) 

21700  FORMAT!'  I'tTSX/ 

1'  1 6 NONCIRCULAR  CROSS  SECTION 

SIO 
TUCL 
LSLCP 

SSL 
USLCP 

TUSL 


SUBROUTINE  LIST 
C— 

C LIST  QUICK  INPUT  FILE 
C 


2'LELL 

lELLl 

PlEC 

BLCL 

3'UELL 

ZELL  I 

PlEC 

SID 

*'LBSL 

3LINE 

PlEC 

BLSL 

SLELL'i 

lix/ 

6'LSSL 

*L1NE 

PlEC 

LSLCP 

7'USSL 

5LINE 

PlEC 

SSL 

8UELL'i 

lix/ 

9'UT5L 

6LINE 

PlEC 

USLCP 

A 

' 1 

MAPAX ') 

• FORMAT !'YUSL 

YUCL' 

• 66X/ 

1 

'YLSL 

YUCL*f66X/ 

6 

'YUSLCP 

YSSL*f66X/ 

7 

'YLSLCP 

YSSL' 

• 66X/ 

A 

'ZUSLCP 

ZUSL' 

t66X/ 

8 

'ZLSLCP 

ZLSL' 

t66X/ 

C 

'ZSSL 

ZStO' 

• 66X1 

•t21X/ 

LSCP  'tElX/ 

USCP  'tElX/ 


LELL'f21X/ 


UELL'f21X/ 


COMMON  /RA/X(10)fZll0)>THI10)tH|10)tKll0>tR!10)tBll0><M!10lf 
1 CURVE!10)tDIRI10)<Xl(lU)(X2(10).ZlllOlfZ2(10)>Tl(10)tT2llO>t 


13*500 

13*600 

13*700 

13*800 

13*900 

135000 

135100 

135200 

135300 

135*00 

135500 

135600 

135700 

135800 

135900 

136000 

136100 

136200 

136300 

136*00 

136500 

136600 

136700 

136800 

136900 

137000 

137100 

137200 

137300 

137*00 

137500 

137600 

137700 

137800 

137900 

13BOOO 

138100 

138200 

13B300 

138*00 

138500 

138600 

138700 

138800 

138900 

139000 

139100 

139200 

139300 

139*00 

139500 

139600 

139700 

139800 

139900 

1*0000 
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2 KSllI0)<ETA<10itRH0ll0)tXail|l,Z0llll<TO(ll«2ttKRV|3I.CL(Z)<  IMIOO 

3 KVI2l«VltO«3.3>.SLI3>*RHOSL(10.3).CRVSL(IOi3>  140200 

4*VBI10*3i2«S>«CRVB  (lOfOl  tXMAPAX(IO)  fZMAPAX  (10)  tXSIOdOl  tZSlDllOl  • 140300 

B BLB(B)  140400 

COHNON  /IA/ICASEI10>«JCASE(30)«NCASE(r>iNPI3)(NS(3liNBiai  140S00 

COMMON  /RS/KltMUKJiMJtMlJtLlNEteLNKtRINrtERRMAXtOPONiELLXt  140000 

1 YUCLtZEROiFOREiAFTEtTOPStBOTSiOPUP.OOlVtPlOlBOtXJfZJiXIiZItHli  140700 

2 RI(BI«THIfBJtAltA2fA3>HlJtTNJtSiTHIJtStJiAIfHJtRJtA4tDX  140000 

3*TTPC«ASrMiUPF0fUPAFtUPT0iUPB0tDNF0«DNAFtDNT0»0NB0  140000 

COMMON  /lS/IR*lH*IF.IP<lPRINTtNSEOtNSE6PltlSTOPtJSTOP(KSTOPt  141000 

1 ItJfL<IPltLMlttliNPTS(12tISL«LBtNtNNltI3  141100 

COHNON  /LA/SLOPC(10)«tOCN(aO>tCARO(OOI  141200 

COMMON  /LS/TIiT  141300 

REAL  KtMiKSItRIiMIiKJtNJiMIJtLINEtKVtKRV  141400 

LOBICAL*!  SLOPEiIDENiCAMOtTl.T  141500 

IFdF.EO.lPlOO  TO  OOAO  141000 

9000  CONTINUE  141700 

REMIND  IF  I41BOO 

MRITEdM»9910)lF  141900 

9010  FORMAT  COL  I ST  INO  OF  QUICK  INPUT  DECK  ON  UNIT  •tI2t'CI  142000 

HRITEdHt200TI  142100 

9020  RCADdFt900OltENOa9040ll0EN  142200 

HRITEdHt20O9)IOEN  142300 

9030  FORMAT C •iSOAll  142400 

00  TO  9020  142500 

-O  9040  MRITEdM«200T>  142600 

M 9050  FORMAT!'  ••132C>'II  142700 

HRITEdHt90S5)  142000 

9055  FORMAT  Cl  •)  142900 

9060  RETURN  143000 

90001  FORMAT laOAl I 143100 

2007  FORMAT C '.OZC-'I)  143200 

2009  FORMAT C C.SOAli'CI  143300 

END  143400 

143500 

SUBROUTINE  SLCE  143600 

C— — 143700 
C COMPUTE  SLICE  GEOMETRY  143000 

C— .......  143900 

common  /RA/X(10>iZdOI<THdO>>MdO)iKdO)tRdO)tBdO>>MdOlt  144000 

1 CURVCdOI<OIRdO>iXldUltX2dOI«ZldO)>Z2dO)iTldO)fT2dO)»  144100 

2 KSId0)fETAI10l>RHO(10ltXa(ll)«ZO(MliTOdli2)tKRV|3|«CL(2»>  144200 

3 KV(2ltVdOf3t3)>SL(3l<HHOSLdOt3)>CRVSLdO«3)  144300 

4tVB(10t3i2«a)  tCRVB  dOtOI  tXNAPAXIlO)  (ZMAPAXIIO)  tXSIOdO)  •ZSIDIIO)  • 144400 

S BLBIO)  144500 

COMMON  /lA/ICASECIO) t JCASE(30) iNCASECT) •NP(3I  tNS(3> •NB(B)  144600 

COMMON  /RS/KI»MI»KJ.HU'MIUtLINEtBLNKtRtNFtERRMAX«OPDN<CLLXt  144700 

1 YUCLiZEROfFOREiAFTEtTOPStBOTStOPUPiOGIViPlOIBOtXJtZJtXIfZltHIi  144000 

2 RItBI>THItajfAltA2.A3tBIJtTHJ»SfTHIJtSIJtAItHJ,RJ«A4«DX  144900 

3*TYPEiASTMiUPFOiUPAFiUPTOtUPBO«ONFO<DNAFiDNTOtDNBO  145000 

COMMON  /IS/lRtlMi|FiIPilPNINT.NSEG.NSEGPltlSTOP«JSTOPfKSTOPt  145100 

I ItJtLdPULMlf IliNPTStl2<lSLtLBtN«NM|iI3  145200 

COMMON  /LA/SLOPEdO)tIOEN(ao>»CARD(aOI  145300  ' 

COMMON  /LS/TItT  145400 

REAL  K*H.KSl>KI«HItKJ.MJtMlJ,LlNEtKV«KRV  145500 

logical*!  SLOPEtIDEN.CARDtTIiT  145600 
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10000  CONTINUE 

IF(TYPEtNE.ASYM|00  TO  10006 
C SET  UP  XO  ANO  ZO  for  ASYH  CASE 
NSEOPlst 

XQIl)«VB(N8(ll*l<ltl<l) 

AIbO, 

DO  lOOOS  LaliS 
11bNB<L>*1 
DO  10005  laUIl 
AlaAHAXKAl.VBIItZtltD) 

10005  CONTINUE 

zoniaAi 

10006  CONTINUE 

C DETERNINE  IF  SLICES  ARC  PRESENT 
tSLaO 

DO  10010  lal.lO 
DO  10010  Jal,3 
00  10010  LaltJ 

IFIVdiJtLliGT.RINFlOO  TO  10010 
ISLal 

00  TO  10015 
10010  CONTINUE 

XFITYPE.EO.ASVH.AND.ISLiEO.OIGO  to  10015 
lFIISL>Ea.O|RETURN 

C DETERMINE  NUMBER  OF  SEGMENTS  IN  EACH  SLICE  NSILiaNP(L)-l 
10015  00  10020  LaltS 
NPILIaO 

DO  10020  laltlO 
00  10020  Jal.3 
1F(V(1«J.L) .LT.RlNF)NPIL)at 
10020  CONTINUE 

00  10030  Lalt3 
NS(LlaNP(L>«l 
10030  CONTINUE 

C CHECK  INPUT  FOR  SUFFICIENCY 
DO  10060  Lai. 3 
llBNPIL) 

IFdl.CO.OlGO  TO  10060 
DO  10055  lal.Il 
12a0 

DO  10060  Jal,3 
lF(VdtJ*LI.LT.RlNF>12aI2*l 
10060  CONTINUE 

IFd2.GE.2IG0  TO  10055 
ISTOPal 

HRITEIlH.lOOSOlLtl 

10050  FORMAT!'  COND  CODE  OOSOl  SLICE'. 12. ••  POINT'. 13.'.  DATA  INSUFFICIE 
INT'l 

lOOSS  CONTINUE 
10060  CONTINUE 

IFdSTOP.EO.OlGO  TO  lOOTO 
WRITE  I IM. 20351 

2035  FORMAT!'  EXECUTION  HALTED' I 
STOP  0050 

C COMPUTE  REMAINING  SLICE  GEOMETRY 
lOOTO  00  10160  Lai. 3 


16ST00 

165800 

165000 

166000 

166100 

166200 

166300 

166600 

166500 

166600 

166700 

166800 

166900 

16T000 

16T100 

16T200 

16T300 

16T600 

167500 

16T600 

167700 

16TBOO 

167900 

160000 

168100 

168200 

168300 

160600 

168500 

168600 

168700 

160800 

168900 

169000 

169100 

169200 

169300 

169600 

169500 

169600 

169700 

169800 

169900 

150000 

150100 

150200 

150300 

150600 

150500 

150600 

150700 

150800 

150900 

ISIOOO 

151100 

151200 
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II«NPIL)  151300 

IF(Il«Ea.0IG0  TO  10160  151600 

00  10150  U2.I1  151500 

IHUUl  151600 

DO  10000  Jal.3  151700 

1PIV(I.J.L>.6T.R1NF)00  70(10100,10130.10160) .J  151800 

10080  CONTINUE  151900 

«la(VII.2.L>-V(IM1.2,L)>/(V(I.l.LI-V<IMl.l.L>)  152000 

«1bATAN(A1)/PI0180  152100 

lF(ABS(Al-Vlt.3,LI).LE.ERRNAX)00  TO  10150  152200 

ISTOPal  152300 

HRITC(IM.10090)k.I.ERRHAX  152600 

10090  FORMAT!'  COND  CODE  00901  SLICE', 12.'.  POINT'. 13.'.  REDUNDANT  INPUT  152500 
1 DATA  NOT  SELF  CONSISTENT  TO  WITHIN  '.1PE15.7)  152600 

00  TO  10150  152700 

10100  lF(V(1.3,L).EO.0.)OO  TO  lOIlO  152800 

V(I.l.L|a(V(I.2.L)-V(IH1.2.L) )/TAN(P10180*V(I.3.LI )*V(IM1.1.LI  152900 

60  TO  10150  153000 

lOllO  WRITEIIH.10120IL.1  153100 

10120  FORMAT!'  COND  CODE  01201  SLICE'. 12.'.  POINT'. 13.'.  X MUST  BE  INPUT  153200 
1 WHEN  SLOPE  IS  ZERO')  153300 

WRITE  I IH, 2035)  153600 

STOP  0120  153500 

10130  V!I.2.L)aV!IM1.2.L)*!V!I.l.L)>V!lMl.l.L))OTAN!PI0180*V!1.3.L))  153000 

60  TO  10150  153700 

10160  V!I.3.L)aArAN!!V(I.2.L)>V!IMl,2.L))/!V!l.l.L)-V!lNl,l,L)))/PI0180  153000 

^ 10150  CONTINUE  153900 

10160  CONTINUE  156000 

lF(IST0P.Ea.0)60  TO  10170  156100 

WRITE !1W. 2035)  156200 

STOP  0090  156300 

C extrapolate  beginning  OF  FIRST  SEGMENT  TO  XaO  156600 

10170  DO  10180  Lai. 3 156500 

IUNP!L)  156600 

IFdl.EO.OlOO  TO  10180  156700 

IF!V!1.1.L).LE.0.)  go  to  10180  156800 

V!l,2.L)«V!1.2.L)-V!l.l.L)«TAN!P10180aV!2.3.L))  156900 

Vll.l.DaO,  155000 

10180  CONTINUE  155100 

C extrapolate  end  of  last  segment  to  XQINSEOPD  155200 

00  10190  Lai. 3 155300 

llaNPIL)  155600 

IF!11.EO.O)60  TO  10190  155500 

tF!V!ll.l.L).6E.X0!NSE6Pl))G0  TO  10190  155600 

V!I1.2.L)aV|tl,2.L)*lxa(NSEGPl)aV!Il.l.L))*TAN(P101SO>V!I1.3.L))  155700 
V!Il.l.L)axO!NSE6Pl)  155800 

10190  CONTINUE  155900 

C PLACE  NON  EXISTENT  SLICES  OUTSIDE  OF  VEHICLE  156000 

AUO.  156100 

DO  10200  Ul.NSEGPl  156200 

AUAMAXKAl.ABSIZO!!)))  156300 

10200  CONTINUE  156600 

AU2.PA1  156500 

DO  10210  Lal,3  156600 

tF!NPIL).NE.0)60  TO  10210  156700 

NP!L)p2  156800 
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ON 


,k: 


NSCL>al 

VdtltDaO."  - 

V(lt2«L>BAl 

V<2tltL>aX0(NSe6Pll 

V(2i2tL)-Al 

V(l«3tL)aO. 

VI2«3fLlaO. 

10210  CONTINUE 

C BE  SURE  V(lt3tL)  IS  DEFINED 
DO  10220  Lai >3 

lFIV(1.3iLi.GT.RINF|V(lt3iL>aV<2i3tLI 
10220  CONTINUE 

C PRINT  RESULTS  OF  SLICE  6E0METRT  CALCULATIONS 
URITE(lHil0230) 

10230  FORNATCOSLICE  SEOHETRTl  •/•0'«21X>'UPPER'<39Xi«LOHER'«39X<<SIOE  • 
1/  9Xt«X'il3Xi'Z'  (nXi'Tif  15Xi'X«<13Xi'Z>fl3Xt'T>ilSXi*X>fl3Xt 
2'Y»fl3Xt*T»/»  'tlSZl*-')) 

DO  102S0  lal.lO 

NRITE(Uil0ZA.0U.IIV(I«J»L>tJal,3ULali3) 

102AO  FORHATC  ••I2<2(lP3ElA.bt2X)  (IPSEIA.O) 

102SO  CONTINUE 

C CHECK  VALIDITY  OF  COORDINATES 
DO  10280  Lali3 
IlaNPILI 
DO  10270  lalfll 
IF(Vllt2iL».6E.0.)O0  TO  10270 
ISTOPal 

HRITC(IH«10260)L»I 

10200  FORHATC  COND  CODE  02601  SLlCE'flZt'i  P0INT'tI3f't  Z VALUE  < 0>  SL 
IICES  MAY  NOT  CROSS  VEHICLE  CL.M 
10270  CONTINUE 
10280  CONTINUE 

IF(ISTOP.EO.O)GO  TO  10290 
HRITEIlHt2035) 

STOP  0260 

C RESIGN  LSL  FOR  QUICK  INPUT 
10290  llaNP<2l 

IFITYPE.EG.ASYMIIIETURN 
DO  10300  lal.Il 
DO  10300  Ja2<3 
V(I.J.2|a-VII»J«2l 
10300  CONTINUE 
RETURN 
END 

SUBROUTINE  RF5L 

C.....— 

C INSERT  ROUNDS  OR  FILLETS  BEThEEN  SEQUENTIAL  SLICING  PLANES 

COMMON  /RA/XllO)fZIIOI>THaO»iH|10l(K(IOI<RllOI«B(lo7fMllO)< 

1 CURVE(10l*OIR<10)tXlU0)iX2(10I.Zlll0).Z2ll0).Tt(10l«T2(10)i 

2 KSl(10ltCTA(10)«RHOI10>iXQ(lll.ZO(ll)fTO(llt2).KRV(3ltCL(2)i 

3 KV(2ttV<10t3>3)«SLI3)«RMOSL(10»3)tCRVSLUOt3) 

AfVB(IO*3i2*8)  iCRVB  llOiB)  iXNAPAXIlOl  •ZMAPAX(IO)  •XSIDIlOt  tZSIOdOl  > 
5 BLB(BI 

COMMON  /IA/ICASE(10).JCASE(30>>NCASEI7|,NP(3|,nSI3>«NB(8) 


1S6900 

1S7000 

157100 

157200 

157300 

157600 

157S00 

157600 

157700 

157800 

157900 

158000 

158100 

158200 

158300 

158600 

150500 

158600 

158700 

158800 

158900 

159000 

159100 

159200 

159300 

159600 

159500 

159600 

159700 

159800 

159900 

160000 

160100 

160200 

160300 

160600 

160500 

160600 

160700 

160800 

160900 

161000 

161100 

161200 

161300 

161600 

161500 

161600 

161700 

161800 

161900 

162000 

162100 

162200 

162300 

162600 
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common  /il5/KItMI,KJ<MJiMIJ.LINE<eLNK<RINF>ERRMAXfOPON«ELLXt 

1 YUCUtZEROtPORE>AFTCfTOPStBOTStOPUP«06IV(PI01BO,XJfZJ<XItZliHIi 

2 RIiBItTHItBJtAUAZtAStBIJtTHJtStTHIJiSIJtAIiHJ.RJtAAtDX 
StTYPEtASVMiUPFOtUPAFfUPTOiUPBOtONFOiONAFtONTOtONBO 

COMMON  /IS/IR«IW<IFilPilPRINT«NSEO<NSEGPl«lSTOP«JSTOP<KSTOPt 
I liJtLtlPltLMltllfNPTStlZtlSLtLBtNtNNltia 
COMMON  /LA/SLOPE  1 1 0 1 1 1 OEN ( BO 1 1 C ARD  < BO I 
COMMON  /LS/TIiT 

REAU  KtNiKSIiKIfNItKJtHJfMtJtLlNE»KVtKRV 
LOGICAL*!  SLOPEtlDENtCARDiTI.T 
11000  CONTINUE 
tSKIPaO 

IFIISL.EO.OIRETURN 

c determine  if  rounds  or  fillets  are  present 

IRFSL  *0 
DO  llOOS  L*lt3 
DO  11005  lal.lO 

IFIRNOSLdiLl.OT.RINFtGO  TO  llOOS 

1RFSL*1 

GO  TO  11006 


11005  CONTINUE 
RETURN 

11006  HRITEIlHdlOOTI 

HOOT  FORMAT  COSLICE  ROUND/FILLET  CENTERS  AND  INTUSECTIONSI  t 

1 /'O'*  9Xi'BIttllXi<BJitllXt(A3itllX«</U>tllXfiXl>tllX»«Zlt, 

1 llXt'Tl'illXi'XZi.llXt'ZZitllXt'TZ'/*  •. 

2 132('-'i) 

C COMPUTE  ROUNO/FILLET  CENTERS  AND  INTERSECTIONS 
DO  11130  Lalt3 
DO  11030  I*2t9 
1P1-I*1 

lFlRHOSL(lfLl«OT.RINF>00  TO  11030 
ISKIP-l 

IF(ABS(V(Ii3tLI-V<IPU3tL)l*OE«ERRHAXIOO  TO  11020 
WRITEIlW.llOlOINIPl 


11010  FORMAT!'  COND  CODE  10101  FILLETS/ROUNDS  MAY  NOT 
IN  LINE  SEGMENTS  WITH  SAME  SLOPEi  SEGMENTS' t2I3l 


BE  INSERTED  BETHEE 


WRITE  UHi2035i 


2035  FORMAT!'  EXECUTION  HALTED'! 


STOP  1010 


11020  Salt 


lF!VllPli3tLI>0T.V!If3>LIISa-l. 
NIaTAN!PI01SO*V!t«3iLI ) 
HJ*TAN!PI0ia0*V!tP1.3tLii 
AlBV!t.2«L!-MI*YIIfl.Li 
62aV!lPlt2«LI>MJ*V!IPliltLI 
BIaAl«S*RHOSL ! I iLI /COS !P10160*V ! I t S'Ll! 
BJaA2>S*RH0SL I I (LI /COS !PIOiao*V I IPl tStLI I 
A3a!MJ*BIaMI*BJ>/!MJ-Hli 
A4a- !BJ-BI I / INJ-MI I 

Xl!I|aA6-$*RHOSLIItLI*SlN!P101SO*V!It3«LII 
X2 ! I I ba«-S*RHOSL I I tLI •sin !P10iaO*V ! IPl (StLI ) 
Z1 ! I )aA3*S*RH0SL i I tLI *00$ !P10iao*V ! I i3tL 1 1 
22 ! I I aaatSvRHOSL ! I tLI pCOS !P10180*V I IPl t3 tLI I 

Tl!I|aV(|,3,|.| 

T2lllaV!IPlt3tLI 


162500 

162600 

162700 

162600 

162900 

163000 

163100 

163200 

163300 

163400 

163500 

163600 

163700 

163000 

163900 

164000 

164100 

164200 

164300 

164400 

164500 

164600 

164T00 

164800 

164900 

165000 

165100 

165200 

165300 

165400 

165500 

165600 

16ST00 

165000 

165900 

166000 

166100 

166200 

166300 

166400 

166S00 

166600 

166T00 

166000 

166900 

16T000 

16T100 

16T200 

16T300 

16T400 

167500 

16T600 

16TTOO 

167000 

16T900 

160000 
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'j 

00 


HRITEIlHill025)L>ItBltBJt*3tA«fXl(t)*ZI<n<TllI)iX2(I)>Z2II)«T2Cl) 
11025  F0RM*T(lXi2lltlP10E13.5l 
11030  CONTINUE 

C INSERT  ROUNDS/FILLETS  INTO  UUICK  ARRAYS 
DO  11035  laltO 
tPlat*! 

RHOSL ( I iL) aRHOSL ( IPl (L ) 

XllllaXlllPlI 

ZIIII-ZKIPII 

T1II>-T1«IP1) 

XZIllaXZlIPlI 

Z2(l|aZ2<IPl) 

T2CllaT2<IPl) 

11035  CONTINUE 

RH0SLI10f|.>al«ET0 

JmO 

00  11110  laltlO 

JaJ*l 

IFIRHOSLIItLI.OT.RlNFIGO  TO  11110 
JaJ*l 

NS(LlaNSILI*l 

NP(L>aNPa>*l 

IF(NP(LI.LE. 10160  TO  11050 
HRITEIiNtllOAO)! 

IIOAO  FORMATC  CONO  CODE  lOAOl  INSERTION  OF  R0UND/FIUET'tI3t'  EXCEEOS  0 
IIHENSIONS  HO)*) 

HRITEIINi2035l 
STOP  10 AO 
11050  Nall 
11060  NaN'l 

NN1bN-1 

IFIN.LE.JIGO  TO  11080 
DO  llOTO  I2al,3 
11070  VINi|2tL)aVINM1.12tL> 

CRVSLINtD-CRVSLINMltLi 
60  TO  11060 
11080  CONTINUE 

V(J>l>L)aXl(I) 

V(Ji2tL)aZl(I> 

V(Ji3iL)aTlll> 

CRVSL(J*ltL>aCRVSLIJ»L> 

V<J*liltLlaX2(ll 

VIJ*li2fLlaZ2<l> 

VIJ*l>3iL)aT2(l> 

CRVSLIJtLiaELLX 
11110  CONTINUE 

HRITEIIatlllZOI 
11120  FORMAT  I • •! 

Ill JO  CONTINUE 
C PRINT  FINAL  QUICK  ARRAYS 
mRITEIIH«I0230I 

10230  format  I 'OF INAL  QUICK  ARRAYS  FOR  ROUNDEO/FILLETTEO  SLICES! •• 

1 /'O'tZlXt'UPPER'tSOXiiLOWER'iSOXf 'SIDE  • 

1/  9X>'X!«l3Xt'Z«  fl3X.'T*<lSXt'X>il3X«'Zi.l3Xt!T'fl5Xi'X«>13Xt 
2<Y'tl3Xi'T'/!  ••132(>-'ll 
DO  10250  tBl,10 


168100 

168200 

16B300 

16BAOO 

168500 

168600 

168700 

168800 

168900 

169000 

169100 

169200 

169300 

169A00 

169500 

169600 

169700 

169BOO 

169900 

ITOOOO 

170100 

170200 

170300 

170A0O 

170500 

170600 

170700 

170800 

170000 

171000 

171100 

171200 

171300 

171A00 

171500 

171600 

171700 

171800 

171900 

172000 

172100 

172200 

172300 

1T2A00 

172500 

172600 

172700 

1T2B00 

172900 

173000 

173100 

173200 

173300 

1T3AOO 

173500 

173600 
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-j 

VO 


HRlTE(lH>10240II>((V(I«J<LI>Jal.3)»Lal>3l 
10240  FORMATC  <tl2i2(tR3E14.6<2XI  tlPSElO.OI 
10250  CONTINUE 
RETURN 
END 

SUBROUTINE  COPY 

C--...-.— .... ....... ..... 

C COPY  ORTA  INTO  VB  ARRAY  FOR  0K30 

C............... ........... 

COMMON  /RA/Xll0i«Z(10l<TH(10)tH(tO).Ka0liRllOI.BI10t.M(lO)t 

1 CURVe(10)fOIR(10)tXl<10l<X2ll0>tZI(IOIf22ll0i.Tlll0»iT2ll0ti 

2 KSin0)tETA<10liRHO(lO)«Xailll,ZQIlll.Ta(lli2).KRV(3)fCL(2l« 

3 KV(2) tVlt0i3t3l«SL(3>«R -OSLIlOiSl tCRVSL(10t3l 

4«VB(14*3i2»B)>CRVB  llOtOt  (XHAPAX(IO)  tZHAPAXdOi  .XSIOIIO)  .ZSIOdOl  t 
5 BLBIS) 

COMMON  /lA/ICASEdO)  tJCASEOO)  fNCASElT)  tNP(3l  iNSCSl  tNBlBI 
COMMON  /RS/Kt«NItKJ»MJtMlJ<LINEiBI.NK<RINFtCRRMAX«OPON«ELLXt 

1 YUCLiZEROtFOREtAFTEtTOPSfBOTSiOPUPfOGIVfPlOlOOtXJtZJtXItZItHIt 

2 RI>BltTHI.BJtAl«A2tA3tBIUfTHJiSiTHIJtSIJtAliHJtRJtA4(0X 
3tTYPE«ASYMtUPF0«UPAF«UPT0iUPBOi0NF0<0NAFt0NT0>DN60 

COMMON  /IS/IRt IM< IF t IP< IPRINT tNSEGtNSEBPl . ISTOP<USTOP>KSTOP> 

I I(JiLtIPlfLM|(Il«NPTS>I2tISLtLB<NtNMl(I3 
COMMON  /LA/SLOPE  d 0 > t IDEM ( GO  > t CARD I BO ) 

COMMON  /LS/TItT 

REAL  K.M.KSItKItNItKJfMJ.MIJtLlNEtKViKRV 
LOGICAL*!  SLOPEiIOENfCARDiTItT 
dimension  EClOiTI 
EQUIVALENCE  (EdflltX  dll 

1 t(Edi2>iZ  dM 

2 *(Edt3)tTMd») 

3 *IEdt4liH  dl) 

4 tlEdf5)tK  d)l 

5 *(Ed<At«R  d>) 

A >(EdiT)fB  d>) 

DATA  RlNFl/l.ETO/ 

LOGICAL*!  F/tpi/ 

12000  CONTINUE 

C 1234S6TB 

GO  TOd20!0»!20!Of!20!Of!20SOi!2050i!20SOt!2080t!20BO 
II  (LB 

C LB-l  ZUCLi  LB«2  ZLCLt  LBs3  YSIO 
12010  NBILBI*NSEG 

00  12020  lalilO 

VBdtltlfLBIaXOdl 

VBdt2*ltLB|aZOd» 

VBdt3*l«LBI«Tedtll 
VBdt3*2(LB|aTadt2l 
CRVB  dtLBIsCURVEdI 
12020  CONTINUE 

DO  12040  t-lflO 

XQdlaRINFl 

ZQdlaRINFl 

TOdtllaRlNFl 

TQdtElaRINFl 

DO  12030  JaltT 


1T3T00 

lT3aoo 

173900 

174000 

174100 

174200 

174300 

174400 

174500 

174600 

174700 

174600 

174900 

175000 

175100 

175200 

175300 

175400 

175500 

175600 

175700 

175600 

175900 

176000 

176100 

176200 

176300 

176400 

176500 

176600 

176700 

176600 

176900 

177000 

177100 

177200 

177300 

177400 

177500 

177600 

177700 

177600 

177900 

17B000 

17B100 

17B200 

17B300 

17S400 

I7BS00 

1TB600 

17B700 

17S600 

176900 

179000 

179100 

179200 
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ElItJIaRINFl  1T9300 

12030  CONTINUE  179*00 

SLOPEIDaF  179500 

RHOlIlaRINFl  179600 

XllllaRlNFl  179700 

X2(I|aRINFl  179800 

ZKD-RINFI  179900 

Z2Ct>«RINFl  ISOOOO 

T1CII«RINF1  180100 

T2I1)«RINF1  180200 

KSICllaRINFl'  180300 

ET*(I»bRINF1  180*00 

CURVEIDaBLNK  180SOO 

DIRdl-OPON  180600 

120*0  CONTINUE  180700 

RETURN  180800 

C LBa*  ZUSL>  LBa5  ZLSLt  LBbS  YSSL  180900 

12090  ItaL8-3  181000 

NBILBl-NSfll)  181100 

DO  12070  laltlO  181200 

DO  12060  Jalt3  181300 

VBdtJtliLBlBVdvJiIll  181*00 

CRVBdtLBlaCRVSLdf  111  181500 

12060  CONTINUE  181600 

lFd.LT.10)VBdi3i2fL8lBVd*lt3.Il>  181700 

00  12070  CONTINUE  181800 

O RETURN  181900 

12080  DO  12090  IsltlO  182000 

VBdtldiTlaXSlOdl  182100 

VBd<2>lt7)aZSlDdl  182200 

VBdtltltBlaXMAPAXdi  182300 

VBdt2>li8)«ZHAPAXdl  182*00 

12090  CONTINUE  182500 

RETURN  182600 

END  102700 

C**«*»***«***«*»***»>»»»*«»*«»»*>»******»»»*«**»*»<**»»»*»»  182800 

SUBROUTINE  BLCK  182900 

C..-.. 183000 

C CHECK  ALL  BODY  LINE  COORDINATES  FOR  VALIDITY  183100 

C COMPUTE  default  VALUES  FOR  ZSID  AND  ZHAPAX  183200 

C — — 183300 

COMMON  /RA/XdO)«ZdOliTHdO)tH(lO),KdO)iR(10)<BdO>tMdOI  t 183*00 

1 CURVEdO)«OIRdOliXldO)»X2dOliZldO>iZ2(10lfTldO>«T2dOli  183500 

2 KSldO)  iETAd0)*RH0(10)  iXOdl ) .ZUd  1 1 .TOdl  tZl  .KRV  (3>  .CL  (2l  i 183600 

3 KVI2I.Vd0«3.3ltSLI3).RHOSLT10i31tCRVSLd0.3)  183700 

*.VBdO. 3.2.8)  .CRVB  (10.8)  .XMAPAXdO)  .ZNAPAX(IO)  .XSlDdO)  .ZSIDdO)  . 183800 

5 8LB(8)  183900 

COMMON  /IA/lCASEd0).JCASE(30)  .NCASE(T)  .NP(3)  .NS(3)  .NB(8)  18*000 

COMMON  /RS/KI.MI.KJ.MJ.MIJ.LINE.BLNK.RINF.ERRMAX.OPDN.ELLX.  18*100 

1 YUCL.ZEROiFORE.AFTE. TUPS. SOTS. OPUP.OGIV.P101BO.XJ.ZJ.XI.ZI. Ml.  18*200 

2 HI.BliTHI.BJ.*l.A2.*3.8lJ.THJ.s.THIJ.SIJ.*I.HJ.HJ.A*.DX  18*300 

3.TVPE.A5VM.UPFO.UPAF.UPTO.UPBO.ONFO.DNAP.ONTO.DNBO  18**00 

COMMON  /IS/lR.IH.tF.IP.lPRtNT.NSE6.NSE6Pl.IST0P.JST0P.K5T0P.  18*500 

1 I.J.L.IP1.LM1.11.NPTS.12.ISL.LB.N.NH1.13  18*600 

COMMON  /L*/SLOPEdO).IDEN)BO).CARO(60l  18*700 

COMMON  /LSFTI.T  18*800 
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REAL  KtMtKSI>KI,MItKJtMJiH|J«LINEtKVtKRV  1S4900 

L061C*L*1  SLOPEtlOEN.CAROtTltT  ISSOOO 

13000  CONTINUE  105100 

C ZSID  AND  ZNAPAX  DEFAULT  185Z00 

DO  13060  LaTiO  ISS300 

MB  CL I *0  105600 

DO  13010  laltlO  105500 

IFIVBCIiltl.Ll.LT.RlNF.AND.VBCIiZiltLI.LT.RINFINBCDal  105600 

13010  CONTINUE  105700 

NBILIaNB(L)>l  105000 

IFINBCLI.GT.OIGO  to  13060  IBSBOO 

IFINBCLI.EO.-liOO  TO  13030  106000 

WRlTE(lHtl30Z0IL  106100 

13020  FORHATC  CONO  CODE  30ZOI  BOOT  L|NE*tI2«'  REQUIRES  X AND  Z VALUES  F 106200 
lOR  THO  OR  MORE  POINTS')  106300 

HRITE(lNi2035l  106600 

2035  FORMAT (•  EXECUTION  HALTED')  106500 

STOP  3020  106600 

13030  1FIL«NE.T)G0  TO  13060  106700 

NBIT)b1  106000 

VBIl«l*li7|a0a  106900 

VBIlt2»ltT|aO.  107000 

VBI2>lfl(7)8X0CNSE0Pl)  107100 

VB(2<2tl<7|sO.  107200 

GO  TO  13060  107300 

13060  IlaNBIT)*!  107600 

NBI0)aNBI7)  107500 

DO  13050  lal.Il  107600 

DO  13050  JaliS  107700 

VB(ItJ*l<B)sVB(lf Jtl«7>  107000 

13050  CONTINUE  107900 

13060  CONTINUE  100000 

C VALIDITY  CHECK  100100 

Al8VBINB<l)*ltlil<ll  100200 

00  13000  La2>0  1B0300 

IlaNBCL)*!  100600 

1F(ABS<VB(I1>1(1iL)-A1).LE.ERRMAX)G0  to  13000  100500 

IFIVB<tltl(l»L).eE.AliAND.L«eEi6)60  TO  13000  100600 

MRITEIlHtl3OT0)ERRHAXtL  IBB700 

13070  FORMAT!'  CONO  CODE  30701  ALL  BODY  LINES  MUST  BE  THE  SAME  LENGTH  TO  100800 
1 WITHIN  ERRMAXB'tlPElA.Tt'  BL'iIZ)  100900 

HRITECIH.ZOSB)  109000 

STOP  3070  109100 

13000  CONTINUE  109200 

DO  13100  LaliB  109300 

IFIABS(VB(ltltltL)).LE.CRRMAX)00  TO  13100  109600 

HRITEIlNil3090ICRRMAX<LtVB(l>ltltL)  109500 

13090  FORMAT!'  CONO  CODE  30901  FIRST  X OF  EACH  BOOT  LINE  MUST  BE  ZERO  IT  109600 

10  WITHIN  ERRMAXB',lPE16.7f ')  BL'tIZi'  Xb',1PE16.T)  109700 

WRITE ClWtZOSS)  109800 

STOP  3090  109900 

13100  CONTINUE  190000 

C RESIGN  ZLCL  AND  ZLSL  190100 

00  13110  La2iSi3  190200 

IlaNBIL)*!  190300 

DO  13110  lalill  190600 
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00  13110  jB2t3 
00  13110  NalfZ 

invBdtJtNtD.GTtRINFlOO  TO  13110 
VBdiJiNiLla-VBdtJtN.LI 
13110  CONTINUE 

C SET  UP  ANOLES  FOR  ZHAPAX  AND  ZSID 
00  14026  LaTtS 
IIbNBIL) 

00  14025  Ultll 

AlBATAN(IVBd*l<2iltL)-VBdt2tltLI)/(VBd*ltUl(LI-VBdiltltLlil 

Al-Al/PIOIBO 

VBdiSiltLIaAl 

VBdi3>2tL>-Al 

14025  CONTINUE 

14026  CONTINUE 
RETURN 
END 

SUBROUTINE  0K3D 

C— ....... ........ 

C SET  UP  QUICK  MODEL  FOR  30  CASE 
C............................... 

COMMON  /RA/XdO)tZdO)iTH(tO)fH(10>.K(10)>RdO),B(10l*NdOlt 

1 CURVE  ( 10)  lOIRI  10)  tXl  (10)  iXZdOl  >21  (101  <2200)  iTl  (101  *12(101  • 

2 KSl(10)*ETA(10)*RHO(10)iXO(ll)*ZO(ll)*Ta(tlf2)tKRV(3)*CL(2)t 

3 KV  (2)  •¥(  10.3*3)  «SL  (3)  *R-<OSL(  10.3)  •CRVSL(  10*3) 

4*  VB(  10*3*2*8)  *CRVB  (lO.B)  .XMAPAXdO)  .ZMAPAXdO)  .XSlOdO)  .ZSIDdO)  * 
S BLB(B) 

COMMON  /lA/lCASEdO)  *JCASE(30)  *NCASE(T)  *NP(3)  *NS(3)  *NB(B) 

COMMON  /RS/KI*M1*KJ*MJ*N1J*LINE*BLNK*R1NF*ERHMAX*0PDN.ELLX* 

1 YUCL*ZERO*FORE*AFTE*TOP5«aOTS*OPUP*OOIV*PI01BO*XJ*ZU*Xt*ZI*HI* 

2 R1*B1*THI*BJ*A1,A2*A3*BIJ«THJ*S*THIJ*SIJ*AI*HJ*RJ*A4*DX 
3*TYPE*ASYM*UPFO*UPAF*UPTO*UPBO*DNFO*DNAF*ONTO*DN80 

common  /IS/lR*tW*IF*IP*IPRlNT*NSEG*NSEGPl*lSTOP*JSTOP*KSTOP* 

1 I*J*L*1P1*LM1*11*NPTS*12*ISL*LB*N*NH1*I3 
COMMON  /LA/SLOPE (10) *IOEN (BO) *CARO)BO) 

COMMON  /LS/TI.T 

REAL  K*N*KSI*KI*MI*KJ*HJ*MIJ*LINE*KV*KRV 
LOGICAL*!  SLOPE.IOEN.CARO.TI.T 
14000  CONTINUE 

lFdF.NE.IP)REMIND  IF 
C CS  MODELS 

HRITE(IF*90001>IOEN 

WRITCdFfZITOOl 

11*1 

WRITE (IF.90004) II *I I *ZERO*VB(NB(I ) *1 *1* 1 *1 I 
C SL  MODELS 

WRITE) IF. 9000S)ruCL 
NRITE(IF*90006)I1*KRV(1)*KV(1) 

WRITE(IF*90007>ZERO*ZEHO*VB(NB(1)*1*1*1*1)*ZERO*ZERO*ZERO 
WRITE (IF*90009> 

00  14020  L-1.8 
11*NB(L) 

WRITE  dF*90005)BLB(L) 

00  14010  141.I1 
1P1>I*1 


190SOO 

190600 

190700 

190800 

190900 

191000 

191100 

191200 

191300 

191400 

191500 

191600 

191700 

191800 

191900 

192000 

192100 

192200 

192300 

192400 

192500 

192600 

192700 

192800 

192900 

193000 

193100 

193200 

193300 

193400 

193500 

193600 

193700 

193800 

193900 

194000 

194100 

194200 

194300 

194400 

194500 

194600 

194700 

194B00 

194900 

195000 

195100 

195200 

195300 

195400 

195500 

195600 

195700 

195800 

195900 

196000 
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oo 

OJ 


12-1 

IFICRVB(ItLI.EO.£LI.X)I2>2 

WRITE ( IFf 900061 1 tCRVB lift)  tKV.1 121 

WRITE(IF«9000TlVBII«liliL) tVBIl<2.liL) tVBlIPlf ltltL>(VBIIPIt2il«L> 
ltVB(Ii3il«LI<VBIIt3«2(U) 


14010  CONTINUE 

WRITE(IF<90009I 
14020  CONTINUE 

WRITEIIFt21900) 

WRITE<IFi9001U 
IFdPRINT.EO.OlGO  TO  14030 
C SET  UP  CARDS  TO  MAKE  QUICK  EXERCISE  MOOEU 
OXB.OOl 
AlaOX 

A2«V8(NB(l)*ItltliI)>DX 

A3a(A2*Al>».0I 

A2-A2*2.*0X 

11*1 

I2a2 

WRtTE(IFi90012l IltI2<AlfA2<A3 
14030  WRITE«IF«90013I 
RETURN 

90001  FORMAT IBOAl) 

90004  FORMAT (212. 6Xi2F10.S) 

90005  FORMAT (A4I 

90006  FORNATII2.1X.A4.I  PlEC  '.AAI 
9000T  FORMAT (3F10.5.2(F9.S. 'A* I .FlO.Sl 
90009  FORMAT (••!•) 

90011  FORMAT (•  •) 

90012  FORMAT(I2.I3.SX.OF10.S) 

90013  FORMAT  I • O'l 

90014  format (3F10.5.F9. 51 
21700  FORMAT (•  I'.TBX/ 


1*  1 6 

NONCIRCULAR  CROSS  SECTION 

2>LELL 

lELLI 

PIEC  BLCL 

SID 

LSCP 

3>UELL 

2ELL1 

PIEC  SIO 

TUCL 

USCP 

4'LBSL 

3LINE 

PIEC  BLSL 

LSLCP 

5LELL<. 

lix/ 

O'LSSL 

4L1NE 

PIEC  LSLCP 

SSL 

7>USSL 

SLINE 

PIEC  SSL 

USLCP 

BUELL* . 

llx/ 

9«UTSL 

6LINE 

PIEC  USLCP 

TUSL 

A 

• 1 

NAP  •> 

1 FORMAT  CYUSL 

YUCL'.66X/ 

1 

•TLSL 

YUCL'f66X/ 

2 

•VLCL 

YUCL'.66X/ 

3 

•VMAP 

YUCL'*66X/ 

4 

•YUSCP 

YS1D'.66X/ 

5 

•YLSCP 

YS10«.66X/ 

6 

•YUSLCP 

YSSL*i66X/ 

7 

•YLSLCP 

YSSL'.66X/ 

8 

•ZUSCP 

ZUCL't66X/ 

9 

•ZLSCP 

ZLCL'iOOX/ 

A 

•ZUSLCP 

ZUSL'.66X/ 

8 

•ZLSLCP 

ZLSL>.66X/ 

C 

•ZSSL 

ZSI0't66XI 

•t21X/ 

•.21X/ 

••21X/ 


LELL't21X/ 


UELL'.21X/ 


196100 

196200 

196300 

196400 

196500 

196600 

196700 

196800 

196900 

197000 

197100 

197200 

197300 

197400 

197S00 

197600 

197700 

197800 

197900 

loaooo 

190100 

198200 

198300 

198400 

198500 

198600 

198700 

198800 

198900 

199000 

199100 

199200 

199300 

199400 

199500 

199600 

199700 

199800 

199900 

200000 

200100 

200200 

200300 

200400 

200S00 

200600 

200700 

200000 

200900 

201000 

201100 

201200 

201300 

201400 

201500 

201600 


6B-E/.-1110a3V 


Table  8.  Concluded 


00 

A 


END 

C*********************************************- ••••••••••••••••••••■«••• 

BLOCKOATk 

COMMON  /RA/XllOltZ(10ltTH(|O)<H(lO)»KllO>tRllO)fBll0l*M<lOI« 

1 CURVE  (10)  tOIRdO)  tXl  IIO)  tX2<10)  tZl  (101  iZZ  (101  (Tl  1101*72(101. 

2 KSl(10)*ET«(10l*RHO(10)*Xa(lll*ZQ(ll)*Ta(llf2i*KRV(3t*CL(2>* 

3 KV(2)tV(IOt3t3I.SL(3)fRHOSL(10t3l*CRVSL(10.3l 

«iVB(10*3*2*ai  .CRVB  (10.8>  .XMAPAXdO)  tZMAPAXdO)  tXSIOdO)  .ZSlDdOl  t 
5 BLB(a) 

COMMON  /lA/ICASE (101 1 JCASE (30) .NCASE (Tl  tNP(3l tNS(3l  .NB(B) 

COMMON  /RS/KI*MItKJ*MJ*NIJ*LINE*BLNKiRINFtERRMAXtOPON*CLLX* 

1 VUCL*ZERO.FORE.AFTE*TOPS*aOTS*OPUPt06IV*P101BO.XJ.ZJ.Xl*ZI*HI. 

2 RI*Bl*THItaJtAl*A2.A3*BlJtTHJ.S*THIJ*$IJ*AI*HJ*RJ.A4*0X 
3fTrPEtASVN*UPFO»UPAFtUPTO*UPBOtDNFO*ONAF*ONTO*DNBO 

COMMON  /IS/IR* IV* IF* IP* IPRINT.NSCG.NSEOPl * ISTOP* JSTOP.KSTOP* 

1 I*J*t*IPl*LMl*Il*NPTS*12*ISL*LB*N*NMl*I3 
COMMON  /LA/SLOPE ( 1 0 ) * lOtN ( 80  > *C ARD ( BO ) 

COMMON  /LS/TI.T 

REAL  K*N*KSI*KI*MI.KU*MJ.MIU*LINE*KV*KRV 
LOGICAL*!  SL0PE*1DEN*CAH0*TI*T 

Data  krv/'line>*'Ogiv>*>ellxi/*line/«line'/*cl/'Zucl'*'Zlcl'/* 

1 KV/iKV5  •••KVO  ‘/.SLOPE/lO^tFl/.T/iT*/ 

DATA  RHO.Xl *X2*Z1 *Z2* Tl *T2*K5I *ETA/90*1 .ETO/ 

DATA  X*Z*  R*H*K*N*B*TH/  80*1. ETO/ 

DATA  CURVE/10**  '/.BLNK/*  '/.RINF/l.EAO/* 

1 ERRMAX/1.E-5/*DIR/10*'OPON'/*OPDN/iOPDN'/*ELLX/iELLX«/* 

2 VUCL/'YUCL'/*ZERO/0./*FORE/  'FORL'/.AFTE/'AFTE'/.TOPS/'TOPS*/* 

3 BOTS/'BOTS</*OPDN/'OPON'/*OGIV/iOOIV*/*OPUP/*OPUP'/ 

DATA  xa*ZO*TO/44*l.E70/ 

DATA  V/  90*1.ET0/*SL/'ZUSL'*'ZLSL'*'VSSL'/ 

DATA  RHOSL/30*1.ETO/*CHVSL/30*>LINE>/ 

DATA  VB/480*1.ET0/*CRVB  /80*'L1NE'/*TVPE/'SYMM*/*ASVM/*A5YM*/ 

DATA  BLB/'ZUCL>*<ZLCL'**YSID'*«ZUSL>*>ZLSL**'YSSL'*'ZSID'*'ZMAP'/ 
DATA  XNAPAX*ZM*PAX*XSID*ZSIO/40*1.ETO/ 

DATA  UPFO  * UPAF  * UPTO  * UPBO  * DNFO  * DNAF  * DNTO  * DNBO 
I /'UPFO'* 'UPAF'* 'UPTO'* 'UPBO*. •ONFO'**DNAFi*iDNTO***DNBOi/ 

DATA  JCASE/582*TT4*834*&26*624*006*010*012*136*590* 

1 626*828*632*014*138* 140*142*638*832*262* 

2 838* 

-9*0/ 

DATA  NCASE/  1*  2*  3*  4*  6* 

1 7*  8 / 

END 


201700 

201800 

201900 

202000 

202100 

202200 

202300 

202400 

202SOO 

202600 

202700 

202B00 

202900 

203000 

203100 

203200 

203300 

203400 

203500 

203600 

203700 

203000 

203900 

204000 

204100 

204200 

204300 

204400 

204500 

204600 

204700 

204000 

204900 

205000 

20S100 

205200 

205300 

205400 

205500 

205600 

205700 

205800 

205900 
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Subroutine 

Name 

Main  Program 

INIT 

INPT 

CASE 

SEGE 

INTR 

QARY 

PRNT 

R((>FL 

QUCK 

LIST 

SLCE 

RFSL 

C(()PY 


Table  9.  Subroutine  Description 


Function 


Overall  program  control 

Initializes  certain  program  variables 

Reads  and  prints  input  data 

Determines  the  case  number  ICASE{I)  of  the 
input  option  for  each  segment  (see  Table  2) 

Determines  the  equation  of  each  segment  and 
computes  certain  intersection  information 

Determines  intersection  points  of  succes- 
sive segments  for  which  only  segment  equa- 
tions are  known  to  that  point  in  code 

Assembles  coordinate  and  slope  data  in  con- 
venient arrays  for  building  QUICK  input 
deck  (QUCK) 

Prints  the  QUICK  input  arrays  assembled  by 
QARY 

Determines  the  equations  of  rounds  and 
fillets  and  computes  intersections  with 
body  line  segments,  prints  results 

Builds  the  QUICK  input  deck  for  the 
symmetrical  case  with  asymmetrical  slices 

Prints  the  QUICK  input  deck  (if  unit  IF  is 
not  the  card  punch)  for  checking  by  user 

Performs  geometry  computation  for  slicing 
planes,  prints  results 

Determines  the  equations  of  rounds  and 
fillets  and  computes  intersections  with 
slicing  plane  segments,  prints  results 

Assembles  coordinate  and  slope  data  in  con- 
venient arrays  for  building  QUICK  input 
deck  for  asymmetrical  case  (QK3D) 
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BLCK 

QK3D 

BL(|)CKDATA 


Table  9.  Concluded 

Computes  default  values  of  certain  body 
line  coordinates  and  checks  validity  of 
previously  computed  coordinates  for 
asymmetrical  case 

Builds  the  QUICK  input  deck  for  the 
asymmetrical  case  with  slices 

Initializes  certain  program  constants 
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Table  10.  Variable  Description 


Variable 

Name 

Type* 

Description 

B 

RA 

z intercept  of  straight  body  line 
segment 

H 

RA 

X coordinate  of  ogive  center 

I 

IS 

Intersection  subscript 

J 

IS 

Scratch  integer 

K 

RA 

Z coordinate  of  ogive  center 

L 

IS 

Scratch  integer 

M 

RA 

Slope  of  body  line  segment 

N 

IS 

Scratch  integer 

R 

RA 

Radius  of  ogival  segment 

S 

RS 

Variable  used  in  selection  from  multi- 
ple intersections,  S=±1. 

T 

LS 

= 'T' 

V(I,J,L) 

RA 

Slice  geometry  array:  J=1,  x coordi- 
nate; J=2,  Z (or  Y)  coordinate;  J=3, 
slope  angle  (deg);  L=1,  upper  slice; 

L=2,  lower  slice;  L=3,  side  slice 

X 

RA 

Input  array  for  X coordinate  of  body 
lines 

z 

RA 

Input  array  for  Z (or  Y)  coordinate  of 
body  lines 

Ai 

RS 

Scratch  variables  i=1,2,... 

BI 

RS 

Z intercept  of  Ith  segment 

BJ 

RS 

Z intercept  of  I- 1st  segment 

CL 

RA 

Body  line  names  used  in  symmetrical  case 
to  assemble  QUICK  input  deck 
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Table  10.  Continued 


DX 

RS 

Round-off  compensation  variable 

HI 

RS 

X coordinate  of  ogive  center  for  Ith 
segment 

HJ 

RS 

X coordinate  of  ogive  center  for  I- 1st 
segment 

IF 

IS 

Unit  number  (data  set  reference  number) 
of  QUICK  file  written  by  QUCK  and  QK3D 

IP 

IS 

Unit  number  of  facility  card  punch 

IR 

IS 

Unit  number  for  input  data 

IW 

IS 

Unit  number  for  facility  printer 

li 

IS 

Scratch  integers,  i=1,2,... 

KI 

RS 

Z coordinate  of  ogive  enter  for  Ith 
segment 

KJ 

RS 

Z coordinate  of  ogive  center  for  I- 1st 
segment 

KV 

RA 

Segment  type  (line  or  ogive) 

LB 

IS 

Overall  D())  variable  tor  asymmetrical 
case,  LB=1,  upper  centerline;  LB=2, 
lower  centerline;  LB=3,  side  (see  Fig.  7) 

MI 

RS 

Slope  of  Ith  segment  (straight  line) 

MJ 

RS 

Slope  of  I- 1st  segment  (straight  line) 

NB{LB) 

lA 

Number  of  segments  in  the  body  line  seg- 
ments in  asymmetrical  case 

NP 

lA 

NS+1 

NS  (L) 

lA 

Number  of  segments  in  the  three  slicing 
planes,  L defined  as  for  V above 

RI 

RS 

Radius  of  ogive  for  Ith  segment 

RJ 

RS 

Radius  of  ogive  for  I- 1st  segment 
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Table  10.  Continued 


SL 

RA 

Body  line  names  for  slicing  planes 

TH 

RA 

Input  array  for  local  slopes  of  body 
lines  at  end  of  segment 

TI 

LS 

=SL(|)PE(I) 

TQ(I,J) 

RA 

Local  segment  slope  (degrees):  J=1 , at 
beginning  of  Ith  segment;  J=2  at  end  of 
Ith  segment;  input  to  subroutine  QUCK 

Ti 

RA 

Slope  (degrees)  of  fillet  or  round:  i=1 
beginning  of  fillet;  i=2,  end  of  fillet 

VB(I,J,N,L) 

RA 

Coordinates  and  slopes  of  the  eight  body 
lines  of  the  asymmetrical  case: 

J= 1 , X coordinate 

J=2,  Z coordinate 

J=3,  slope  (deg);  N=1,  beginning  of 
Ith  segment,  N=2  end  of  Ith  segment 
(N=2  used  only  with  J=3) 

L=1  Upper  centerline  (ZUCL) 

Ij=2  Lower  centerline  (ZLCL) 

L=3  Side  (YSID) 

L=4  Upper  slicing  plane  (ZUSL) 

L=5  Lower  slicing  plane  (ZLSL) 

L=6  Side  slicing  plane  (YSSL) 

L=7  Side  (ZSID) 

L=8  Map  axis  (ZMAPAX) 

(see  Fig.  8) 


XI 

RS 

X coordinate  of  end  of  Ith  segment 

XJ 

RS 

X coordinate  of  end  of  I- 1st  segment 

XQ 

RA 

X coordinate  of  beginning  of  Ith  segment 
(I=1,NSEG)  and  end  of  last  segment 
(I=NSEGP1),  input  to  QUCK 

Xi 

RA 

X coordinate  at  intersections  of  fillet 
or  round  with  body  lines,  i as  above  in  ' 

Zl 

RS 

Z coordinate  at  end  of  Ith  segment 

ZJ 

RS 

Z coordinate  at  end  of  I-lst  segment 
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Table  10.  Continued 


ZQ 

RA 

Z coordinate  of  beginning  of  Ith  segment 
(I=1,NSEG)  and  end  of  last  segment 
(I=NSEGP1),  input  to  QUCK 

Zi 

RA 

Z coordinate  at  intersections  of  fillet 
or  round  with  body  lines,  i as  above  in  Ti 

BIJ 

RS 

Z intercept 

BLB(L) 

RA 

Body  line  names  for  QK3D,  L as  above  for  V] 

DIR 

RA 

Selection  indicator  variable  for  multiple 
root  intersections 

ETA 

RA 

Z coordinate  of  round  or  fillet  center  at 
end  of  Ith  segment 

IP1 

IS 

Scratch  integer,  =1+1 

ISL 

IS 

Slice  indicator:  ISI,=0  no  slices 

ISL=1  slices  present 

KRV 

RA 

Segment  type,  line  or  ogive 

KSI 

RA 

X coordinate  of  round  or  fillet  center  at 
end  of  Ith  segment 

LM1 

IS 

scratch  integer,  =L-1 

MU 

RS 

Tangent  of  straight  line  slope  angle 

NM1 

IS 

Scratch  integer,  =N-1 

RH(|> 

RA  ' 

Input  variable  for  radius  of  round  or 
fillet  at  end  of  Ith  segment 

THI 

RS 

Slope  (deg)  at  end  of  Ith  segment 

THJ 

RS 

Slope  (deg)  at  end  of  I-lst  segment 

AFTE 

RA 

='AFTE',  possible  DIR  value  (aft  solution) 

ASYM 

RA 

='ASYM',  possible  TYPE  value 

BLNK 

RS 

= 'bbbb'  (four  blank  characters) 

B(|)TS 

RS 

= 'B((iTS',  possible  DIR  value  (bottom 
solution) 
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Table  10.  Continued 


CARD (80) 

LA 

Card  image  array 

CRVB(I,L) 

RA 

Curve  type  indicator  for  QUICK,  CRVB 
='LINE'  or  'ELLX’ 

ELLX 

RS 

='ELLX' 

F())RE 

RS 

= 'F(|>RE',  possible  DIR  value  (fore 
solution) 

IDEN(80) 

LA 

Input  variable  for  identification  data 

LINE 

RS 

='LINE',  possible  CURVE  value 

NPTS 

IS 

Number  of  stations  at  which  QUICK  will 
print  out  geometry  data  if  IPRINT/0 

NSEG 

IS 

Number  of  segments  for  symmetrical  case 

(j)GIV 

RS 

= '4)GIV',  possible  CURVE  value 

<|)PDN 

RS 

='^PDN',  possible  DIR  value  (open  down) 

(J)PUP 

RS 

= '<()PUP',  possible  DIR  value,  (open  up) 

RINF 

RS 

Real  infinity,  = 1.E69 

THU 

RS 

Straight-line  slope  (deg) 

TLSL  ] 
TSSL 
TUSL  J 

1 

RA 

Input  arrays  for  local  slopes  of  the., 
lower  slice,  side  slice,  and  upper  slice 
segments,  respectively 

XLSL  ] 
XSSL 
XUSL  J 

1 

1 

RA 

Input  arrays  of  X coordinate  of  end  of 
lower,  side,  and  upper  slicing  plane 
segments,  respectively 

ZLSL  ] 
YSSL 
ZUSL  J 

1 

1 

RA 

Input  arrays  of  Z coordinate  of  end  of 
lower,  side,  and  upper  slicing  plane 
segments,  respectively 

T({)PS 

RS 

= 'T()>PS',  possible  DIR  value  (top  solution) 

TYPE 

RS 

Type  body;  TYPE  =’SYMM',  symmetrical 

='ASYM',  assymetrical - 

XSID 

RA 

X coordinate  of  side  body  line  (Z) 
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Table  10.  Continued 

YUCL 

RS 

’YUCL' 

ZER(j) 

RS 

Real  zero  = n.O 

ZSID 

RA 

Z coordinate  of  side  body  line 

CRVSL{I,L) 

RA 

Curve  type  (='LINE'  or  'ELLX')  for 
slicing  planes,  L as  in  V above 

CURVE 

RA 

Curve  type  for  Ith  segment,  ='LINE'  or 
' (|)GIV’ 

I CASE 

lA 

Integer  for  identifying  input  option 
(see  Table  2) 

INPUT 

Namelist  name  for  input  data 

IST4)P 

IS 

Error  indicator  to  halt  execution  when 
set  different  from  zero 

JCASE 

lA 

Integer  array  used  to  determine  ICASE 
(see  Table  2) 

JST(|)P 

IS 

Error  indicator  similar  to  IST(t)P 

KST4)P 

IS 

Error  indicator  similar  to  ISTcfiP 

NCASE 

lA 

ICASE  values  which  input  options  may 
not  follow  ICASE  = 5 or  9 

RH<t)LS 

RA 

Input  array  for  lower  slice  fillet/round 
radii 

RH(|)SL(I,L) 

RA 

Array  of  fillet/round  radii,  L as  in  V 
above 

RH4)SS 

RA 

Input  array  for  side  slice  fillet/round 
radii 

RH(t)US 

RA 

Input  array  for  upper  slice  fillet/round 
radii 

SUPE 

LA 

Variable  to  indicate  Ith  segment  is  to 
be  tangent  to  end  of  1-1 st  seqment 

ERRMAX 

RS 

Computational  error  limit,  = 5.E-5 
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Table  10.  Concluded 


IPRINT 

IS 

Print  indicator;  IPRINiyo 
KWIKN(}>SE  to  generate  cards 
exercise  the  completed  math 

causes 

to  make  QUICK 
model 

NSEGPI 

IS 

=NSEG+1 

PI(ti180 

RS 

Tr/180 

XHAPAX 

RA 

X coordinate  of  map  axis  on 

input 

ZMAPAX 

RA 

Z coordinate  of  map  axis  on 

input 

*R=real, 

I=integer, 

L= logical,  A=array,  S=scalar 
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